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HUMAN METABOLIC MODELS AND METHODS 
BACKGROUND OF THE INVENTION 

This invention relates generally to analysis 
of the activity of chemical reaction networks and, more 
5 specif ically, to computational methods for simulating 
and predicting the activity of Homo sapiens reaction 
networks. 

Therapeutic agents, including drugs and 
gene-based agents, are being rapidly developed by the 

10 pharmaceutical industry with the goal of preventing or 
treating human disease. Dietary supplements, including 
herbal products, vitamins and amino acids, are also 
being developed and marketed by the nutraceutical 
industry. Because , of the complexity of the biochemical 

15 reaction networks in and between human cells, even 

relatively minor perturbations caused by a therapeutic 
agent or a dietary component in the abundance or 
activity of a particular target, such as a metabolite, 
gene or protein, can affect hundreds of biochemical 

20 'reactions. These perturbations can lead to desirable 
therapeutic effects, "sucK as cell stasis or cell Heath 
in the case of cancer cells or other pathologically 
hyperproliferative cells. However, these perturbations 
can also lead to undesirable side effects, such as 

25 production of toxic byproducts, if the systemic effects 
of the perturbations are not taken into account. 

Current approaches to drug and nutraceutical 
development do not take into account the effect of a 
perturbation in a molecular target on systemic cellular 
30 behavior. In order to design effective methods of 
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repairing, engineering or disabling cellular 
activities, it is essential to understand human 
cellular behavior from an integrated perspective. 

Cellular metabolism, which is an example of a 
5 process involving a highly integrated network of 
biochemical reactions, is fundamental to all normal 
cellular or physiological processes, including 
homeostatis, proliferation, differentiation, programmed 
cell death (apoptosis) and motility. Alterations in 
10 cellular metabolism characterize a vast number of human 
diseases. For example, tissue injury is often 
characterized by increased catabolism of glucose, fatty 
acids and amino acids, which, if persistent, can lead 
to organ dysfunction. Conditions of low oxygen supply 
15 (hypoxia) and nutrient supply, such as occur in solid 
tumors, result in a myriad of adaptive metabolic 
changes including activation of glycolysis and 
neovascularization. Metabolic dysfunctions also 
contribute to neurodegenerative diseases, 
20 cardiovascular disease, neuromuscular diseases, obesity 
and diabetes. Currently, despite the importance of 
cellular metabolism to normal and pathological 
processes, a detailed systemic understanding of 
cellular metabolism in human cells is currently 
25 lacking. 

Thus, there exists a need for models that 
describe Homo -sapiens reaction networks, including core 
metabolic reaction networks and metabolic reaction 
networks in specialized cell types, which can be used 
30 to simulate different aspects of human cellular 
behavior under physiological, pathological and 
• therapeutic conditions. The present invention 
satisfies this need, and provides related advantages as 
well. 
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SUMMARY OF THE INVENTION 

The invention provides a computer readable 
medium or media, including: (a) a data structure 
relating a plurality of Homo sapiens reactants to a 
5 plurality of Homo sapiens reactions, wherein each of 
the Homo sapiens reactions includes a reactant 
identified as a substrate of the reaction, a reactant 
identified as a product of the reaction and a 
stoichiometric coefficient relating the substrate and 

10 the product, (b) a constraint set for the plurality of 
Homo sapiens reactions, and (c) commands for 
determining at least one flux distribution that 
minimizes or maximizes an objective function when the 
constraint set is applied to the data representation, 

15 wherein the at least one flux distribution is 

predictive of a Homo sapiens physiological function. 
In one embodiment, at least one of the Homo sapiens 
reactions in the data structure is annotated to 
indicate an associated gene and the computer readable 

20 medium or media further includes a gene database 

including information characterizing the associated 
gene. In another embodiment, at least one of the Homo 
sapiens reactions is a regulated reaction and the 
computer readable medium or media further includes a 

25 constraint set for the plurality of Homo sapiens 
reactions, wherein the constraint set includes a 
variable constraint for the regulated reaction. 

The invention provides a method for 
predicting a Homo sapiens physiological function, 
30 including: (a) providing a data structure relating a 
plurality of Homo sapiens reactants to a plurality of 
Homo sapiens reactions, wherein each of the Homo 
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sapiens reactions includes a reactant identified as a 
substrate of the reaction, a reactant identified as a 
product of the reaction and a stoichiometric 
coefficient relating the substrate and the product; (b) 
5 providing a constraint set for the plurality of Homo 
sapiens reactions; (c) providing an objective 
function, and (d) determining at least one flux 
distribution that minimizes or maximizes the objective 
function when the constraint set is applied to the data 

10 structure, thereby predicting a Homo sapiens 

physiological function. In one embodiment , at least 
one of the Homo sapiens reactions in the data structure 
is annotated to indicate an associated gene and the 
method predicts a Homo sapiens physiological function 

15 related to the gene. 

The invention provides a method for 
predicting a Homo sapiens physiological function, 
including: (a) providing a data structure relating a 
plurality of Homo sapiens reactants to a plurality of 

20 Homo sapiens reactions, wherein each of the Homo 

sapiens reactions includes a reactant identified as a 
substrate of the reaction, a reactant identified as a 
product of the reaction and a stoichiometric 
coefficient relating the substrate and the product, 

25 wherein at least one of the Homo sapiens reactions is a 
regulated reaction; (b) providing a constraint set for 
the plurality of Homo sapiens reactions, wherein the 
constraint set includes a variable constraint for the 
regulated reaction; (c) providing a condition-dependent 

30 value to the variable constraint; (d) providing an 
objective function, and (e) determining at least one 
flux distribution that minimizes or maximizes the 
objective function when the constraint set is applied 
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to the data structure, thereby predicting a Homo 
sapiens physiological function. 

Also provided by the invention is a method 
for making a data structure relating a plurality of 
5 Homo sapiens reactants to a plurality of Homo sapiens 
reactions in a computer readable medium or media, 
including: (a) identifying a plurality of Homo sapiens 
reactions and a plurality of Homo sapiens reactants 
that are substrates and products of the Homo sapiens 

10 reactions; (b) relating the plurality of Homo sapiens 
reactants to the plurality of Homo. sapiens reactions in 
a data structure, wherein each of the Homo sapiens 
reactions includes a reactant identified as a substrate 
of the reaction, a reactant identified as a product of 

15 the reaction and a stoichiometric coefficient relating 
the substrate and the product; (c) determining a 
constraint set, for the plurality of Homo sapiens 
reactions; (d) providing an objective function; (e) 
determining at least one flux distribution that 

20 minimizes or maximizes the objective function when the 
constraint set is applied to the data structure, and 
(f) if the at least one flux distribution is not 
predictive of a Homo sapiens physiological function, 
then adding a reaction to or deleting a reaction from 

25 the data structure and repeating step (e) , if the at 
least one flux distribution is predictive of a Homo 
sapiens physiological function, then storing the data 
structure in a computer readable medium or media. The 
invention further provides a data structure relating a 

30 plurality of Homo sapiens reactants to a plurality of 
Homo sapiens reactions, wherein the data structure is 
produced by the method. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure -1 shows a schematic representation of 
a hypothetical metabolic network. 

Figure 2 shows mass balance constraints and 
5 flux constraints (reversibility constraints) that can' 
be placed on the hypothetical metabolic network shown 
in Figure 1. 

Figure 3 shows the stoichiometric matrix (S) 
for the hypothetical metabolic network shown in Figure 
10 1. 

Figure 4 shows , in Panel A, an exemplary 
biochemical reaction network and in Panel B, an 
exemplary regulatory control structure for the reaction 
network in panel A. 

15 DETAILED DESCRIPTION OF THE INVENTION 

The present invention provides in silico 
models that describe the interconnections between genes 
in the Homo sapiens genome and their associated 
reactions and reactants. The models can be used to 

20 simulate different aspects of the cellular behavior of 
human cells under different normal, pathological and 
therapeutic conditions, thereby providing valuable 
information for therapeutic, diagnostic and research 
applications. An advantage of the models of the 

25 invention is that they provide a holistic approach to 
simulating and predicting the activity of Homo sapiens 
cells. The models and methods can also be extended to 
simulate the activity of multiple interacting cells, 
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including organs, physiological systems and whole body 
metabolism. 

As an example , the Homo sapiens metabolic 
models of the invention can be used to determine the 
5 effects of changes from aerobic to anaerobic 

conditions, such as occurs in skeletal muscles during 
exercise or in tumors, or to determine the effect of 
various dietary changes. The Homo sapiens metabolic 
models can also be used to determine the consequences 
10 of genetic defects, such as deficiencies in metabolic 
enzymes such as phosphof ructokinase, phosphoglycerate 
kinase, phosphoglycerate mutase, lactate dehydrogenase 
and adenosine deaminase. 

The Homo sapiens metabolic models can also be 
15 used to choose appropriate targets for drug design. 
Such targets include genes, proteins or reactants, 
which when modulated positively or negatively in a 
simulation produce a desired therapeutic result. The 
models and methods of the invention can also be used to 
20 predict the effects of a therapeutic agent or dietary 
supplement on a cellular function of interest. 
Likewise, the models and methods can be used to predict 
both desirable and undesirable side effects of the 
therapeutic agent on an interrelated cellular function 
25 in the target cell, as well as the desirable and 

undesirable effects that may occur in other cell types. 
Thus, the models and methods of the invention can make 
the drug development process more rapid and cost 
effective than is currently possible. 

30 - The Homo sapiens metabolic models can also be 

used to predict or validate the assignment of 
particular biochemical reactions to the enzyme-encoding 
genes found in the genome, and to identify the presence 
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of reactions or pathways not indicated by current 
genomic data. Thus, the models can be used to guide 
the research and discovery process, potentially leading 
to the identification of new enzymes, medicines or 
5 metabolites of clinical importance. 

The models of the invention are based on a 
data structure relating a plurality of Homo sapiens 
reactants to a plurality of Homo sapiens reactions, 
wherein each of the Homo sapiens reactions includes a 

10 reactant identified as a substrate of the reaction, a 
reactant identified as a product of the reaction and a 
stoichiometric coefficient relating the .substrate and 
the product. The reactions included in the data 
structure can be those that are common to all or most 

15 Homo sapiens cells, such as core metabolic reactions, 
or reactions specific for one or more given cell type. 

As used herein, the term "Homo sapiens 
reaction" is intended to mean a conversion that 
consumes a substrate or forms a product that occurs in 

20 or by a Homo sapiens cell. The term can include a 
conversion that occurs due to the activity of one or 
more enzymes that are genetically encoded by a Homo 
sapiens genome. The term can also include a conversion 
that occurs spontaneously in a Homo sapiens cell. 

25 Conversions included in the term include, for example, 
changes in chemical composition such as those due to 
nucleophilic or electrophilic addition, nucleophilic or 
electrophilic substitution, elimination, isomerization, 
deamination, phosphorylation, methylation, reduction, 

30 oxidation or changes in location such as those that 

occur due to a transport reaction that moves a reactant 
from one cellular compartment to another. In the case 
of a transport reaction, the substrate and product of 
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the reaction can be chemically the same and the 
substrate and product can be differentiated according 
to location in a particular cellular compartment. 
Thus, a reaction that transports a chemically unchanged 
reactant from a first compartment to a second 
compartment has as its substrate the reactant in the 
first compartment and as its- product the reactant in 
the second compartment. It will be understood that 
when used in reference to an in siliao model or data 
structure, a reaction is intended to be a 
representation of a chemical conversion that consumes a 
substrate or produces a product. 

As used herein, the term "Homo sapiens 
reactant" is intended to mean a chemical that is a 
15 substrate or a product of a reaction that occurs in or 
by a Homo sapiens cell. The term can include 
substrates or products of reactions performed by one or 
more enzymes encoded by a Homo sapiens genome, 
reactions occurring in Homo sapiens that are performed 
20 by one or more non-genet ically encoded macromolecule, 
protein or enzyme, or reactions that occur 
spontaneously in a Homo sapiens cell. Metabolites are 
understood to be reactants within the meaning of the 
term. It will be understood that when used in 
25 reference to an in silico model or" data structure , a 
reactant is intended to be a representation of a 
chemical that is a substrate or a product of a reaction 
that occurs in or by a Homo sapiens cell. 



5 



10 



As used herein the term "substrate" is 
30 intended to mean a reactant that can be converted to 
one or more products by a reaction. The term can 
include, for example,, a reactant that is to be 
chemically changed due to nucleophilic or electrophilic 
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addition, nucleophilic or electrophilic substitution, 
elimination, isomerization, deamination, 
phosphorylation, methylation, reduction, oxidation or 
that is to change location such as by being transported 
5 across a membrane or to a different compartment- 

As used herein, the term "product" is 
intended to mean a reactant that results from a 
reaction with one or more substrates. The term can 
include, for example, a reactant that has been 

10 chemically changed due to nucleophilic or electrophilic 
addition, nucleophilic or electrophilic substitution, 
elimination, isomerization, deamination, 
phosphorylation, methylation, reduction or oxidation or 
that has changed location such as by being transported 

15 across a membrane or to a different compartment. 

As used herein, the term "stoichiometric 
coefficient" is intended to mean a numerical constant 
correlating the number of one or more reactants and the 
number of one or more products in a chemical reaction. 

20 Typically, the numbers are integers as they denote the 
number of molecules of each reactant in an elementally 
balanced chemical equation that describes the 
corresponding conversion. However, in some cases the 
numbers can take on non-integer values, for example, 

25 when used in a lumped reaction or to reflect empirical 
data. 

As used herein, the term "plurality, " when 
used in reference to Homo sapiens reactions or 
reactants, is intended to mean at least 2 reactions or 
30 reactants. The term can include any number of Homo 
sapiens reactions or reactants in the range from 2 to 
the number of naturally occurring reactants or 
reactions for a particular of Homo sapiens cell. Thus, 
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the term can include,, for example, at least 10, 20, 30, 
50, 100, 150, 200, 300, 400, 500, 600 or more reactions 
or reactants. The number of reactions or reactants can 
be expressed as a portion of the total number of 
5 naturally occurring reactions for a particular Homo 

sapiens cell, such as at least 20%, 30%, 50%, 60%, 75%, 
90%, 95% or 98% of the total number of naturally 
occurring reactions that occur in a particular Homo 
sapiens cell. 

10 As used herein, the term "data structure" is 

intended to mean a physical or logical relationship 
among data elements, designed to support specific data 
manipulation functions. The term can include, for 
example, a list of data elements that can be added 

15 combined or otherwise manipulated such as a list of 

representations for reactions from which reactants can 
be related in a matrix or network. The term can also 
include a matrix that correlates data elements from two 
or more lists of information such as a matrix that 

20 correlates reactants to reactions. Information 

included in the term can represent j for example, a 
substrate, or product of a chemical reaction, a chemical 
reaction relating one or more substrates to one or more 
products, a constraint placed on a reaction, or a 

25 stoichiometric coefficient. 



As used herein, the term "constraint" is 
intended to mean an upper or lower boundary for a 
reaction. A boundary can specify a minimum or maximum 
flow of mass, electrons or energy through a reaction. 
30 A boundary can further specify directionality of a 

reaction. A boundary can be a constant value such as 
zero, infinity, or a numerical value such as an 
integer. Alternatively, a boundary can be a variable 
boundary value as set forth below. 
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As used herein, the term "variable," when 
used in reference to a constraint is intended to mean 
capable of assuming any of a set of values in response 
to being acted upon by a constraint function. The term 

5 "function, " when used .in the context of a constraint, 
is intended to be consistent with the meaning of the 
term as it is understood in the computer and 
mathematical arts. A function can be binary such that 
changes correspond to a reaction being off or on. 

10 Alternatively, continuous functions can be used such 

that changes in boundary values correspond to increases 
or decreases in activity. Such increases or decreases 
can also be binned or effectively digitized by a 
function capable of converting sets of values to 

15 discreet integer values. A function included in the 
term can correlate a boundary value with the presence, 
absence or amount of a biochemical reaction network 
participant such as a reactant, reaction, enzyme or 
gene. A function included in the term can correlate a 

20 boundary value with an outcome of at least one reaction 
in a reaction network that includes the reaction that 
is constrained by the boundary limit. A function 
included in the term can also correlate a boundary 
value with an environmental condition such as time, pH, 

25 temperature or redox potential. 

' As used herein, the term "activity, " when 
used in reference to a reaction, is intended to mean 
the amount of product produced by the reaction, the 
amount of substrate consumed- by the reaction or the 

30 rate at which a product is produced or a substrate is 
consumed. The amount of product produced by the 
reaction, the amount of substrate consumed by the 
reaction or the rate at which a product is produced or 
a substrate is consumed can also be referred to as the 

35 flux for the reaction. 
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As used herein, the term "activity," when 
used in reference to a Homo sapiens cell, is intended 
to mean the magnitude or rate of a change from an 
initial state to a final state. The term can include, 
5 for example, the amount of a chemical consumed or 
produced by a cell, the rate at which a chemical is 
consumed or produced by a cell,* the amount or rate of 
growth of a cell or the amount of or rate at which 
energy, mass or electrons flow through a particular 
10 subset of reactions. 



The invention provides a computer readable 
medium, having a data structure relating a plurality of 
Homo sapiens reactants to a plurality of Homo sapiens 
reactions, wherein each of the Homo sapiens reactions 
15 includes a reactant identified as a substrate of the 
reaction, a reactant identified as a product of the 
reaction and a stoichiometric coefficient relating the 
substrate and the product. 



Depending on the application, the plurality 
20 of Homo sapiens reactions can include reactions 

selected from core metabolic reactions or peripheral 
metabolic reactions. As used herein, the term "core," 
when used in reference to a metabolic pathway, is 
intended to mean a metabolic pathway selected from 
25 glycolysis/gluconeogenesis, the pentose phosphate 
pathway (PPP) , the tricarboxylic acid (TCA) cycle, 
glycogen storage, electron transfer system (ETS) , the 
malate/aspartate shuttle, the glycerol phosphate 
shuttle, and plasma and mitochondrial membrane 
30 transporters. As used herein, the term "peripheral," 
when used in reference to a metabolic pathway, is 
intended to mean a metabolic pathway that includes one 
or more reactions that are not a part of a core 
metabolic pathway. 
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A plurality of Homo sapiens reactants can be 
related to a plurality of Homo sapiens reactions in any 
data structure that represents , for each reactant, the 
reactions by which it is consumed or produced. Thus, 
the data structure, which is referred to herein as a 
"reaction network data structure/' serves as a 
representation of a biological reaction network or 
system. An example of a reaction network that can be 
represented in a reaction network data structure of the 
invention is the collection of reactions that 
constitute the core metabolic reactions of Homo 
sapiens, or the metabolic reactions of a skeletal 
muscle cell, as shown in the Examples. 

The choice of reactions to include in a 
particular reaction network .data structure, from among 
all the possible reactions that can occur in human 
cells, depends on the cell type or types and the 
physiological, pathological or therapeutic condition 
being modeled, and can be determined experimentally or 
from the literature, as described further below. 

The reactions to be included in a particular 
network data structure of Homo sapiens can be 
determined experimentally using, for example, gene or 
protein expression profiles, where the molecular 
characteristics of the cell can be correlated to the 
expression levels. The expression or lack of . 
expression of genes or. proteins in a cell type can be 
used in determining whether a reaction is included in 
the model by association to the expressed gene(s) and 
or protein (s). Thus, it is possible to use 
experimental technologies to determine which genes 
and/or proteins are expressed in a specific cell type, 
and to further use this information to. determine which 
reactions are present in the cell type of interest. In 
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this way a subset of reactions from all of those 
reactions that can occur in human cells are selected to 
comprise the set of reactions that represent a specific 
cell type. cDNA expression profiles have been 
5 demonstrated to be useful , for example, for 

classification of breast cancer cells (Sorlie et al., 
Proc. Natl. Acad. Sci. U.S.A. 98 (19) : 10869-10874 
(2001)). 

The methods and models of the invention can 
10 be applied to any Homo sapiens cell type at any stage 
of differentiation, including, for example, embryonic 
stem cells, hematopoietic stem cells, differentiated 
hematopoietic cells, skeletal muscle cells, cardiac 
muscle cells, smooth muscle cells, skin cells, nerve 
15 cells, kidney cells, pulmonary cells, liver cells, 

adipocytes and endocrine cells (e.g. beta islet cells 
of the pancreas, mammary gland cells, adrenal cells, 
and other specialized hormone secreting cells) . 

The methods and models of the invention can 
20 be applied to normal cells or pathological cells. 

Normal cells that exhibit a variety of physiological 
activities of interest, including homeostasis, 
proliferation, differentiation, apoptosis, contraction 
and motility, can be modeled. -Pathological— cells can 
25 also be modeled, including cells that reflect genetic 
or developmental abnormalities, nutritional 
deficiencies, environmental assaults, infection (such 
as by bacteria, viral, protozoan or fungal agents), 
neoplasia, aging, altered immune or endocrine function, 
30 tissue damage, or any combination of these factors. 
The pathological cells can be representative of any 
type of human pathology, including, for example, 
various metabolic disorders of carbohydrate, lipid or 
protein metabolism, obesity, diabetes, cardiovascular 
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disease, fibrosis, various cancers, kidney failure, 
immune pathologies, neurodegenerative diseases, and 
various monogenetic metabolic diseases described in the 
Online Mendelian Inheritance in Man database (Center 
5 for Medical Genetics, Johns Hopkins University 

(Baltimore, MD) and National Center for Biotechnology 
Information, National Library of Medicine (Bethesda, 
MD) ) . 

The methods and models of the invention can 

10 also be applied to cells undergoing therapeutic 

perturbations, such as cells treated with drugs that 
target participants in a reaction network, cells 
treated with gene-based therapeutics that increase 'or 
decrease expression of an encoded protein, and cells 

15 treated with radiation. As used herein, the term 
"drug" refers to a compound of any molecular nature 
with a known or proposed therapeutic function, 
including, for example, small molecule compounds, 
peptides and other macromolecules, peptidomimetics and 

20 antibodies, any of which can optionally be tagged with 
cytostatic, targeting or detectable moieties. The term 
"gene-based therapeutic" refers- -to -nucleic acid 
therapeutics, including, for example, expressible genes 
with normal or altered protein activity, antisense 

25 compounds, ribozymes, DNAzymes, RNA interference 

compounds (RNAi) and the like. The therapeutics can 
target any reaction network participant, in any 
cellular location, including participants in 
extracellular, cell surface, cytoplasmic, mitochondrial 

30 and nuclear locations. Experimental data that are 
gathered on the response of cells to therapeutic 
treatment, such as alterations in gene or protein 
expression profiles, can be used to tailor a network 
for a pathological state of a particular cell type. 
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The methods and models of the invention can 
be applied to Homo sapiens cells as they exist in any 
form, such as in primary cell isolates or in 
established cell lines, or in the whole body, in intact 
5 organs or in tissue explants. Accordingly, the methods 
and models can take into account intercellular 
communications and/or inter-organ communications, the 
effect of adhesion to a substrate or neighboring cells 
(such as a stem cell interacting with mesenchymal cells 
10 or a cancer cell interacting with its tissue 

microenvironment, or beta-islet cells without normal 
stroma) , and other interactions relevant to 
multicellular systems. 

The reactants to be used in a reaction 

15 network data structure of the invention can be obtained 
from or stored in a compound database. As used herein, 
the term "compound database'' is intended to mean a 
computer readable medium or media containing a 
plurality of molecules that includes substrates and 

20 products of biological reactions. The plurality of 
molecules can include molecules found in multiple 
organisms, thereby constituting a universal compound 
database. Alternatively, the plurality of molecules 
can be limited to those that occur in a particular 

25 organism,, thereby constituting an organism-specific 
compound database. Each reactant in a compound 
database can be identified according to the chemical 
species and the cellular compartment in which it is 
present. Thus, for example, a distinction can be made 

30 between glucose in the extracellular compartment versus 
glucose in the cytosol. Additionally each of the 
reactants can be specified as a metabolite of a primary 
or secondary metabolic pathway. Although 
identification of a reactant as a metabolite of a 

35 primary or secondary metabolic pathway does not 
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indicate any chemical distinction between the reactants 
in a reaction, such a designation can assist in visual 
representations of large networks of reactions. 

As used herein, the term M compartment " is 
intended to mean a subdivided region containing at 
least one reactant, such that the reactant is separated 
from at least one other reactant in a second region. A 
subdivided region included in the term can be 
correlated with a subdivided region of a cell. Thus, a 
subdivided region included in the term can be, for 
example, the intracellular space of a cell; the 
extracellular space around a cell; the periplasmic 
space, the interior space of an organelle such as a 
mitochondrium, endoplasmic reticulum, Golgi apparatus, 
vacuole or nucleus; or any subcellular space that is 
separated from another by a membrane or other physical 
barrier. Subdivided regions can also be made in order 
to create virtual boundaries in a reaction network that 
are not correlated with physical barriers. Virtual 
boundaries can be made for the purpose of segmenting 
the reactions in a network into different compartments 
or substructures . 

As used herein, the term ^substructure" is 
intended to mean a portion of the information in a data 
structure that is separated from other information in 
the data structure such that the portion of information 
can be separately manipulated or analyzed. The term 
can include portions subdivided according to a 
biological function including, for example, information 
relevant to a particular metabolic pathway such as an 
internal flux pathway, exchange flux pathway, central 
metabolic pathway, peripheral metabolic pathway, or 
secondary metabolic pathway. The term can include 
portions subdivided according to computational or 
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mathematical principles that allow for a particular 
type of analysis or manipulation of the data structure. 

The reactions included in a reaction network 
data structure can be obtained from a metabolic 
reaction database that includes the substrates, 
products, and stoichiometry of a plurality of metabolic 
reactions of Homo sapiens. The reactants in a reaction 
network data structure can be designated as either 
substrates or products of a particular reaction, each 
with a stoichiometric coefficient assigned to it to 
describe the chemical conversion taking place in the 
reaction. Each reaction is also described as occurring 
in either a reversible or irreversible direction. 
Reversible reactions can either be represented as one 
reaction that operates in both the forward and reverse 
direction or be decomposed into two irreversible 
reactions, one corresponding to the forward reaction 
and the other corresponding to the backward reaction. 

Reactions included in a reaction network data 
structure can include intra^system or exchange 
reactions. Intra-system reactions are the chemically 
and electrically balanced interconversions of chemical 
species and transport processes, which serve to 
replenish or drain the relative amounts of certain 
metabolites. These intra-system reactions can be 
classified as either being transformations or 
translocations. A transformation is a reaction that 
contains distinct sets of compounds as substrates and 
products, while a translocation contains reactants 
located in different compartments. Thus a reaction 
that simply transports a metabolite from the 
extracellular environment to the cytosol, without 
changing its chemical composition is solely classified 
as a translocation, while a reaction that takes an 
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extracellular substrate and converts it into a 
cytosolic product is both a translocation and a 
transformation. 

Exchange reactions are those which constitute 
5 sources and sinks, allowing the passage of metabolites 
into and out of a compartment or across a hypothetical 
system boundary. These reactions are included in a 
model for simulation purposes and represent the 
metabolic demands placed on Homo sapiens. While they 
10 may be chemically balanced in certain cases, they are 
typically not balanced and can often have only a single 
substrate or product. As a matter of convention the 
exchange reactions are further classified into demand 
exchange and input/output exchange reactions. 

15 The metabolic demands placed on the Homo 

sapiens metabolic reaction network can be readily 
determined from the dry weight composition of the cell 
which is available in the published literature or which 
can be determined experimentally. The uptake rates and 

20 maintenance reguirements for Homo sapiens cells can 
also be obtained from the. published literature or 
determined experimentally. . 

Input /output exchange reactions are used to 
allow extracellular reactants to enter or exit the 

25 reaction network represented by a model of the 

invention. For each of the extracellular metabolites a 
corresponding input /output exchange reaction can be 
created. These reactions are always reversible with 
the metabolite indicated as a substrate with a 

30 stoichiometric coefficient of one and no products 

produced by the reaction. This particular convention 
is adopted to allow the reaction to take on a positive 
flux value (activity level) when the metabolite is 
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being produced or removed from the reaction network and 
a negative flux value when the metabolite is being 
consumed or introduced into the reaction network. 
These reactions will be further constrained during the 
5 course of a simulation to specify exactly which 

metabolites are available to the cell and which can be 
excreted by the cell. 

A demand exchange reaction is always 
specified as an irreversible reaction containing at 

10 least one substrate. These reactions are typically 
formulated to represent the production of an 
intracellular metabolite by the metabolic network or 
the aggregate production of many reactants in balanced 
ratios such as in the representation of a reaction that 

15 leads to biomass formation, also referred to as growth . 

A demand exchange reactions can be introduced 
for any metabolite in a model of the invention. Most 
commonly these reactions are introduced for metabolites 
that are required to be produced by the cell for the 

20 purposes of creating a new cell such as amino acids, 
nucleotides, phospholipids, and other biomass 
constituents, or metabolites that are to be produced 
for alternative purposes. Once these metabolites are 
identified, a demand exchange reaction that is 

25 irreversible and specifies the metabolite as a 

substrate with a stoichiometric coefficient of unity 
can be created. With these specifications, if the 
reaction is active it leads to the net production of 
the metabolite by the system meeting potential 

30 production demands. Examples of processes that can be 
represented as a demand exchange reaction in a reaction 
network data structure and analyzed by the methods of 
the invention include, for example, production or 
secretion of an individual protein; production or 
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secretion of an individual metabolite such as an amino 
acid, vitamin, nucleoside, antibiotic or surfactant; 
production of ATP for extraneous energy requiring 
processes such as locomotion; or formation of biomass 
5 constituents. 

In addition to these demand exchange 
reactions that are placed on individual metabolites, 
demand exchange reactions that utilize multiple 
metabolites in defined stoichiometric ratios can be 

10 introduced. These reactions are referred to as 

aggregate demand exchange reactions. An example of an 
aggregate demand reaction is a reaction used to 
simulate the concurrent growth demands or production 
requirements associated with cell growth that are 

15 placed on a cell, for example, by simulating the 
formation of multiple biomass constituents 
simultaneously at a particular cellular growth rate. 

A hypothetical reaction network is provided 
in Figure 1 to exemplify the above -de scribed reactions 

20 and their interactions. The reactions can be 

represented in the exemplary data structure shown in 
Figure 3 as set forth below. The reaction network, 
shown in Figure 1, includes intrasystem reactions that 
occur entirely within the compartment indicated by the 

25 shaded oval such as reversible reaction R 2 which acts 
on reactants B and G and reaction R 3 which converts one 
equivalent of B to 2 equivalents of F. The reaction 
network shown in Figure 1 also contains exchange 
reactions such as input/output exchange reactions A^ 

30 and E xt , and the demand exchange reaction, V grcwth , which 
represents growth in response to the one equivalent of 
D and one equivalent of F. Other intrasystem reactions 
include R a which is a translocation and transformation 
reaction that translocates reactant A into the 
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compartment and transforms it to reactant G and 
reaction R 6 which is a transport reaction that 
translocates reactant E out of the compartment, 

A reaction network can be represented as a 
5 set of linear algebraic equations which can be 

presented as a stoichiometric matrix S, with S being an 
m x n matrix where m corresponds to the number of 
reactants or metabolites and n corresponds to the 
number of reactions taking place in the network. An 

10 example of a stoichiometric matrix representing the 

reaction network of Figure 1 is shown in Figure 3. As 
shown in Figure 3, each column in the matrix 
corresponds to a particular reaction n, each row 
corresponds to a particular reactant m, and each 

15 element corresponds to the stoichiometric coefficient 
of the reactant m in the reaction denoted n. The 
stoichiometric matrix includes intra-system reactions 
such as R 2 and R 3 which are related to reactants that 
participate in the respective reactions according to a 

20 stoichiometric coefficient having a sign indicative of 
whether the reactant is a substrate or product of the 
reaction and a value correlated with the number of 
equivalents of the reactant consumed or produced by the 
reaction. Exchange reactions such as -E xt and -A xt are 

25 similarly correlated with a stoichiometric coefficient. 
As exemplified by reactant E, the same compound can be 
treated separately as an internal reactant (E) and an ' 
external reactant (E external ) such that an exchange 
reaction (R 6 ) exporting the compound is correlated by 

30 stoichiometric coefficients of -1 and 1, respectively. 
However, because the compound is treated as a separate 
reactant by virtue of its compartmental location, a 
reaction, such as R 5 , which produces the internal 
reactant (E) but does not act on the external reactant 

35 (E externa i) is correlated by stoichiometric coefficients 
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of 1 and 0, respectively. Demand reactions such as 
Vgrowth can also be included in the stoichiometric matrix 
being correlated with substrates by an appropriate 
stoichiometric coefficient. 

5 

As set forth in further detail below, a 
stoichiometric matrix provides a convenient format for 
representing and analyzing a reaction network because 
it can be readily manipulated and used to compute 

10 network properties, for example, by using linear 

programming or general convex analysis. A reaction 
network data structure can take on a variety of formats 
so long as it is capable of relating reactants and 
reactions in the manner exemplified above for a 

15 stoichiometric matrix and in a manner that can be 
manipulated to determine an activity of one or more 
reactions using methods such as those exemplified 
below. Other examples of reaction network data 
structures that are useful in. the invention include a 

20 connected graph, list of chemical reactions or a table 
of reaction equations. 



A reaction network data structure can be 
constructed to include all reactions that are involved 

25 in Homo sapiens metabolism or any portion thereof. A 
portion of Homo sapiens metabolic reactions that can be 
included in a reaction network data structure of the 
invention includes, for example, a central metabolic 
pathway such as glycolysis, the TCA cycle, the PPP or 

30 ETS; or a peripheral metabolic pathway such as amino 
acid biosynthesis, amino acid degradation, purine 
biosynthesis, pyrimidine biosynthesis, lipid 
biosynthesis, fatty acid metabolism, vitamin or 
cof actor biosynthesis, transport processes and 

35 alternative carbon source catabolism. Examples of 
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individual pathways within the peripheral pathways are 
set forth in Table 1. 

Depending upon a particular application, a 
reaction network data structure can include a plurality 
5 of Homo sapiens reactions including any or all of the 
reactions listed in Table 1. 

For some applications , it can be advantageous 
to use a reaction network data structure that includes 
a minimal number of reactions to achieve a particular 

10 Homo sapiens activity under a particular set of 

environmental conditions. A reaction network data 
structure having a minimal number of reactions can be 
identified by performing the simulation methods 
described below in an iterative fashion where different 

15 reactions or sets of reactions are systematically 
removed ahd the effects observed. Accordingly, the 
invention provides a computer readable medium, 
containing a data structure relating a plurality of 
Homo sapiens reactants to a plurality of Homo sapiens 

20 reactions, wherein the plurality of Homo sapiens 

reactions contains at least 65 reactions. For example, 
the core metabolic reaction database shown in Tables 2 
and 3 contains 65 reactions, and is sufficient to 
simulate aerobic and anaerobic metabolism on a number 

25 of carbon sources, including glucose. 

Depending upon the particular cell type or 
types, the physiological, pathological or therapeutic 
conditions being tested and the desired activity, a 
reaction network data structure can contain smaller 
30 numbers of reactions such as at least 200, 150, 100 or 
50 reactions. A reaction network data structure having 
relatively few reactions can provide the advantage of 
reducing computation time and resources required to 
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perform a simulation. When desired, a reaction network 
data structure having a particular subset of reactions 
can be made or used in which reactions that are not 
relevant to the particular simulation are omitted. 
5 Alternatively, larger numbers of reactions can be 

included in order to increase the accuracy or molecular 
detail of the methods of the invention or to suit a 
particular application. Thus, a reaction network data 
structure can contain at least 300, 350, 400, 450, 500, 

10 550, 600 or more reactions up to the number of 

reactions that occur in or by Homo sapiens or that are 
desired to simulate the activity of the full set of 
reactions occurring in Homo sapiens. A reaction 
network data structure that is substantially complete 

15 with respect to the metabolic reactions of Homo sapiens 
provides the advantage of being relevant to a wide 
range of conditions to be simulated, whereas those with 
smaller numbers of metabolic reactions are limited to a 
particular subset of conditions to be simulated. 

20 A Homo sapiens reaction network data 

structure can include one or more reactions that occur 
in or by Homo sapiens and that do not occur, either 
naturally or following manipulation, in or by another 
organism, such as Saccharomyces cerevisiae. It is 

25 understood that a Homo sapiens reaction network data 
structure of a particular cell type can also include 
one or more reactions that occur in another cell type. 
Addition of such heterologous reactions to a reaction 
network data structure of the invention can be used in 

30 methods to predict the consequences of heterologous 

gene transfer and protein expression, for example, when 
designing in vivo and ex vivo gene therapy approaches. 
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The reactions included in a reaction network 
data structure of the invention can be metabolic 
reactions. A reaction network data structure can also 
be constructed to include other types of reactions such 
5 as regulatory reactions , signal transduction reactions , 
cell cycle reactions , reactions controlling 
developmental processes, reactions involved in 
apoptosis, reactions involved in responses to hypoxia, 
reactions involved in responses to. cell-cell or cell- 
10 substrate interactions, reactions involved in protein 
synthesis and regulation thereof, reactions involved in 
gene transcription and translation, and regulation 
thereof, and reactions involved in assembly of a cell 
and its subcellular components. 

15 A reaction network data structure or index of 

reactions used in the data structure such as that 
available in a metabolic reaction database, as 
described above, can be annotated to include 
information about a particular reaction. A reaction 

20 can be annotated to indicate, for example, assignment 
of the reaction to a protein, macromolecule or enzyme 
that performs the reaction, assignment of a gene(s) 
that codes for the protein, macromolecule or enzyme, 
the Enzyme Commission (EC) number of the particular 

25 metabolic reaction, a subset of reactions to which the 
reaction belongs, citations to references from which 
information was obtained, or a level of confidence with 
which a reaction is believed to occur in Homo sapiens. 
A computer readable medium or media of the invention 

30 can include a gene database containing annotated 

reactions. Such information can be obtained during the 
course of building a metabolic reaction database or 
model of the invention as described below. 
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As used herein, the term "gene database" is 
intended to mean a computer readable medium or media 
that contains at least one reaction that is annotated 
to assign a reaction to one or more macromolecules that 
5 perform the reaction or to assign one or more nucleic 
acid that encodes the one or more macromolecules that 
perform the reaction. A gene database can contain a 
plurality of reactions, some or all of which are 
annotated. An annotation can include, for example, a 

10 name for a macromolecule; assignment of a function to a 
macromolecule; assignment of an organism that contains 
the macromolecule or produces the macromolecule; 
assignment of a subcellular location for the 
macromolecule; assignment of conditions under which a 

15 macromolecule is regulated with respect to performing a 
reaction, being expressed or being degraded; assignment 
of a cellular component that regulates a macromolecule; 
an amino acid or nucleotide sequence for the 
macromolecule; or any other annotation found for a 

20 macromolecule in a genome database such as those that 
can be found in Genbank, a site maintained by the NCBI 
(ncbi.nlm.gov) , the Kyoto Encyclopedia of Genes and 
Genomes (KEGG) (www.genome.ad.jp/kegg/), the protein 
database SWISS-PR0T (ca.expasy.org/sprot/), the 

25 LocusLink database maintained by the NCBI 

(www.ncbi.nlm.nih.gov/LocusLink/), the Enzyme 
Nomenclature database maintained by G.P. Moss of Queen; 
Mary and Westfield College in the United Kingdom 
(www.chem.qmw.ac.uk/iubmb/enzyme/) . 

30 A gene database of the invention can include 

a substantially complete collection of genes or open 
reading frames in Homo sapiens or a substantially 
complete collection of the macromolecules encoded by 
the Homo sapiens genome. Alternatively, a gene 

35 database can include a portion of genes or open reading 
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frames in Homo sapiens or a portion of the 
macromolecules encoded by the Homo sapiens genome, such 
as the portion that includes substantially all 
metabolic genes or macromolecules. The portion can be 
at least 10%, 15%, 20%, 25%, 50%, 75%, 90% or 95% of 
the genes or open reading frames encoded by the Homo 
sapiens genome, or the macromolecules encoded therein. 
A gene database can also include macromolecules encoded 
by at least a portion of the nucleotide sequence for 
the Homo sapiens genome such as at least 10%, 15%, 20%, 
25%, 50%, 75%, 90% or 95% of the Homo sapiens genome. 
Accordingly, a computer readable medium or media of the 
invention can include at least one reaction for each 
macromolecule encoded by a portion of the Homo sapiens 
genome . 

An in silico Homo sapiens model of the 
invention can be built by an iterative process which 
includes gathering information regarding particular 
reactions to be added to a model, representing the 
reactions in a reaction network data structure, and 
performing preliminary simulations wherein a set of 
constraints is placed on the reaction network and the 
output evaluated to identify errors in the network. 
Errors in the network such as gaps that lead to non- 
natural accumulation or consumption of a particular 
metabolite can be identified as described below and 
simulations repeated until a desired performance of the 
model is attained. An exemplary method for iterative 
model construction is provided in Example I. 

Thus, the invention provides a method for 
making a data structure relating a plurality of Homo 
sapiens reactants to a plurality of Homo sapiens 
reactions in a computer readable medium or media. The 
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method includes the steps of: (a) identifying a 
plurality of Homo sapiens reactions and a plurality of 
Homo sapiens reactants that are substrates and products 
of the Homo sapiens reactions; (b) relating the 
5 plurality of Homo sapiens reactants- to the plurality of 
Homo sapiens reactions in a data structure, wherein 
each of the Homo sapiens reactions includes a reactant 
identified as a substrate of the reaction, a reactant 
identified as a product of the reaction and a 

10 stoichiometric coefficient relating the substrate and 
the product; (c) making a constraint set for the 
plurality of Homo sapiens reactions; (d) providing an 
objective function; (e) determining at least one flux 
distribution that minimizes or maximizes the objective 

15 function when the constraint set is applied to the data 
structure, and (f ) if the at least one flux 
distribution is not predictive of Homo sapiens 
physiology, then adding a reaction to or deleting a 
reaction from the data structure and repeating step 

20 (e) , if the at least one flux distribution is 

predictive of Homo sapiens physiology, then storing the 
data structure in a computer readable medium or media. 

Information to be included in a data 
structure of the invention can be gathered from a 
25 variety of sources including, for example, annotated 

genome sequence information and biochemical literature. 

Sources of annotated human genome sequence 
information include, for example, KEGG, SWISS-PROT, 
LocusLink, the Enzyme Nomenclature database, the 
3Q International Human Genome Sequencing Consortium and 
commercial databases. KEGG contains a broad range of 
information, including a substantial amount of 
metabolic reconstruction. The genomes of 63 organisms 
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can be accessed here, with gene products grouped by 
coordinated functions, often represented by a map 
(e.g., the enzymes involved in glycolysis would be 
grouped together) . The maps are biochemical pathway 
5 templates which show enzymes connecting metabolites for 
various parts of metabolism. These general pathway 
templates are customized for a given organism by 
highlighting enzymes on a given template which have 
been identified in the genome of the organism. Enzymes 
10 and metabolites are active and yield useful information 
about stoichiometry, structure, alternative names and 
the like, when accessed. 

SWISS-PROT contains detailed information 
about protein function. Accessible information 
15 includes alternate gene and gene product names, 

function, structure and sequence information, relevant 
literature references , and the like. 

LocusLink contains general information about 
the locus where the gene is located and, of relevance, 
20 tissue specificity, cellular location, and implication 
of the gene product in various disease states. 

The Enzyme Nomenclature database can be used 
to compare the gene products of two organisms. Often 
the gene names for genes with similar functions in two 

25 or more organisms are unrelated. When this is the 

case, the E.C. (Enzyme Commission) numbers can be used 
as unambiguous indicators of gene product function. 
The information in the Enzyme Nomenclature database is 
also published in Enzyme Nomenclature (Academic Press, 

30 San Diego, California, 1992) with 5 supplements to 

date, all found in the European Journal of Biochemistry 
(Blackwell Science, Maiden, MA) . 
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Sources of biochemical information include, 
for example, general resources relating to metabolism, 
resources relating specifically to human metabolism, 
and resources relating to the biochemistry, physiology 
5 and pathology of specific human cell types. 

Sources of general information relating to 
metabolism, which were used to generate the human 
reaction databases and models described herein, were 
J.G. Salway, Metabolism at a Glance , 2 nd ed., Blackwell 

10 Science, Maiden, MA (1999) and T.M. Devlin, ed., 

Textbook of Biochemistry with Clinical Correlations , 
4 th ed., John Wiley and Sons, New York, NY (1997). 
Human metabolism-specific resources included J.R. 
Bronk, Human Metabolism: Functional Diversity and 

15 Integration , Addison Wesley Longman, Essex, England 
(1999) . 

The literature used in conjunction with the 
skeletal muscle metabolic models and simulations 
described herein included R. Maughan et al., 

20 Biochemistry of Exercise and Training , Oxford 

University Press, Oxford, England (1997), as well as 
references on muscle pathology such as S. Carpenter et 
al., Pathology of Skeletal Muscle , 2 nd ed. , Oxford 
University Press, Oxford, England (2001), and more 

25 specific articles on muscle metabolism as may be found 
in the Journal of Physiology (Cambridge University 
Press, Cambridge, England) . 

In the course of developing an in silico 
30 model of Homo sapiens metabolism, the types of data 
that can be considered include, for example, 
biochemical information which is information related to 
the experimental characterization of a chemical 
reaction, often directly indicating a protein (s) 
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associated with a reaction and the stoichiometry of the 
reaction or indirectly demonstrating the existence of a 
reaction occurring within a cellular extract; genetic 
information, which is information related to the 
5 experimental identification and genetic 

characterization of a gene(s) shown to code for a 
particular protein (s) implicated in carrying out a 
biochemical event; genomic information, which is 
information related to the identification of an open 

10 reading frame and functional assignment, through 

computational sequence analysis, that is then linked to 
a protein performing a biochemical event; physiological 
information, which is information related to overall 
cellular physiology, fitness characteristics, substrate 

15 utilization, and phenotyping results, which provide 
evidence of the assimilation or dissimilation of a 
compound used to infer the presence of specific 
biochemical event (in particular translocations); and 
modeling information, which is information generated 

20 through the course of simulating activity of Homo 
sapiens cells using methods such as those described 
herein which lead to predictions regarding the status 
of a reaction such as whether or not the reaction is 
required to fulfill certain demands placed on a ■■ 

25 metabolic network. Additional information relevant to 
multicellular organisms that can be- considered includes 
cell type-specif ic or condition-specific gene 
expression information, which can be determined 
experimentally, such as by gene array analysis or from 

30 expressed sequence tag (EST) analysis, or obtained from 
the biochemical and physiological literature. 

The majority of the reactions occurring in 
Homo sapiens reaction networks are catalyzed by 
enzymes/proteins, which are created through the 
35 transcription and translation of the genes found within 
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the chromosome in the cell. The remaining reactions 
occur either spontaneously or through non-enzymatic 
processes. Furthermore, a reaction network data 
structure can contain reactions that add or delete 
5 steps to or from a particular reaction pathway. For 
example, reactions can be added to optimize or improve 
performance of a Homo sapiens model in view of 
empirically observed activity. Alternatively, 
reactions can be deleted to remove intermediate steps 

10 in a pathway when the intermediate steps are not 
necessary to model flux through the pathway. For 
example, if a pathway contains 3 nonbranched steps, the 
reactions can be combined or added together to give a 
net reaction, thereby reducing memory required to store 

15 the reaction network data structure and the 

computational resources required for manipulation of 
the data structure. 

The reactions that occur due to the activity 
of gene-encoded enzymes can be obtained from a genome 

20 database which lists genes identified from genome 

sequencing and subsequent genome annotation. Genome 
annotation consists of the locations of open reading 
frames and assignment of function from homology to 
other known genes or empirically determined activity. 

25 Such a genome database can be acquired through public 
or private databases containing annotated Homo sapiens 
nucleic acid or protein sequences. If desired, a model 
developer can perform a network reconstruction and 
establish the model content associations between the 

30 genes, proteins, and reactions as described, for 
example, in Covert et al. Trends in Biochemical 
Sciences 26:179-186 (2001) and Palsson, WO 00/46405. 



As reactions are added to a reaction network 
35 data structure or metabolic reaction database, those 
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having known or putative associations to the 
proteins/enzymes which enable/catalyze the reaction and 
the associated genes that code for these proteins can 
be identified by annotation. Accordingly, the 
5 appropriate associations for all of the reactions to 
their related proteins or genes or both can be 
assigned. These associations can be used to capture 
the non-linear relationship between the genes and 
proteins as well as between proteins and reactions. In 

10 some cases one gene codes for one protein which then 
perform one reaction. However/ often there are 
multiple genes which are required to create an active 
enzyme complex and often there are multiple reactions 
that can be carried out by one protein or multiple 

15 proteins that can carry out the same reaction. These 
associations capture the logic (i.e. AND or OR 
relationships) within the associations. Annotating a 
metabolic reaction database with these associations can 
allow the methods to be used to determine the effects 

20 of adding or eliminating a particular reaction not only 
at the reaction level, but at the genetic or protein 
level in the context of running a simulation or 
predicting Homo sapiens activity. 

A reaction network data structure of the 
25 invention can be used to determine the activity of one 
or more reactions in a plurality of Homo sapiens 
reactions independent of any knowledge or annotation of 
the identity of the protein that performs the reaction 
or the gene encoding the protein. A model that is 
30 annotated with gene or protein identities can include 
reactions for which a protein or encoding gene is not 
assigned. While a large portion of the reactions in a 
cellular metabolic network are associated with genes in 
the organism' s genome, there are also a substantial 
35 number of reactions included in a model for which there 
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are no known genetic associations. Such reactions can 
be added to a reaction database based upon other 
information that is not necessarily related to genetics 
such as biochemical or cell based measurements or 
5 theoretical considerations based on observed 

biochemical or cellular activity. For example, there 
are many reactions that can either occur spontaneously 
or are not protein-enabled reactions. Furthermore, the 
occurrence of a particular reaction in a cell for which 
10 no associated proteins or genetics have been currently 
identified can be indicated during the course of model 
building by the iterative model building methods of the 
invention. 

The reactions in a reaction network data 

15 structure or reaction database can be assigned to 

subsystems by annotation, if desired. The reactions 
can be subdivided according to biological criteria, 
such as according to traditionally identified metabolic 
pathways (glycolysis, amino acid metabolism and the 

20 like) or according to mathematical or computational 
criteria that facilitate manipulation of a model that 
incorporates or manipulates the reactions. Methods and 
criteria for subdviding a reaction database are 
described in further detail in Schilling et al., CL. 

25 Theor. Biol. 203:249-283 (2000), and in Schuster et 
al., Bioinformatics 18:351-361 (2002). The use of 
subsystems can be advantageous for a number of analysis 
methods, such as extreme pathway analysis, and can make 
the management of model content easier. Although 

30 assigning reactions to subsystems can be achieved 
without affecting the use of the entire model for 
simulation, assigning reactions to subsystems can allow 
a user to search for reactions in a particular 
subsystem which may be useful in performing various 

35 types of analyses. Therefore, a reaction network data 
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structure can include any number of desired subsystems 
including, for example, 2 or more subsystems, 5 or more 
subsystems, 10 or more subsystems, 25 or more 
subsystems or 50 or more subsystems. 

5 The reactions in a 'reaction network data 

structure or metabolic reaction database can be 
annotated with a value indicating the confidence with 
which the reaction is believed to occur in the Homo 
sapiens cell. The level of confidence can be, for 
10 example, a function of the amount and form of 

supporting data that is available. This data can come 
in various forms including published literature, 
documented experimental results, or results of 
computational analyses. Furthermore, the data can 

15 provide direct or indirect evidence for the existence 
of a chemical reaction in a cell based on genetic, 
biochemical, and/or physiological data. 

The invention further provides a computer 
readable medium, containing (a) a data structure 

20 relating a plurality of Homo sapiens reactants to a 
plurality of Homo sapiens reactions, wherein each of 
the Homo sapiens reactions includes a reactant 
identified as a substrate of the reaction, a reactant 
identified as a product of the reaction and a 

25 stoichiometric coefficient relating the substrate and 
the product, and (b) a constraint set for the plurality 
of Homo sapiens reactions. 

Constraints can be placed on the value of any 
of the fluxes in the metabolic network using a 
30 constraint set. These constraints can be 

representative of a minimum or maximum allowable flux 
through a given reaction, possibly resulting from a 
limited amount of an enzyme present. Additionally, the 
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constraints can determine the direction or 
reversibility of any of the reactions or transport 
fluxes in the reaction network data structure. Based 
on the in vivo environment where Homo sapiens lives the 
5 metabolic resources available to the cell for 
biosynthesis of essential molecules for can be 
determined. Allowing the corresponding transport 
fluxes to be active provides the in silico Homo sapiens 
with inputs and outputs for substrates and by-products 
10 produced by the metabolic network. 

Returning to the hypothetical reaction 
network shown in Figure 1, constraints can be placed on 
each reaction in the exemplary format shown in Figure 
2, as follows. The constraints are provided in a 
15 format that can be used to constrain the reactions of 
the stoichiometric matrix shown in Figure 3. The 
format for the constraints used for a matrix or in 
linear programming can be conveniently represented as a 
linear inequality such as 

20 bj <, Vj £ aj :j = l....n (Eq. 1) 

where Vj is the metabolic flux vector, bj is the 
minimum flux value and a^ is the maximum flux value. 
Thus, aj can take on a finite value representing a 
maximum allowable flux through a given reaction or bj 

25 can take on a finite value representing minimum 

allowable flux through a given reaction. Additionally, 
if one chooses to leave certain reversible reactions* or 
transport fluxes to operate in a forward and reverse 
manner the flux may remain unconstrained by setting bj 

30 to negative infinity and aj to positive infinity as 
shown for reaction R 2 in Figure 2. If reactions 
proceed only in the forward reaction bj is set to zero 
while aj is set to positive infinity as shown for 
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reactions R lf R 3 , R 4 , R 5 , and R 6 in Figure 2. As an 
example, to simulate the event of a genetic deletion or 
non-expression of a particular protein, the flux 
through all of the corresponding metabolic reactions 
5 related to the gene or protein in question are reduced 
to zero by setting and bj to be zero. Furthermore, 
if one wishes to simulate the absence of a particular 
growth substrate one can simply constrain the 
corresponding transport fluxes that allow the 
10 metabolite to enter the cell to be zero by setting a^ 
and bj to be zero. On the other hand if a substrate is 
only allowed to enter or exit the cell via transport 
mechanisms, the corresponding fluxes can be properly 
constrained to reflect this scenario. 

15 The ability of a reaction to be actively 

occurring is dependent on a large number of additional 
factors beyond just the availability of substrates. 
These factors , which can be represented as variable 
constraints in the models and methods of the invention 

20 include, for example, the presence of cof actors 

necessary to stabilize the protein/enzyme, the presence 
or absence of enzymatic inhibition and activation 
factors , the active formation of the protein/enzyme 
through translation of the corresponding mRNA 

25 transcript, the transcription of the associated gene(s) 
or the presence of chemical signals and/or proteins 
that assist in controlling these processes that 
ultimately determine whether a chemical reaction is 
capable of being carried out within an organism. Of 

30 particular importance in the regulation of human cell 
types is the implementation of paracrine and endocrine 
signaling pathways to control cellular activities. In 
these cases a cell secretes signaling molecules that 
may be carried far afield to act on distant targets 

35 (endocrine signaling) , or act as local mediators 
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(paracrine signaling) . Examples of endocrine signaling 
molecules include hormones such as insulin, while 
examples of paracrine signaling molecules include 
neurotransmitters such as acetylcholine. These 
5 molecules induce cellular responses through signaling 
cascades that affect the activity of biochemical 
reactions in the cell. Regulation can be represented 
in an in silico Homo sapiens model by providing a 
variable constraint as set forth below. 

10 

Thus, the invention provides a computer 
readable medium or media, including (a) a data 
structure relating a plurality of Homo sapiens 
reactants to a plurality of Homo sapiens reactions, 

15 wherein each of the reactions includes a reactant 

identified as a substrate of the reaction, a reactant 
identified as a product of the reaction and a 
stoichiometric coefficient relating the substrate and 
the product, and wherein at least one of the reactions 

20 is a regulated reaction; and (b) a constraint set for 
the plurality of reactions, wherein the constraint set 
includes a variable constraint for the regulated 
reaction. 

As used herein, the term "regulated, " when 
25 used in reference to a reaction in a data structure, is 
intended to mean a reaction that experiences an altered 
flux due to a change in the value of a constraint or a 
reaction that has a variable constraint. 

As used herein, the term "regulatory 
30 reaction" is intended to mean a chemical conversion or 
interaction that alters the activity of a protein, 
macromolecule or enzyme. A chemical conversion or 
interaction can directly alter the activity of a 
protein, macromolecule or enzyme such as occurs when 
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the protein, macromolecule or enzyme is 
post-translationally modified or can indirectly alter 
the activity of a protein, macromolecule or enzyme such 
as occurs when a chemical conversion or binding event 
5 leads to altered expression of the protein, 

macromolecule or enzyme. Thus, transcriptional or 
translational regulatory pathways can indirectly alter 
a protein, macromolecule or enzyme or an associated 
reaction. Similarly, indirect regulatory reactions can 

10 include reactions that occur due to downstream 

components or participants in a regulatory reaction 
network. When used in reference to a data structure or 
in silico Homo sapiens model, the term is intended to 
mean a first reaction that is related to a second 

15 reaction by a function that alters the flux through the 
second reaction by changing the value of a constraint 
on the second reaction. 

As used herein, the term "regulatory data 
structure" is intended to mean a representation of an 

20 event, reaction or network of reactions that activate 
or inhibit a reaction, the representation being in a 
format that can be manipulated or analyzed. An event 
that activates a reaction can be an event that 
initiates the reaction or an event that increases the 

25 rate or level of activity for the reaction. An event" 
that inhibits a reaction can be an event that stops the 
reaction or an event that decreases the rate or level 
of activity for the reaction. Reactions that can be 
represented in a regulatory data structure include, for 

30 example, reactions that control expression of a 

macromolecule that in turn, performs a reaction such as 
transcription and translation reactions, reactions that 
lead to post translational modification of a protein or 
enzyme such as phophorylation, dephosphorylation, 

35 prenylation, methylation, oxidation or covalent 



WO 03/082214 PCT/US03/09751 

42 

modification, reactions that process a protein or 
enzyme such as removal of a pre- or pro-sequence, 
reactions that degrade a protein or enzyme or reactions 
that lead to assembly of a protein or enzyme. 

5 As used herein, the term "regulatory event" 

is intended to mean a modifier of the flux through a 
reaction that is independent of the amount of reactants 
available to the reaction. A modification included in 
the term can be a change in the presence, absence, or 

10 amount of an enzyme that performs a reaction. A 
modifier included in the term can be a regulatory 
reaction such as a signal transduction reaction or an 
environmental condition such as a change in pH,* 
temperature, redox potential or time. It will be 

15 understood that when used in reference to an in silica 
Homo sapiens model or data structure a regulatory event 
is intended to be a representation of a modifier of the 
flux through a Hojno sapiens reaction that is 
independent of the amount of reactants available to the 

20 reaction. 

The effects of regulation on one or more 
reactions that occur in Homo sapiens can be predicted 
using an in silico Homo sapiens model of the invention. 

25 Regulation can be taken into consideration in the 
context of a particular condition being examined by 
providing a variable constraint for the reaction in an 
in silico Homo sapiens model. Such constraints 
constitute condition-dependent constraints. A data 

30 structure can represent regulatory reactions as Boolean 
logic statements (Reg-reaction) . The variable takes on 
a value of 1 when the reaction is available for use in 
the reaction network and will take on a value of 0 if 
the reaction is restrained due to some regulatory 

35 feature. A series of Boolean statements can then be 
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introduced to mathematically represent the regulatory 
network as described for example in Covert et al. J. 
Theor; Biol. 213:73-88 (2001). For example, in the 
case of a transport reaction (A_in) that imports 
metabolite A, where metabolite A inhibits reaction R2 
as shown in Figure 4, a Boolean rule can state that: 

Reg-R2 = IF NOT (A_in) . (Eq. 2) 

This statement indicates that reaction R2 can occur if 
reaction A_in is not occurring (i.e. if metabolite A is 
not present) . Similarly, it is possible to assign the 
regulation to a variable A which would indicate an 
amount of A above or below a threshold that leads to 
the inhibition of reaction R2. Any function that 
provides values for variables corresponding to each of 
the reactions in the biochemical reaction network can 
be used to represent a regulatory reaction or set of 
regulatory reactions in a regulatory data structure. 
Such functions can include, for example, fuzzy logic, 
heuristic rule-based descriptions, differential 
equations or kinetic equations detailing system 
dynamics . 

A reaction constraint placed on a reaction 
can be incorporated into -an in silico Homo sapiens 
model using the following general equation: 

(Reg-Reaction) *b.j <, v 3 z a^* (Reg-Reaction) 

: (Eq. 3) 

j = 1 .... n 

For the example of reaction R2 this equation is written 
as follows: 



(0)*Reg-R2 £ R2 s <eo)*Reg-R2. (Eq. 4) 
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Thus, during the course of a simulation, depending upon 
the presence or absence of metabolite A in the interior 
of the cell where reaction R2 occurs, the value for the 
upper boundary of flux for reaction R2 will change from 
0 to infinity, respectively. 

With the effects of a regulatory event or 
network taken into consideration by a constraint 
function and the condition-dependent constraints set to 
an initial relevant value, the behavior of the Homo 
sapiens reaction network can be simulated for the 
conditions considered as set forth below. 

Although regulation has been exemplified 
above for the case where a variable constraint is 
dependent upon the outcome of a reaction in the data 
structure, a plurality of variable constraints can be 
included in an in silico Homo sapiens model to 
represent regulation of a plurality of reactions. 
Furthermore, in the exemplary case set forth above, the 
regulatory structure includes a general control stating 
that a reaction is inhibited by a particular 
environmental condition. Using a general control of 
this type, it is possible to incorporate molecular 
mechanisms and additional detail into the regulatory 
structure that is responsible for determining the 
active nature of a particular chemical reaction within 
an organism. 

Regulation can also be simulated by a model 
of the invention and used to predict a Homo sapiens 
physiological function without knowledge of the precise 
molecular mechanisms involved in the reaction network 
being modeled. Thus, the model can be used to predict, 
in silico r overall regulatory events or causal 
relationships that are not apparent from in vivo 
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observation of any one reaction in a network or whose 

in vivo effects on a particular reaction are not known. 
Such overall regulatory effects can include those that 

result from overall environmental conditions such as 

changes in pH r temperature, redox potential, or the 
passage of time. 

The in silico Homo sapiens model and methods 
described herein can be implemented on any conventional 
host computer system, such as those based on Intel. RTM. 
microprocessors and running Microsoft Windows operating 
systems. Other systems, such as those using the UNIX or 
LINUX operating system and based on IBM. RTM. , DEC. RTM. 
or Motorola. RTM. microprocessors are also contemplated. 
The systems and methods described herein can also be 
implemented to run on client-server systems and 
wide-area networks, such as the Internet. 

Software to implement a method or model of 
the invention can be written in any well-known computer 
language, such as Java, C, C++, Visual Basic, FORTRAN 
or COBOL and compiled using any well-known compatible 
compiler. The software of the invention normally runs 
from instructions stored in a memory on a host computer 
system. A memory or computer readable medium can be a . 
hard disk, floppy disc, compact disc, magneto-optical 
disc, Random Access Memory, Read Only Memory or Flash 
Memory. The memory or computer readable medium used in 
the invention can be contained within a single computer 
or distributed in a network. A network can be any of a 
number of conventional network systems known in the art 
such as a local area network (LAN) or a wide area 
network (WAN) . Client-server environments, database 
servers and networks that can be used in the invention 
are well known in the art. For example, the database 
server can run on an operating system such as UNIX, 
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running a relational database management system, a 
World Wide Web application and a World Wide Web server. 
Other types of memories and computer readable media are 
also contemplated to function within the scope of the 
5 invention. 

A database or data structure of the invention 
can be represented in a markup language format 
including, for example, Standard Generalized Markup 
Language (SGML), Hypertext markup language (HTML) or 

10 Extensible Markup language (XML) . Markup languages can 
be used to tag the information stored in a database or 
data structure of the invention, thereby providing 
convenient annotation and transfer of data between 
databases and data structures. In particular, an XML 

15 format can be useful for structuring the data 

representation of reactions, reactants and their 
annotations; for exchanging database contents, for 
example, over a network or internet; for updating 
individual elements using the document object model; or 

20 ' for providing differential access to multiple users for 
different information content of a data base or data 
structure of the invention. XML programming methods 
and editors for writing XML code are known in the art 
as described, for example, in Ray, "Learning XML" 

25 O'Reilly and Associates, Sebastopol, CA (2001). 

A set of constraints can be applied to a 
reaction network data structure to simulate the flux of 
mass through the reaction network under a particular 
set of environmental conditions specified by a 
30 constraints set. Because the time constants 

characterizing metabolic transients and/or metabolic 
reactions are typically very rapid, on the order of 
milli-seconds to seconds, compared to the time 
constants of cell growth on the order of hours to days, 
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the transient mass balances can be simplified to only 
consider the steady state behavior. Referring now to an 
example where the reaction network data structure is a 
stoichiometric matrix, the steady state mass balances 
5 can be applied using the following system of linear 
equations 

S • v = 0 (Eq. 5) 

where S is the stoichiometric matrix as defined above 
and v is the flux vector. This equation defines the 

10 mass, energy, and redox potential constraints placed on 
the metabolic network as a result of stoichiometry . 
Together Equations 1 and 5 representing the reaction 
constraints and mass balances, respectively, 
effectively define the capabilities and constraints of 

15 the metabolic genotype and the organism's metabolic 

potential. All vectors, v, that satisfy Equation 5 are 
said to occur in the mathematical nullspace of S. 
Thus, the null space defines steady-state metabolic 
flux distributions that do not violate the mass, 

20 energy, or redox balance constraints. Typically, the 
number of fluxes is greater than the number of mass 
balance constraints, thus a plurality of flux 
distributions satisfy the mass balance constraints and 
occupy the null space. -The" null" space, which defines 

25 the feasible set of metabolic flux distributions, is 
further reduced in size by applying the reaction 
constraints set forth in Equation 1 leading to a 
defined solution space. A point in this space 
represents a flux distribution and hence a metabolic 

30 phenotype for the network. An optimal solution within 
the set of all solutions can be determined using 
mathematical optimization methods when provided with a 
stated objective and a constraint set. The calculation 
of any solution constitutes a simulation of the model.- 
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Objectives for activity of a human cell can 
be chosen. While the overall objective of a 
multi-cellular organism may be growth or reproduction, 
individual human cell types generally have much more 
5 complex objectives, even to the seemingly extreme 

objective of apoptosis (programmed cell death), which 
may benefit the organism but certainly not the 
individual cell. For example, certain cell types may 
have the objective of maximizing energy production, 

10 while others have the objective of maximizing the 
production of a particular hormone, extracellular 
matrix component, or a mechanical property such as 
contractile force. In cases where cell reproduction is 
slow, such as human skeletal muscle, growth and its 

15 effects need not be taken into account. In other 
cases, biomass composition and growth rate could be 
incorporated into a "maintenance" type of flux, where 
rather than optimizing for growth, production of 
precursors is set at a level consistent with 

20 experimental knowledge and a different objective is 
optimized. 



Certain cell types, including cancer cells, 
can be viewed as having an objective of maximizing cell 
growth. Growth can be defined in terms of biosynthetic 

25 requirements based on literature values of biomass 

composition or experimentally determined values such as 
those obtained as described above. Thus, biomass 
generation can be defined as an exchange reaction that 
removes intermediate metabolites in the appropriate 

30 ratios and represented as an objective function. In 
addition to draining intermediate metabolites this 
reaction flux can be formed to utilize energy molecules 
such as ATP, NADH and NADPH so as to incorporate any 
maintenance requirement that must be met. This new 

35 reaction flux then becomes another constraint /balance 
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equation that the system must satisfy as the objective 
function. Using the stoichiometric matrix of Figure 3 
as an example, adding such a constraint is analogous to 
adding the additional column V growth to the 
5 stoichiometric matrix to represent fluxes to describe 
the production demands placed on the metabolic system. 
Setting this new flux as the objective function and 
asking the system to maximize the value of this flux 
for a given set of constraints on all the other fluxes 
10 is then a method to simulate the growth of the 
organism. 

Continuing with the example of the 
stoichiometric matrix applying a constraint set to a 
reaction network data structure can be illustrated as 
15 follows. The solution to equation 5 can be formulated 
as an optimization problem, in which the flux 
distribution that minimizes a particular objective is 
found. Mathematically, this optimization problem can 
be stated as: 

20 Minimize Z (Eq. 6) 

where z = £ c /* v i (Eq- D 



where Z is the objective which is represented as a 
25 linear combination of metabolic fluxes v* using the 
weights c A in this linear combination. The 
optimization problem can also be stated as the 
equivalent maximization problem; i.e. by changing the 
sign on Z. Any commands for solving the optimazation 
30 problem can be used including, for example, linear 
programming commands. 
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A computer system of the invention can 
further include a user interface capable of receiving a 
representation of one or more reactions. A user 
interface of the invention can also be capable of 
sending at least one command for modifying the data 
structure, the constraint set or the commands for 
applying the constraint set to the data representation, 
or a combination thereof. The interface can be a 
graphic user interface having graphical means for 
making selections such as menus or dialog boxes. The 
interface can be arranged with layered screens 
accessible by making selections from a main screen. 
The user interface can provide access to other 
databases useful in the invention such as a metabolic 
reaction database or links to other databases having 
information relevant to the reactions or reactants in 
the reaction network data structure or to Homo sapiens 
physiology. Also, the user interface can display a 
graphical representation of a reaction network or the 
results of a simulation using a model of the invention. 

Once an initial reaction network data 
structure and set of constraints has been created, this 
model can be tested by preliminary simulation. During 
preliminary simulation, gaps in the network or 
"dead-ends" in which a metabolite can be produced but 
not consumed or where a metabolite can be consumed but 
not produced can be identified. Based on the results 
of preliminary simulations areas of the metabolic 
reconstruction that require an additional reaction can 
be identified. The determination of these gaps can be 
readily calculated through appropriate queries of the 
reaction network data structure and need not require 
the use of simulation strategies, however, simulation 
would be an alternative approach to locating such gaps. 
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In the preliminary simulation testing and 
model content refinement stage the- existing model is 
subjected to a series of functional tests to determine 
if it can perform basic requirements such as the 
5 ability to produce the required biomass constituents 
and- generate predictions concerning the basic 
physiological characteristics of the particular cell 
type being modeled. The more preliminary testing that 
is conducted the higher the quality of the model that 
10 will be generated. Typically, the majority of the 

simulations used in this stage of development will be 
single optimizations. A single optimization can be 
used to calculate a single flux distribution 
demonstrating how metabolic resources are routed 
15 determined from the solution to one optimization 

problem. An optimization problem can be solved using 
linear programming as demonstrated in the Examples 
below. The result can be viewed as a display of a flux 
distribution on a reaction map. Temporary reactions 
20 can be added to the network to determine if they should 
be included into the model based on modeling/simulation 
requirements. 

Once a model of the invention is sufficiently 
complete with respect to the content of the reaction 

25 network data structure according to the criteria set 

forth above, the model can be used to simulate activity 
of one or more reactions in a reaction network. The 
results of a simulation can be displayed in a variety 
of formats including, for example, a table, graph, 

30 reaction network, flux distribution map or a phenotypic 
phase plane graph. 



Thus, the invention provides a method for 
predicting a Homo sapiens physiological function. The 
method includes the steps of (a) providing a data 
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structure relating a plurality of Homo sapiens 
reactants to a plurality of Homo sapiens reactions, 
wherein each of the Homo sapiens reactions includes a 
reactant identified as a substrate of the reaction, a 
5 reactant identified as a product of the reaction and a 
stoichiometric coefficient relating said substrate and 
said product; (b) providing a constraint set for the 
plurality of Homo sapiens reactions; (c) providing an 
objective function, and (d) determining at least one 
10 flux distribution that minimizes or maximizes the 

objective function when the constraint set is applied 
to the data structure, thereby predicting a Homo 
sapiens physiological function. 

A method for predicting a Homo sapiens 
15 physiological function can include the steps of (a) 

providing a data structure relating a plurality of Homo 
■ sapiens reactants to a plurality of Homo sapiens 
reactions, wherein each of the Homo sapiens reactions 
includes a reactant identified as a substrate of the 
20 reaction, a reactant identified as a product of the 

reaction and a stoichiometric coefficient relating the 
substrate and the product, and wherein at least' one of 
the reactions is a regulated reaction; (b) providing a 
constraint set for the plurality of reactions, wherein 
25 the constraint set includes a variable constraint for 
the regulated reaction; (c) providing a 
condition-dependent value to the variable constraint; 
(d) providing an objective function, and (e) 
determining at least one flux distribution that 
30 minimizes or maximizes the objective function when the 
constraint set is applied to the data structure, 
thereby predicting a Homo sapiens physiological 
function. 
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As used herein, the term "physiological 
function/' when used in reference to Homo sapiens, is 
intended to mean an activity of a Homo sapiens cell as 
a whole. An activity included in the term can be the 
5 magnitude or rate of a change from an initial state of 
a Homo sapiens cell to a final state of the Homo 
sapiens cell. An activity included in the term can be, 
for example, growth, energy production, redox 
equivalent production, biomass production, development, 
10 or consumption of carbon nitrogen, sulfur, phosphate, 
hydrogen or oxygen. An activity can also be an output 
of a particular reaction that is determined or 
predicted in the context of substantially all of the 
reactions that affect the particular reaction in a Homo 
15 sapiens cell or substantially all of the reactions that 
occur in a Homo sapiens cell (e.g. muscle contraction). 
Examples of a particular reaction included in the term 
are production of biomass precursors, production of a 
protein, production of an amino acid, production of a 
20 purine, production of a pyrimidine, production of a 
lipid, production of a fatty acid, production of a 
cof actor or transport of a metabolite. A physiological 
function can include an emergent property which emerges 
from the whole but not from the sum of parts where the 
25 parts are observed in isolation (see for example, 
Palsson, Nat. Biotech 18 : 1147-1150 - (2000) ) . 

A physiological function of Homo sapiens 
reactions can be determined using phase plane analysis 
of flux distributions. Phase planes are 
30 representations of the feasible set which can be 

presented in two or three dimensions. As an example, 
two parameters that describe the growth conditions such 
as substrate and oxygen uptake rates can be defined as 
two axes of a two-dimensional space. The optimal flux 
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distribution can be calculated from a reaction network 
data structure and a set of constraints as set forth 
above for all points in this plane by repeatedly 
solving the linear programming problem while adjusting 
5 the exchange fluxes defining the two-dimensional space. 
A finite number of qualitatively different metabolic 
pathway utilization patterns can be identified in such 
a plane, and lines can be drawn to demarcate these 
regions. The demarcations defining the regions can be 
10 determined using shadow prices of linear optimization 
as described, for example inChvatal, Linear 
Programming New York, W.H. Freeman and Co. (1983). The 
regions are referred to as regions of constant shadow 
price structure. The shadow prices define the 
15 intrinsic value of each reactant toward the objective 
function as a number that is either negative, zero, or 
positive and are graphed according to the uptake rates 
represented by the x and y axes. When the shadow 
prices become zero as the value of the uptake rates are 
changed there is a qualitative shift in the optimal 
reaction network. 



20 



One demarcation line in the phenotype phase 
plane is defined as the line of optimality (LO) . This, 
line represents the optimal relation between respective 
25 metabolic fluxes. The LO can be- identified by varying 
the x-axis flux and calculating the optimal y-axis flux 
with the objective function defined as the growth flux 
. From the phenotype phase plane analysis the 
conditions under which a desired activity is optimal 
can be determined. The maximal uptake rates lead to 
the definition of a finite area of the plot that is the 
predicted outcome of a reaction network within the 
environmental conditions represented by the constraint 
set. Similar analyses can be performed in multiple 
35 dimensions where each dimension on the plot corresponds 



30 
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to a different uptake rate. These and other methods 
for using phase plane analysis, such as those described 
in Edwards et al., giofcech, Bioenq. 77:27-36(2002), can 
be used to analyze the results of a simulation using an 
5 in silico Homo sapiens model of the invention. 

A physiological function of Homo sapiens can 
also be determined using a reaction map to display a 
flux distribution. A reaction map of Homo sapiens can 
be used to view reaction networks at a variety of 
10 levels. In the case of a cellular metabolic reaction 
network a reaction map can contain the entire reaction 
complement representing a global perspective. 
Alternatively, a reaction map can focus on a particular 
region of metabolism such as a region corresponding to 
15 a reaction subsystem described above or even on an 
individual pathway or reaction. 

Thus, the invention provides an apparatus 
that produces a representation of a Homo sapiens 
physiological function, wherein the representation is 

20 produced by a process including the steps of: (a) 

providing a data structure relating a plurality of Homo 
sapiens reactants to a plurality of Homo sapiens 
reactions, wherein each of the Homo sapiens reactions 
includes a reactant identified as a substrate of the 

25 reaction, a reactant identified as a product of the 

reaction and a stoichiometric coefficient relating said 
substrate and said product; (b) providing a constraint 
set ..for the plurality of Homo sapiens reactions; (c) 
providing an objective function; (d) determining at 

30 least one flux distribution that minimizes or maximizes 
the objective function when the constraint set is 
applied to the data structure, thereby predicting a 
Homo sapiens physiological function, and (e) producing 
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a representation of the activity of the one or more 
Homo sapiens reactions. 

The methods of the invention can be used to 
determine the activity of a plurality of Homo sapiens 
5 reactions including, for example, biosynthesis of an 
amino acid, degradation of an amino acid, biosynthesis 
of a purine, biosynthesis of a pyrimidine, biosynthesis 
of a lipid, metabolism of a fatty acid, biosynthesis of 
a cof actor, transport of a metabolite and metabolism of 
10 an alternative carbon source. In addition, the methods 
can be used to determine the activity of one or more of 
the reactions described above or listed in Table 1. 

The methods of the invention can be used to 
determine a phenotype of a Homo sapiens mutant. The 
15 activity of one or more Homo sapiens reactions can be 
determined using the methods described above, wherein 
the reaction network data structure lacks one or more 
gene-associated reactions that occur in Homo sapiens. 
Alternatively, the methods can be used to determine the 
20 activity of one or more Homo sapiens reactions when a 
reaction that does not naturally occur in Homo sapiens 
is added to the reaction network data structure. 
Deletion of a gene can also be represented in a model 
of the invention by constraining the flux through the 
25 reaction to zero, thereby allowing the reaction to 
remain within the data structure. Thus, simulations 
can be made to predict the effects of adding or 
removing genes to or from Homo sapiens. The methods 
can be particularly useful for determining the effects 
30 of adding or deleting a gene that encodes for a gene 
product that performs a reaction in a peripheral 
metabolic pathway. 
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A drug target or target for any other agent 
that affects Homo sapiens function can be predicted 
using the methods of the invention. Such predictions 
can be made by removing a reaction to simulate total 
5 inhibition or prevention by a drug or agent. 

Alternatively, partial inhibition or reduction in the 
activity a particular reaction can be predicted by 
performing the methods with altered constraints. For 
example, reduced activity can be introduced into a 
10 model of the invention by altering the a 3 or b 3 values 
for the metabolic flux vector of a target reaction to 
reflect a finite maximum or minimum flux value 
corresponding to the level of inhibition.' Similarly, 
the effects of activating a reaction, by initiating or 
15 increasing the activity of the reaction, can be 

predicted by performing the methods with a reaction 
network data structure lacking a particular reaction or 
by altering the a d or b 3 values for the metabolic flux 
vector of a target reaction to reflect a maximum or 
20 minimum flux value corresponding to the level of 

activation. The methods can be particularly useful for 
identifying a target in a peripheral metabolic pathway. 

Once a reaction has been identified for which 
activation or inhibition produces a desired effect on 

25 Homo sapiens function, an enzyme or macromolecule that 
performs the reaction in Homo sapiens or a gene that 
expresses the enzyme or macromolecule can be identified 
as a target for a drug or other agent. A candidate 
compound for a target identified by the methods of the 

30 invention can be isolated or synthesized using known 
methods. Such methods for isolating or synthesizing 
compounds can include, for example, rational design 
based on known properties of the target (see, for 
example, DeCamp et al., Protein Engineering Principles 

35 and Practice , Ed. Cleland and Craik, Wiley-Liss, New 
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York, pp. 467-506 (1996)), screening the target against 
combinatorial libraries of compounds (see for example, 
Houghten et al., Nature , 354, 84-86 (1991); Dooley et 
al., Science, 266, 2019-2022 (1994), which describe an 

5 iterative approach, or R. Houghten et al. 

PCT/US91/08694 and U.S. Patent 5,556,762 which describe 
the positional-scanning approach) , or a combination of 
both to obtain focused libraries. Those skilled in the 
art will know or will be able to routinely determine 

10 assay conditions to be used in a screen based on 

properties of the target or activity assays known in 
the art. 

A candidate drug or agent, whether identified 
by the methods described above or by other methods 
15 known in the art, can be validated using an in silico 
Homo sapiens model or method of the invention. The 
effect of a candidate drug or agent on Homo sapiens 
physiological function can be predicted based on the 
activity for a target in the presence of the candidate 
20 drug or agent measured in vitro or in vivo. This 
activity can be represented in an in silico Homo 
sapiens model by adding a reaction, to the model, 
removing a reaction from the model or adjusting a 
constraint for a reaction in the model to reflect the 
25 measured effect of the candidate drug or agent on the 
activity of the reaction. By running a simulation , 
under these conditions the holistic effect of the 
candidate drug or agent on Homo sapiens physiological 
function can be predicted. 

30 The methods of the invention can be used to 

determine the effects of one or more environmental 
components or conditions on an activity of a Homo 
sapiens cell. As set forth above an exchange reaction 
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can be added to a reaction network data structure 
corresponding to uptake of an environmental component, 
release of a component to the environment, or other 
environmental demand. The effect of the environmental 
5 component or condition can be further investigated by 
running simulations with adjusted a., or bj values for 
the metabolic flux vector of the exchange reaction 
target reaction to reflect a finite maximum or minimum 
flux value corresponding to the effect of the 

10 environmental component or condition. The 

environmental component can be, for example an 
alternative carbon source or a metabolite that when 
added to the environment of a Homo sapiens cell can be 
taken up and metabolized. The environmental component 

15 can also be a combination of components present for 
example in a minimal medium composition. Thus, the 
methods can be used to determine an optimal or minimal 
medium composition that is capable of supporting a 
particular activity of Homo sapiens. 

20 The invention further provides a method for 

determining a set of environmental components to 
achieve a desired activity for Homo sapiens. The 
method includes the steps of (a) providing a data 
structure relating a plurality of Homo sapiens 

25 reactants to a plurality of Homo sapiens reactions, 

wherein each of the Homo sapiens reactions includes a 
reactant identified as a substrate of the reaction, a 
reactant identified as a product of the reaction and a 
stoichiometric coefficient relating the substrate and 

30 the product; (b) providing a constraint set for the 
.plurality of Homo sapiens reactions; (c) applying the 
constraint set to the data representation, thereby 
determining the activity of one or more Homo sapiens 
reactions (d) determining the activity of one or more 
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Homo sapiens reactions according to steps (a) through 
(c) , wherein the constraint set includes an upper or 
lower bound on the amount of an environmental component 
and (e) repeating steps (a) through (c) with a changed 
constraint set, wherein the activity determined in step 
(e) is improved compared to the activity determined in 
step (d) . 

The following examples are intended to 
illustrate but not limit the present invention. 

EXAMPLE I 



This example shows the construction of a 
universal Homo sapiens metabolic reaction database, a 
Homo sapiens core metabolic reaction database and a 
Homo sapiens muscle cell metabolic reaction database. 
15 This example also shows the iterative model building 
process used to generate a Homo sapiens core metabolic 
model and a Homo sapiens muscle cell metabolic model. 

A universal Homo sapiens reaction database 
was prepared from the genome databases and biochemical 
20 literature. The reaction database shown in Table 1 
contains the following information: 

Locus ID - the locus number of the gene found 
in the LocusLink website. 

Gene Ab. - various abbreviations which are 

25 used for the gene. 

Reaction Stoichiometry - includes all 
metabolites and direction of the reaction, as well as 

reversibility. 

E.C. - The Enzyme Commission number. 
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Additional information included in the 
universal reaction database, although not shown in 
Table 1, included the chapter of Salway, supra (1999), 
where relevant reactions were found; the cellular 
location, if the reaction primarily occurs in a given 
compartment; the SWISS PROT identifier, which can be 
used to locate the gene record in SWISS PROT; the full 
name of the gene at the given locus; the chromosomal 
location of the gene; the Mendelian Inheritance in Man 
(MIM) data associated with the gene; and the tissue 
type, if the gene is primarily expressed in a certain 
tissue. Overall, 1130 metabolic enzyme- or 
transporter-encoding genes were included in the 
universal reaction database. 

Fifty-nine reactions in the universal 
reaction database were identified and included based on 
biological data as found in Salway supra (1999), 
currently without genome annotation. Ten additional 
reactions, not described in the biochemical literature 
or genome annotation, were subsequently included in the 
reaction database following preliminary simulation 
testing and model content refinement. These 69 
reactions are shown at the end of Table 1. 

From the universal Homo sapiens reaction 
database shown in Table 1, a core metabolic reaction 
database was established, which included core metabolic 
reactions as well as some amino acid and fatty acid 
metabolic reactions, as described in Chapters 1, 3, 4, 
7 , 9, 10, 13, 17, 18 and 44 of J.G. Salway, Metabolism 
at a Glance , 2 nd ed. , Blackwell Science, Maiden, MA 
(1999). The core metabolic reaction database included 
211 unique reactions, accounting for 737 genes in the 
Homo sapiens genome. The core metabolic reaction 
database was used, although not in its entirety, to 
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create the core metabolic model described in Example 
II. 

To allow for the modeling of muscle cells, 
the core reaction database was expanded to include 446 
5 unique reactions, accounting for 889 genes in the Homo 
sapiens genome. This skeletal muscle metabolic 
reaction database was used to create the skeletal 
muscle metabolic model described in Example II. 

Once the core and muscle cell metabolic 
10 reaction databases were compiled, the reactions were 
represented as a metabolic network data structure, or 
"stoichiometric input file." For example, the core 
metabolic network data structure shown in Table 2 
contains 33 reversible reactions, 31 non-reversible 
15 reactions, 97 matrix columns and 52 unique enzymes. 
Each reaction in Table 2 is represented so as to 
indicate the substrate or substrates (a negative 
number) and the product or products (a positive 
number) ; the stoichiometry; the name of each reaction 
20 (the term following the zero) ; and whether the reaction 
is reversible (an R following the reaction name) . A 
metabolite that appears in the mitochondria is 
indicated by an "m, " and a metabolite that appears in 
the extracellular space is indicated by an "ex." 

25 To perform a preliminary simulation or to 

simulate a physiological condition, a set of inputs and 
outputs has to be defined and the network objective 
• function specified. To calculate the maximum ATP 
production of the Homo sapiens core metabolic network 

30 using glucose as a carbon source, a non-zero uptake 
value for glucose was assigned and ATP production was 
maximized as the objective function, using the 
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representation shown in Table 2. The network 1 s 
performance was examined by optimizing for the. given 
objective function and the set of constraints defined 
in the input file, using flux balance analysis methods. 
5 The model was refined in an iterative manner by 
examining the results of the simulation and 
implementing the appropriate changes. 

Using this iterative procedure, two metabolic 
reaction networks were generated, representing human 
10 core metabolism and human skeletal muscle cell 
metabolism. 

EXAMPLE II 

This example shows how human metabolism can 
be accurately simulated using a Homo sapiens core 
15 metabolic model. 

The human core metabolic reaction database 
shown in Table 3 was used in simulations of human core 
metabolism. This reaction database "contains a total of 
20 65 reactions, covering the classic biochemical pathways 
of glycolysis, the pentose phosphate pathway, the 
tricitric acid cycle, oxidative phosphorylation, 
glycogen storage, the malate /aspartate shuttle, the 
glycerol phosphate shuttle, and plasma and 
25 mitochondrial membrane transporters. The reaction 
network was divided into three compartments: the 
cytosol, mitochondria, and the extracellular space. 
The total number of metabolites in the network is 50, 
of which 35 also appear in the mitochondria. This core 
30 metabolic network accounts for 250 human genes. 
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To perform simulations using the core 
metabolic network, network properties such as the P/O 
ratio were specified using Salway, supra (1999) as a 
reference. Oxidation of NADH through the Electron 
5 Transport System (ETS) was set to generate 2.5 ATP 

molecules (i.e. a P/O ratio of 2.5 for NADH), and that 
of FADH 2 was set to 1.5 ATP molecules (i.e. a P/O ratio 
of 1.5 for FADH 2 ) . 

Using the core metabolic network, aerobic and 
10 anaerobic metabolisms were simulated in silico. 

Secretion of metabolic by-products was in agreement 
with the known physiological parameters. Maximum yield 
of all 12 precursor-metabolites, (glucose-6-phosphate, 
fructose-6-phosphate, ribose-5-phosphate, 
15 erythrose-4-phosphate, triose phosphate, 

3-phosphoglycerate, phosphoenolpyruvate, pyruvate, 
acetyl CoA, a-ketoglutarate, succinyl CoA, and 
oxaloacetate) was examined and none found to exceed the 
values of its theoretical yield. 

20 Maximum ATP yield was also examined in the 

cytosol and mitochondria. Salway, supra (1999) reports 
that in the absence of membrane proton-coupled 
transport systems," the energy yield is 38 ATP molecules 
per molecule of glucose and otherwise 31 ATP molecules 

25 per molecule of glucose. The core metabolic model 
demonstrated the same values as described by Salway 
supra (1999) . Energy yield in the mitochondria was 
determined to be 38 molecules of ATP per glucose 
molecule. This is equivalent to production of energy 
30 in the absence of proton-couple transporters across 
mitochondrial membrane since all the protons were 
utilized only in oxidative phosphorylation. In the 
cytosol, energy yield was calculated to be 30.5 
molecules of ATP per glucose molecule. This value 
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reflects the cost of metabolite exchange across the 
mitochondrial membrane as described by Sa^Lway, supra 
(1999). 



EXAMPLE III 



5 This example shows how human muscle cell 

metabolism can be accurately simulated under various 
physiological and pathological conditions using a Homo 
sapiens muscle cell metabolic model. 

As described in Example I, the core metabolic 
10 model was extended to also include all the major 

reactions occurring in the skeletal muscle cell, adding 
new functions to the classical metabolic pathways found 
in the core model , such as fatty acid synthesis and 
(3-oxidation, triacylglycerol and phospholipid 
15 formation, and amino acid metabolism. Simulations were 
performed using the muscle cell reaction database shown 
in Table 4. The biochemical reactions were again 
compartmentalized into cytosolic and mitochondrial 
compartments. 

20 To simulate physiological behavior of human 

skeletal muscle cells, an objective function had to be 
defined. Growth of muscle cells-occurs in time scales 
of several hours to days. The time scale of interest 
in the simulation, however, was in the order of several 
25 to tens of minutes, reflecting the time period of 

metabolic changes during exercise. Thus, contraction 
(defined as, and related to energy production) was 
chosen to be the objective function, and no ■ additional 
constraints were imposed to represent growth demands in 
30 the cell. 
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To study and test the behavior of the 
network, twelve physiological cases (Table 5) and five 
disease cases (Table 6) were examined. The input and 
output of metabolites were specified as indicated in 
Table 5, and maximum energy production and metabolite 
secretions were calculated and taken into account. 



Table 5 



Metabolite 
Exchange 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


Glucose 


I 


I 






I 


I 














02 


I 




I 




T 
1 




T 

■L 




I 




I 




Palmitate 


I 


I 














I 


I 






Glycogen 


I 


I 


I 


I 


















Phosphocrea 
tine 


I 


I 


















I 


I 


Triacylgly- 
cerol 


I 


I 










I 


I 










Isoleucine 


I 


I 






















Valine 


I 


I 






















yrate . 


























Pyruvate 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


Lactate 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


Albumin 


0 


0 


0 


0 


0 


0 


1 0 


0 


0 


0 


0 


0 
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Disease 


Enzyme Deficiency 


Reaction 
Cohs trained 


McArdle's disease 


phosphorylase 


GBEl 


Tarui' s disease 


phosphofructokianse 


PFKL 


Phosphoglycerate 
kinase deficiency 


phosphoglycerate 
kinase 


PGK1R 


Phosphoglycerate 
mutase deficiency 


phosphoglycerate 
mutase 


PGAM3R 


Lactate 

dehydrogenase 

deficiency 


Lactate dehyrogenase 


LDHAR 



10 



The skeletal muscle model was tested for 
utilization of various carbon sources available during 
various stages of exercise and food starvation (Table 
15 5) . The by-product secretion of the network in an 

aerobic to anaerobic shift was qualitatively compared 
to physiological outcome of exercise and found to 
exhibit the same general features such as secretion of 
- fermentative by-products and lowered energy yield. 



20 The network behavior was also examined for 

five disease cases (Table 6) . The test cases were 
chosen based on their physiological relevance to the 
model's predictive capabilities. In brief, McArdle's 
disease is marked by the impairment of glycogen 

25 breakdown. Tarui 1 s disease is characterized by a 
deficiency in phosphof ructokinase . The remaining 
diseases examined are marked by a deficiency of 
metabolic enzymes phosphoglycerate kinase, 
phosphoglycerate mutase, and lactate dehydrogenase. 

30 each case, the changes in flux and by-product secretion 
of metabolites were examined for an aerobic to 
anaerobic- metabolic shift with glycogen and 
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phosphocreatine as the sole carbon sources to the 
network and pyruvate, lactate, and albumin as the only 
metabolic by-products allowed to leave the system. To 
simulate the disease cases, the corresponding deficient 
5 enzyme was constrained to zero. In all cases, a severe 
reduction in energy production was demonstrated during 
exercise, representing the state of the disease as seen 
in clinical cases. 

Throughout this application various 
10 publications have been referenced. The disclosures of 
these publications in their entireties are hereby 
incorporated by reference in this application in order 
to more fully describe the state of the art to which 
this invention pertains . 

15 t Although the invention has been described 

with reference to the examples provided above, it 
should be understood that various modifications can be 
made without departing from the spirit of the 
invention. Accordingly, the invention is limited only 

20 by the claims. 
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Table 1 

Locus ID GeneAb. Reaction Stoichiometry E.C. 

1. Carbohydrate Metabolism 

1.1 Glycolysis / Gluconeogenesis [PATH:hsa00010] 



3098 HK1 


GLC + ATP -> GBP + ADP 


2.7.1.1 


3099 HK2 


GLC + ATP -> G6P + ADP 


2.7.1.1 


3101 HK3 


GLC + ATP ■> G6P + ADP 


2.7.1.1 


2645 GCK, HK4, MODY2, NIDDM 


GLC + ATP -> G6P + ADP 


.2.7.1.2 


2538 G6PC, G6PT 


G6P + H20 -> GLC + PI 


3,1.3.9 


2821 GPI 


G6P <■> F6P 


5.3.1.9 


5211 PFKL 


F6P + ATP -> FDP + ADp 


2.7.1.11 


5213 PFKM 


F6P + ATP -> FDP + ADP 


2.7.1.11 


5214 PFKP, PFK-C 


F6P + ATP -> FDP + ADP 


2 7 111 


5215 PFKX 


F6P + ATP -> FDP + ADP 


2 7 111 


2203 FBP1, FBP 


FDP + H20 -> F6P + PI 


3.1.3.11 


8789 FBP2 


FDP + H20 -> F6P + P) 


3.1.3.11 


226 ALDOA 


FDP <-> T3P2 + T3P1 


4 12 13 


229 ALDOB 


FDP <-> T3P2 + T3P1 


4 12 13 


230 ALDOC 


FDP <-> T3P2 + T3P1 


4 12 13 


7167 TPH 


T3P2 <-> T3P1 


5 3 11 


2597 GAPD GAPDH 


T3P1 + p| + WAD <-> NADH + 13PDG 


1 2 1 12 


26330 GAPDS GAPDH-2 


T3P1 + PI + NAD <-> NADH + 13PDG 


19 119 


5230 PGK1 PGKA 


13PDG + ADP <-> SPG + ATP 


2.7.2.3 


5233 PGK2 1 


1 3PDf5 + ADP <-> ^PR + ATP 


2.7.2.3 


5923 PRAM1 PfiAMA 


13PDG -> 23PDG 


5-4 ,2,4 






3.1.3.13 




oru ^- trvj 


BuSfafci I 


Pf3AM9 PHiAMM 




3.4.2.4 




oqpnfS + M9r» 'IPrs ■+■ PI 


0 4 0 40 
jj 1 .j. 1 3 




Or 0 ^ tru 


J.*f.^. 1 


fifiQ RPPM 
OPS? Dr Olvl 


loruo ^* toruo 


5.4^.4 






j. I . j, 1 0 




SPf^ <_> 9Pft 
oru — ^ ^rva 


*5 4 9 1 


9n93 PWPil PPM PMPill 1 


9Pf5 <-r> PPP + W90 
tro ^«-^ rcr T n^^j 




9n*?fi PMP9 


9Pf5 ppp 4. H9n 


A. 9 i 1 1 
f.Z.l . 1 1 


9tt97 Fwn^ 


9Pfi <-> PPP + H9H 
cry rcr ■ n^u 


4. 9 1 1 1 


9R937 PNP1R 


9Pfi <-> PPP 4- M9H 


A 9 1 11 


i«j 1 rvl— r\, r r\ 1 


PPP 4- AnP -> PYR 4- ATP 
rcr ~ nur r t r\ ~ t\ 1 r 


9 7 1 An 

£..1 .\ .<+U 


5315 PKM9 PK3 THRP1 DIP3 


PPP 4- AnP j> PYR + ATP 
rcr ■ rvur r t r\ ~ r\ 1 r 


9 7 1 An 


Mfifi PDHA1 PWP1A PnHA 


PYRm 4- PHAm 4- iSiAHm -> 4- KJAriHm 4- Pn9m 4- ArfDAm 


1 9 A 1 


51fi1 PDHA9 PnHAl 
P |P 1 run/vi, i UiIrVL. 


PYRm 4- mfim 4- KJADm _i> 4- KlADMrn 4- P09m 4- AfPO A.m 


1 9 A 1 


0 ipz runo 


PYRm 4- rriAm 4- MADm -> 4- KIADI-tm 4- r.D9m 4- AfiPPjAm 


19 4 1 


1737 ni at ni ta pnp-F9 


PYRm 4- rnAm 4- MADm -> 4- NADHm 4- Pn9m 4- APnDAm 


9 1 19 


ftn^n pnvi P3RP 


PYRm 4- rnAm 4- NADm ~> 4- MAnHm + P09m 4« APPOAm 


9 ? 1 19 


3Q3Q I'HHA I DM1 


NAD 4- 1 AP <-> PYR +.KJADH 


1 1 1 97 


3Q45 1 nHR 


NAD 4- 1 AP <-> PYR + MADH 


1 1 1 97 


3948 LDHC LDH3 


NAD + LAC <-> PYR + NADH 


1 1 1 27 


5236 PGM1 


G1P <-> G6P 


54 2 2 


5237TPGM2 


G1P <-> G6P 


5.4.2.2 


5238 PGM3 


G1P <-> G6P 


5.4.2.2 


1738 DLD LAD PHE3 DLDH E3 


DLIPOm + FADm <-> LIPOm + FADH2m 

VUll will » 1 #\W^4 1 1 ^ *^ wil III * ■ r^fc/l \C-\ 1 1 


1814 


124 ADH1 


ETH + NAD <-> ACAL + NADH 


1111 


125 ADH2 


ETH + NAD <-> ACAL + NADH 


1111 


1 26 ADH3 


ETH + NAD <-> ACAL + NADH 


1.1.1.1 


127 ADH4 


ETH + NAD <-> ACAL + NADH 


1.1.1.1 


12S ADH5 


FALD + RGT + NAD <-> FGT + NADH 


1.2.1.1 




ETH + NAD <-> ACAL + NADH 


1.1.1.1 


13Q ADH6 


ETH + NAD <-> ACAL + NADH 


1.1.1.1 


131 ADH7 


ETH + NAD <-> ACAL + NADH 


1.1.1.1 


1P327 AKR1A1, ALR, ALDR1 




1.1.1.2 


37. ACYP1 




3-6.1.7 


38 ACYP2 




3.6.1.7 


1.2 Citrate cycie (TCA cycle) PATH:hsaOO02O 






1431 CS 


ACCOAm + OArn + H20m -> COAm + CiTm 


4.1,3.7 


4S AC01, IREB1, IRP1 


CIT<->ICIT 


4.2.1.3 


5Q AC02 


CtTm <-> ICITm 


4.2.1.3 


MIL IDH1 


ICIT + NADP -> NADPH + C02 + AKG 


1.1.1.42 
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2418 1DH2 
MIS IDH3A 
242Q 1DH3B 
2421 IDH3G 
4967 OGDH 
1243 DLST, DLTS 
2202 SUCLG1, SUCUM 
BB03 SUCLA2 
22Z1 FH 

4130 MDH1 

4131 MDH2 
5031 PC, PCB 

4Z ACLY, ATPCU CLATP 
3657 

5105 PCK1 

51 06 PCK2, PEPCK 

1.3 Pentose phosphate cycle PATH:hsa00030 
2532 G6PD, G6PD1 
3562 H6PD 

252S6 PGLS, 6PGL 
522SPGD 
6120 RPE- 
7086 TKT 

8277 TKTL1 , TKR, TKT2 



6888 T ALP 01 
5621PRPS1,PRSI,PRS,I 
5534 PRPS2 t PRS II, PRS, II 
^£2GDH 

1.4 Pentose and glucuronate tnterconverslons PATH:hsa00040 
221 AKR1 B1 , AR, ALDR1 , ADR 

7359 UG P1 G1P + UTP -> UDPG + PPI 

7360 UGP2, UGPP2 G1P + UTP -> UDPG + PPI 
I25S UGDH, UDPGDH 

jnj2QUGT2B11 

54658 UGT1A1, UGT1A, GNT1, UGT1 

7361 UGT1A, UGT1, UGT1A 
1362 UGT2B, UGT2, UGT2B 

7363 UGT2B4, UGT2B11 

7364 UGT2B7, UGT2B9 
Z365UGT2B10 

7366 UGT2B15, UGT2B8 

7367 UGT2B17 
13AADAC, DAC 

2391 LIPE f LHS, HSL 
1.5 Fnjctose and mannose metabolism PATH:hsa00051 



ICITm + NADPm -> NADPHm + C02m + AKGm 

ICITm + NADm -> C02m + NADHm + AKGm 

ICITm + NADm -> C02m + NADHm + AKGm 

ICITm + NADm -> CQ2m + NADHm + AKGm 

AKGm + NADm + COAm -> C02m + NADHm + SUCCOAm 

AKGm + NADm + COAm -> C02m + NADHm + SUCCOAm 

GTPm + SUCCm + COAm <-> GDPm + Pirn + SUCCOAm 

ATPm + SUCCm + COAm <-> ADPm + Pirn + SUCCOAm 

FUMm + H20m <-> MALm 

MAL + NAD <-> NADH + OA 

MALm + NADm <-> NADHm + OAm 

PYRm + ATPm + C02m -> ADPm + OAm + Pirn 

ATP + C1T + COA + H20 -> ADP + PI + ACCOA + OA 

OA + GTP -> PEP + GDP + C02 

OAm + GTPm -> PEPm + GDPm + C02m 

G6P + NADP <-> D6PGL + NADPH 

D6PGL + H20 -> D6PGC 
D6PGL + H20->D6PGC 
D6PGC + NADP -> NADPH + C02 + RL5P 
RL5P.<->X5P 
R5P + X5P <-> T3P1 + S7P 
X5P + E4P <-> F6P + T3P1 
R5P + X5P <->T3P1 +S7P 
X5P + E4P <-> F6P + T3P1 
T3P1 + S7P <-> E4P + F6P 
R5P + ATP <-> PRPP + AMP 
R5P + ATP <-> PRPP + AMP 



4?51 MPI, PMH 
53Z2PMM1 

52Z3 PMM2, CDG1, CDGS 
2Z62GMDS 
£790 FPGT, GFPP 
52QZ PFKFB1, PFRX 

52QB PFKFB2 

5203 PFKFB3 

5210 PFKFB4 

2Z95KHK 
6252 SORD 

2525 FUT4, FCT3A, FUOTIV 
2523 FUT7 
3Q3S HAS1 , HAS 
2022 HAS2 



MAN6P <-> F6P 
MAN6P <->MAN1P 
MAN6P <-> MAN1P 



ATP + F6P -> ADP + F26P 

F26P -> F6P + PI 

ATP + F6P->ADP + F26P 

F26P -> F6P + PI 

ATP + F6P -> ADP + F26P 

F26P -> F6P + PI 

ATP + F6P->ADP + F26P 

F26P -> F6P + PI 

DSOT + NAD -> FRU + NADH 



1.1.1.42 
1.1.1.41 

1-1-1-41 
.1.1,1.41 
1.2.4.2 
2.3.1.61 

6.2.1.4 
6.2.1.4 
4.2.1.2 

1.1.1.37 
1-1.1-37 
6.4.1.1 
4.1,3.8 

4.1.1.32 
4.1.1.32 

1.1.1. 49 . 

1.1,1.47 

3.1.1.31 

3.1,1.31 

1.1.1.44 

5.1.3.1 

2,2,1. 1 

2.2.1.1 

2.2.1.2 
2.7.6.1 
2.7.6.1 
1-1.1.47 

1- 1-1.21 
27.7.9 
2.7.7.9 
.1.1.1.22 
2,4.1.17 
2.4.1.17 

2- 4.1.17 
2-4.1.17 
2,4.1,17 
2.4.1.17 
2.4.1.17 
2.4.1.17 
2.4.1.17 
3.1 .1 .- 
3.1. 1- 

5.3.1.8 

5.4.2.8 

5.4.2.8 

4,2.1.47 

2.7.7.30 

2.7.1 .105 

3 1.3.46 

2.7.1.105 

2X246 

2.7.1.105 

3.1.3.46 

2.7.1.105 

2J.3.46 

2.7.1.3 

1- 1-1.14 
2.4.1. - 
2A±i 
2.4.1. - 

2- 4.1.- 
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S4Z3 OGT, O-GLCNAC 
51144 LOC51 144 
1.6 Galactose metabolism PATH:hsa00052 

2534 GALK1, GALK GLAC + ATP -> GAL1P + ADP 

2535 GALK2, GK2 GLAC + ATP-> GAL1P + ADP 
25S2 GALT UTP + GAL1 P <-> PPI + UDPGAL 
2582 GALE UDPGAL <-> UDPG 
2Z2QGLB1 

3933 LOT, LAC 

_™ B4GALT1, GGTB2, BETA4GAL-T1 , 
^GTI.GTB 



29Q£ LALBA 
2Z12 GLA, GALA 
2545 GAA 

2594GANAB 

2525 GANC 

B222 MGAM, MG, MGA 



MELI -> GLC + GLAC 
MLT -> 2 GLC 
6DGLC -> GLAC + GLC 
MLT -> 2 GLC 
6DGLC -> GLAC + GLC 
MLT -> 2 GLC 
6DGLC-> GLAC + GLC 
MLT -> 2 GLC 
6DGLC-> GLAC + GLC 



1.7 Ascorbate and eldarate metabolism PATH:hsa00053 

213 ALDH1, PUMB1 ACAL + NAD -> NADH + AC 

2UALDH2 ACALm + NADm -> NADHm + ACm 

219 ALDH5, ALDHX 
223ALDH9, E3 



224 ALDH10, FALDH, SLS 
3354RALDH2 

1591 CYP24 

1592 CYP26A1, P450RA1 

1593 CYP27A1, CTX, CYP27 
CYP27B1, PDDR, VDD1, VDR, CYP1, 

*VDDR, I.P450C1 
1.8 Pyarvate metabolism PATH:hsaOO620 
54988 FU20581 

21 ACACA,ACAC,ACC 



1594 



ATP + AC + COA -> AMP + PPI + ACCOA 
ACCOA + ATP + C02 <-> MALCOA + ADP + PI + H 



32 ACACB, ACCB, HACC275, ACC2 ACCOA + ATP + C02 <-> MALCOA + ADP + PI + H 



2739 GL01, GLYI 
3029 HAGH, GL02 
2223 FDH 

938Q GRHPR, GLXR 
42QQME2 
133Z3ME3 
29S9Z HUMNDME 
4199 ME1 

3fi ACAT1, ACAT, T2. THIL, MAT 
39ACAT2 

1.9 Glyoxylate and dicarboxylate metabolism PATH:hsa00630 
524QPGP 
2253 GLYD 

10797 MTHFD2, NMDMC 



RGT + MTHGXL <-> LGT 
LGT -> RGT + LAC 
FALD + RGT + NAD <-> FGT + NADH 

MALm + NADm -> C02m + NADHm + PYRm 
MALm + NADPm -> C02m + NADPHm + PYRm 
MAL + NADP -> C02 + NADPH + PYR 
MAL + NADP -> C02 + NADPH + PYR 
2 ACCOAm <-> COAm + AACCOAm 
2 ACCOAm <-> COAm + AACCOAm 



4522MTHFD1 



1.10 Propanoate metabolism PATH:hsa00640 
34 ACADM, MCAD 



35 ACADSB 



3HPm + NADHm -> NADm + GLYAm 
METHF <-> FTHF 
METTHF + NAD -> METHF + NADH 
METTHF + NADP <-> METHF + NADPH 
METHF <-> FTHF 

THF + FOR + ATP-> ADP + PI + FTHF 

MBCOAm + FADm -> MCCOAm + FADH2m 
iBCOAm + FADm -> MACOAm + FADH2m 
IVCOAm + FADm -> MCRCOAm + FADH2m 
MBCOAm + FADm -> MCCOAm + FADH2m 



2.4,1. - 
.1.1.1.- 

2.7.1.6 
2,7,1 A 
2.7.7 .1 Q 
5.1.3.2 
3^123. 

3.2.1.108 

2.4.1.90 

2.4.1.38 
2.4.1.22 
2AA2Z 
3. 2,1. 22 
3.2.1.20 

3.2.1.20 

3.2.1.20 

3.2.1.20 

3.2.1.3 

1.2.1.3 

1,2,1.3 

1.2.1.3 

1.2.1.3 

.1.2.1.19 

1.2.1.3 

1.2.1.3 

1,14.-,- 

.1,1 4.-.- 

.1.1 4.-.- 

1.14.-,- 



6.2.1.1 

6.4.1.2 

6,3.4.14 

6.4.1.2 

6.3.4.14 

4.4.1.5 

3.1.2.6 

1.2.1.1 

1.1.1.79 

1X133 

1.1.1.40 

. 1 .1,1.40 
.1.1,140 
2.3.1.9 
2.3.1.9 

3, 1.3.18 

1, 1.1.29 

2<5A9_ 

1.5.1.15 

1.5.1.15 

3.5.4.9 

6.3.4.3 

1.3.99.3 



1.3.99.3 
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iBS2ECHS1,SCEH 
1962 EHHADH 

3030 HADHA, MTPA, GBP 



23417 MLYCD. MCD 
IB ABAT, GABAT 
5Q95PCCA 
5QSSPCCB 
4594 MUT, MCM 
4229 MMSDH 

3523 FACVL1, VLCS, VLACS 
1.11 Butarroate metabolism PATH:hsa00650 

3Q23 HADH2, ERAB 

3032 HADHSC, SCHAD 
35ACADS,SCAD 

7915 ALDH5A1, SSADH, SSDH 
2571 GAD1, GAD, GAD67, GAD25 
25Z2GAD2 
25Z3GAD3 

315ZHMGCS1, HMGCS 
3158 HMGCS2 
3155 HMGCU HL 
5Q1SOXCT 
S22 BDH 

1629 DBT, BCATE2 



1.13 Inositol metabolism PATH:hsa00031 
Z Energy Metabolism 

2.1 Oxidative phosphorylation PATH:hsa00190 
4535 MTND1 
4536MTND2 

4537 MTND3 

4538 MTND4 

4539 MTND4L 

4540 MTNP5 

4541 MTND6 
4SaiNDUFA1 t MWFE 

4695 NDUFA2, B8 

4696 NDUFA3, B9 

4697 NDUFA4, MLRQ 

4333 NDUFA5, UQOR13, B13 

4700 NDUFA6, B14 

4ZQ1 NDUFA7, 614.5a, B14.5A 

4ZQ2 NDUFA8, PGIV 

4104 NDUFA9, NDUFS2L 

4ZC5 NDUFA10 



IBCOAm + FADm -> MACOAm + FADH2m 
IVCOAm + FADm -> MCRCOAm + FADH2m 
MACOAm + H20m -> HIBCOAm 
MCCOAm + H20m -> MHVCOAm 
MHVCOAm + NADm -> MAACOAm + NADHm 
HIBm + NADm -> MMAm + NADHm 
MACOAm + H20m -> HIBCOAm 
MCCOAm + H20m -> MHVCOAm 
MHVCOAm + NADm -> MAACOAm + NADHm 
HIBm + NADm -> MMAm + NADHm 
MACOAm + H20m -> HIBCOAm 
MCCOAm + H20m -> MHVCOAm 

' C160CARm + COAm + FADm + NADm -> FADH2m + NADHm + 
C140COAm+ ACCOAm 

GABA + AKG -> SUCCSAL + GLU 

PROPCOAm + C02m + ATPm -> ADPm + Pirn + DMMCOAm 
PROPCOAm + C02m + ATPm -> ADPm + Plm + DMMCOAm 
LMMCOAm -> SUCCOAm 

MMAm + COAm + NADm -> NADHm + C02m + PROPCOAm 



C140COAm + 7 COAm + 7 FADm + 7 NADm -> 7 FADH2m + 7 
NADHm + 7 ACCOAm 

MBCOAm + FADm -> MCCOAm + FADH2m 
IBCOAm + FADm -> MACOAm + FADH2m 

GLU -> GABA + C02 

GLU -> GABA + C02 

GLU -> GABA + C02 

H3MCOA + COA <-> ACCOA + AACCOA 

H3MCOA + COA <-> ACCOA + AACCOA 

H3MCOAm -> ACCOAm + ACTACm 

3HBm + NADm NADHm + Hm + ACTACm 
OMVALm + COAm + NADm -> MBCOAm + NADHm + C02m 
OIVALm + COAm + NADm -> IBCOAm + NADHm + C02m 
OICAPm + COAm + NADHm -> IVCOAm + NADHm + C02m 



NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm +4 Hm-> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm-> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm +4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H . 



4,2 .1 , 17 
1.1.1.35 
4,2,1.17 
1.1.1.35 
4-2.1.17 

1.1.1-35 

4.2.1.17 

4.1.1.9 

2.6.1.19 

6.4.1.3 

6.4.1.3 

5.4.99.2 

■U.1.27 

6.2. 1.- 

.1.1.1.35 

1.1.1.35 
1.3.99 ,2 

1.2.1.24 

4.1.1.15 

4 . 1 .1. 15 - 

4 .1.1. 15 

4.1.3.5 

4.1.3.5 

4.1,3.4 

2.8.3.5. 

1,1.1,30 

2.3.1. - 



1.6.5.3 

1.6.5.3 

1.6.5.3 

1.6.5.3 

1.6.5.3 

1.6.5.3 

1.6.5.3 

1.6.5.3 

1.6.5.3 

1.6.99.3 

1.6.5.3 

1.6.99.3 

J. 6,5,3 

1.6.99.3 

1.6.5.3 

1.6.99.3 

1.6.5.3 

1.6.99.3 

1.6.5.3 

1.6.99.3 

1.6.5.3 

1-6.99.3 

1.6.5.3 

1.6.99.3 

1,6.5.3 
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MADHm 4- Om + 4 Hm -> QH2m + NADm + 4 H 

|X/*\^f till * 1 1 * *T % llll \*m* ■ i " l^#»*^lll ' • M 


1 ,6.99.3 


47ft.fi Mnl 1FAP.1 SD AP 


NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 


1 .6.5.3 




MADHm + Om + 4 Hm •> OH 2m + NADm + 4 H 

I V/AL^ I 1 1 1 1 * II T ^ 1 111* ^**l * 1 ^Ji4** ■ 1 1 'Til 




4707 WDUFB1 MNI L CI-SGDH 


NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 

1 >• r^i_/ 1 llll ■ Will • T 1 Mil r \*%* II 1 ~ % * 


1.6.5.3 




NADHm + Om + 4 Hm -> QH2m + NADm + 4 H 


^ «6 .99*3 


A7ftA Kim IFR9 Afiflfi 


NADHm 4- Om 4- 4 Hm -> QH2m 4- NADm + 4 H 

linuni 1 1 * veil J ■ ^ 1 nil ^ 1 ifcin * i^#w^iii ■■tii 


1.6.5.3 




MAnHm +• Om 4- 4 Hm -> OM?m 4- NADm 4- 4 H 


1.6.99.3 


A7flQ Mnl IPR^ RIO 

4/ua iNUUrtJo, Diz 


NADHm 4- Om 4- 4 Hm -> OH2m 4- NADm 4- 4 H 

IHrALJ n 1 1 1 ■ Vail 11 ' " lUII ^ Ml Ifctll • Iw^WIII T ~ 1 1 


1 .6.5.3 




MAnHm 4- Om 4- 4 Hm OH?m 4- NADm 4- 4 H 


1.6.99.3 


/74n MniiFRA ric 

4/ ID NUUPD4, DlO 


MADHm 4- Om + 4 Hm -> DH7ni 4- NADm + 4 H 


1.6.5.3 




KiAHMm 4. Dm + 4 Hm -> DM9m 4* MADm + 4 H 


1 6 99 3 


^*74<t Kim i cdc crnu 
4711 NUUrKJ, ooUh 


MAnHm 4- Dm 4- 4 Hm -t> nH9m 4- NADm 4- 4 H 


•j.6.5.3 




KIAnUm a. Om 4- j4 Urn OW9m 4> MADm 4- 4 H 


1 6 99 3 


471 4 NDUrBo, B17 


MAnHm 4. Om 4- 4 Hm flH9m 4- MADm 4- 4 H 


1.6.5.3 




MfinUm Am Om 4- A Mm OWOm 4- KlADm + 4 H 

NAUrtui + um + *t nm -> wnzm + NMum t^n 


1 A QQ O 

. 1 .U.JJ. kJ 


i^TJA I f rM IPB7 Bale 

4713 NDUFB7, BIB 


iMAUrtm + uim ■ *t nm un^rn t inmutti ~ h n 


1.6.5.3 






1 .6.99.3 


4714 NDUFBo, ASHi 


MAnUm -l. Om -J- A Um DW9m ■#- MAnm 4- 4 H 


1 .6.5.3 




klAnUm a. Om 4. A Um -!% OMOm 4- WAHm 4- A H 


1 .0. j j.j 


A^Atr tini icon ■ i^nDHO nn 

4715 NDUFB9, UQORzz, Bzz 


Kl A nUm A- Om 4- A Mm nH9m J. MAnm 4- A M 


16 5 3 




MAnUm j_ Om J. >1 Mm *> OMOm 4- MAHm 4> A M 

NAUnm 4- tjm * 4 nm un/m t oiMum t4 »n 


1 -D.33.iJ 


4-716 NDUFB10, PD5W 


MfinUm a. Om 4* A Mm ^> OH9m 4- MAnm + 4 H 


1 6.5.3 




MAIMJm j_ Om J. jl Um >% OUOm 4- KlADm 4. A U 

NAUnm 4- um + 4 nm -> un/m + NAum + 4 n 


1 CQQO 

I ,D.3J.J 


4717 NDUFC1 , KFYI 


MArtUm A. Om 4. A Um "> OUOm 4- KJAHm 4> A M 

NAUnm 4- vjm + 4 nm -> Un/m t w\Um f 4 n 


1.6.5.3 




klArsUm J. Om U. X Um S. OUOm A. KlfiHm 4- A U 

NAL/rim 4- um + 4 nm -> Un^m + iM/u/m + <(n 


1.6.99.3 


4718 NDUFC2, B 14.5b, B14.5B 


KlAnUm j_ Om 4- A Um OUOm 4- MAnm 4> A M 

NAUnm 4- um + 4 nm -> un^m + iMMum t ^n 


1 6.5.3 




KlAnUm J. Om 4. A Um »» OUOm 4- MAnm 4- A U 

NAUnm 4- um * 4 nm -> un/m ▼ wMJm t< tn 


| .O.jj.ij 


4724 NDUFS4, AQDQ 


NADHm 4- um + 4 Hm -> QHzm 4- NAum + 4n 


^ ceo 
T .P.P. J 




MAfMJm _l_ Om A. A Um OUOm 4- MAnm 4> A U 

NAUnm 4- um + 4- nm -> unzm + wum t- n 


1 .D.3J.O 


4725 NDUFS5 


NADHm 4- Qm + 4 Hm -> unzm + NAum + 4 H 


1 .13. D.J 




KIAHLiiM j_ Om J. it Um ** OUOm X MAnm 4. A U 

NAUHm 4- um + 4 rim -> unzm * iMnum + 4 n 


4 fi QQ "3 


4726 NDUFS6 


kl »r\U« ■ /»^__ ■ * , «^ OUO»w» x MAnm x A U 

NADHm + am 4- 4 Hm -> unzm + NAum *■ 4 n 


4 c c 0 




NADHm + Qm + 4 Hm -> QHzm 4- NADm + 4H 




4731 NDUFV3 


KlAHUm n«» j_ >i 1 I— w /^LJO—, « KIAnm A. A U 

NADHm 4- um + 4 Hm -> unzm + NAUm + 4 n 


,1 .P.O. J 


4727 NDUFS7, PSST 


NADHm 4- Qm 4- 4 Hm -> QH2m 4- NADm + 4rl 


4 C C Q 

1 .P.P. J 




MArMJm t_ Or— 1 _l >t Um ««. OUOm J- MAnm X A U 

NADHm + um + 4 Hm -> unzm + NAum + 4 n 


4 C QQ Q 


4722 NDUFS3 


NADHm 4- Qm + 4 Hm -> QH2m 4- NADm + 4 H 


1 .P.D.J 




klAnUm J- Om -L /I Um *^ OUOm JL MAnm X A U 

NAUHm + Qm + 4. nm -> QHzm + NAum + 4n 


1 K QQ *3 


4720 NDUFS2 


NADHm 4- Qm + 4 Hm -> QH2m 4- NADm + 4 n 


4 Ct K. Q 


4729 NDUFV2 


MAHUm X Om X A 14m -s. OUOm X MAnm X A U 

NAUHm 4- um + 4 Hm *■> unzm + NAum + 4n 


I.P.D.O 




NADHm + Qm + 4 Hm -> QHzm 4- NAum + 4 n 


4 R QQ 9 


4723 NDUFV1 , UQOR1 


MAHUm X Om X / Um *» OUOm X MAnm X A U 

NAUHm 4" um + 4 Hm -> unzm + NAum + 4n 


i ceo 




NADHm + Qm + 4 Hm -> QHzm 4- NADm 4- 4 H 


4 A QQ 'S 


4719 NDUFS1, PRO1304 


MAnUm x Om X A Um *> OUOm X MAnm X A U 

NAUHm + Qm + 4 Hm -> unzm iNAum + in 


4 fi QQ *a 

1 .D.Jj.o 




KlAnUm -t- Aw« _i_ A Um nUOm x MAnm X A 14 

NADHm 4- um + 4 Hm •> QHzm + NAum + 4 n 


< c c 0 
1 r P.O.O 


4728 NDUFS8 


klAHUm X Om X A Um «w OLIOm X KIAnm X A U 

NADHm + Qm 4* 4 Hm -> QHZm NAUm + 4 n 


i ceo 

1 -O.O.J 




KlAnUm X Om X A Um «». OUOm X MAnm X A U 

NADHm 4- Qm + 4 Hm -> QHzm + NAum + 4 n 


4 C QQ O 

1 .0.ag.J 


6391 SDHC 


CM l/^Om X CAOm ^ s CI Uim X C AnUOm 

oUOum + rAum <~> ruwim + rAUnzm 


1 .3.5.1 




CAnUOm X Om ^ -s. CAHm x OUOm 

rAUnzm + um <-> rftum + unzm 




£232 SDHD, CBT1, PGL, PGL1 


Of lOOm X CAn~ ^ «w CI IKJm X CAnUOm 

oUULrm + rAUm <-> rUMm rAunzm 


1 .3.5.1 




1 PAnHOm + Om <"^> Ffinm X OHOm 

1 rMun^in t um rnum f unzm 




fiSBS SDHA, SDHz, SDHF, FP 


CI lOOm X Cflrim ^ v C| iy m a- C&nMOm 

aUUum •+ rAum <•■> ruivim *r rAunzm 


1^3 5 1 




CAnUOm X Om ^ CAhm X OUOm 

rAUnzm um <-> rflum + unzm 




coon crMJQ or\un in onu 
fi22Q bUHB, 5DH1, IP, SDH 


Q1 lOOm 4- FADm <_> Fl llUlm + FAnHOm 


1 .3.5.1 




FfinW9m X Om «r_^ CfiHm X OMOm 

1 nun&in ~ win rttuiii ~ unz.111 




•7150c 1 mrtjcc t Dic-i 
73oo UQORroi, KI51 


OOm 4- 4 FFROm 4- 4 Hm _> 4 FFRlm 4- 7 H?Om + 4 H 


1.10.2.2 


ftplg M 1 U YD 


OOm + 4 FFROm 4- 4 Hm -> 4 FFRlm 4- 0 H20m 4- 4 H 


1.10.2.2 


1537 CYO 1 


OOm 4* 4 FFROm x A Mm A PPPlm X O HOOm 4- 4 H 

uzm ^ 4 rcnuin * 4 nm *f rtrvim ▼ x n&um ~ *+ n 


^.10.2.2 




02m + 4 FEROm + 4 Hm -> 4 FERIm + 2 H20m + 4 H 


j.10,2.2 


,7385 UQCRC2 


02m + 4 FEROm + 4 Hm -> 4 FERIm + 2 H20m 4- 4 H 


1 r 10.2.? 


7388 UQCRH 


02m + 4 FEROm 4- 4 Hm -> 4 FERIm 4- 2 H20m 4- 4 H 


g t .1P,2,2 


iaai UQCRB, QPC, UQBP, QP-C 


02m + 4 FEROm + 4Hm->4 FERIm + 2 H20m + 4 H 


1.10.2^ 


2ZQfiaQP-C 


02m + 4 FEROm + 4 Hm -> 4 FERIm 4- 2 H20m + 4 H 


1.10.^2 


iflSZSUQCR 


02m + 4 FEROm + 4 Hm -> 4 FERIm + 2 H20m + 4 H 


1.10.2,2 


1333 COX5BL4 


QH2m 4- 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 


1-9.3.1 


4514 MTC03 


QH2m + 2 FERIm 4- 4 Hm -> Qm 4- 2 FEROm + 4 H 


j.9.3.1 


4512 MTC01 


QH2m 4- 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 


1.9.3.1 
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4513 MTC02 
1329COX5B 
122ZCOX4 

J332COX6A1.COX6A 

1339 COX6A2 

1340 COX6B 
1345COX6C 

9377 COX5A, COX, VA, COX-VA 

1346 COX7A1, COX7AM, COX7A 

1347 COX7A2, COXVlla-L 
1348CQX7A3 

1349 COX7B 

9167 COX7A2U COX7RP, EB1 

1350 COX7C 
J1351COX8. COX VIII 

4508 MTATP6 

4509 MTATP8 
499 ATP5A2 

507 ATP5BL1. ATPSBL1 

508 ATP5BL2, ATPSBL2 
SIS ATP5H 

532 ATP6S1, ORF, VATPS1, XAP-3 
514ATP5E 
513ATP5D ' 
5Qfi ATP5B, ATPSB 
* 5Q2ATP5C1, ATP5C 



QH2m + 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 
QH2m + 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 
QH2m + 2 FERIm + 4 Hm •> Qm + 2 FEROm + 4 H 
QH2m + 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 
QH2m + 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 
QH2m + 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 
QH2m + 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 
QH2m + 2 FERIm + 4 Hm •> Qm + 2 FEROm + 4 H 
QH2m + 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 
QH2m + 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 
QH2m + 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 
QH2m + 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 
QH2m + 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 
QH2m + 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 
QH2m + 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 
ADPm + Plm + 3H -> ATPm + 3 Hm + H20m 
ADPm + Pim + 3H -> ATPm + 3 Hm + H20m 
ADPm + Pirn + 3 H -> ATPm + 3 Hm + H20m 
ADPm + Pirn + 3 H -> ATPm + 3 Hm + H20m 
ADPm + Pim + 3H -> ATPm + 3 Hm + H20m 
ADPm + Pim + 3H -> ATPm + 3 Hm + H20m 
ADPm + Pim + 3H -> ATPm + 3 Hm + H20m 
ADPm + Pirn + 3 H -> ATPm + 3 Hm + H20m 
ADPm + Pirn + 3 H -> ATPm + 3 Hm + H20m 
ADPm + Pirn + 3 H -> ATPm + 3 Hm + H20m 
ADPm + Pim + 3H -> ATPm + 3 Hm + H20m 



498 ATP5A1, ATP5A, ATPM, OMR, HATP1 ADPm + Pirn + 3 H -> ATPm + 3 Hm + H20m 



539 ATP50, ATPO, OSCP ADPm + Pirn + 3 H 

516 ATP5G1 , ATP5G ADPm + Pirn + 3 H 

517 ATP5G2 ADPm + Pirn + 3 H 

518 ATP5G3 ADPm + Pirn + 3 H 
515 ATP5F1 ADPm + Pirn + 3 H 

521 ATP5I ADPm + Pirn + 3 H 

522 ATP5J, ATP5A, ATPM, ATP5 ADPm + Pirn + 3 H 
9551 ATP5J2, ATP5JL, F1FO-ATPASE ADPm + Pirn + 3 H 

10476 ATP5JD ADPm + Pim + 3H 

10632 ATP5JG ADPm + Pirn + 3 H 

223£ ATP6S14 ADPm + Pirn + 3 H 

528 ATP6D ADPm + Pirn + 3 H 

523 ATP6A1, VPP2 ADPm + Pirn + 3 H 
524ATP6A2.VPP2 ADPm + Pirn + 3H 
525 ATP6B1 . VPP3, VATB ADPm + Pirn + 3 H 
525 ATP6B2, VPP3 ADPm + Pirn + 3 H 
529ATP6E ADPm + Pim + 3H 
527 ATP6C, ATPL ADPm + Plm + 3 H 
523ATP6F ADPm + Pirn + 3 H 



-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 
-> ATPm 



+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 
+ 3Hm 



+ H20m 
+ H20m 
+ H20m 
+ H20m 
+ H20m 
+ H20m 
+ H20m 
+ H20m 
+ H20m 
+ H20m 
+ H20m 
+ H20m 
+ H20m 
+ H20m 
+ H20m 
+ H20m 
+ H20m 
+ H20m 
+ H20m 



inqi2 TCIRG1,TIRC7,00116,OC116kDa, ^p^p^^H -> ATPm + 3 Hm + H20m 
2j ^ u 0C-116KDA, ATP6N1C 



23545 TJ6 
5Q61Z ATP6N1B 

535ATP6N1 
51382 VATD 
8992 ATP6H 
9550 ATP6J 
51606 LOC51606 

4S5ATP4A, ATP6A 

495ATP4B, ATP6B 

4ZBATP1A1 

477ATP1A2 

478 ATP1A3 

479 ATP1AL1 
23432ATP1B4 

451 ATP1 B 1 , ATP1 B 
482ATP1B2. AMOG 
483ATP1B3 
2ZQ32 ATP2C1, ATP2C1A, PMR1 



ADPm + Pirn + 3 H 
ADPm + Pim + 3H 
ADPm +.Pim + 3H 
ADPm + Pim + 3H 



ATPm + 3 Hm + H20m 
ATPm + 3 Hm + H20m 
ATPm + 3 Hm + H20m 
ATPm + 3 Hm + H20m 



ADPm + Pirn + 3 H -> ATPm + 3 Hm + H20m 
ADPm + Pirn + 3 H -> ATPm + 3 Hm + H20m 
ADPm + Pirn + 3 H -> ATPm + 3 Hm + H20m 
ATP + H + Kxt + H20 <-> ADP + PI + Hext + K 
ATP + H + Kxt + H20 <-> ADP + PI + Hext + K 
ATP + 3 NA + 2 Kxt + H20 <-> ADP + 3 NAxl + 
ATP + 3 NA + 2 Kxt + H20 <-> ADP + 3 NAxl + 
ATP + 3 NA + 2 Kxt + H20 <-> ADP + 3 NAxt + 
ATP + 3 NA + 2 Kxt + H 20 <-> AD P + 3 NAxt + 
ATP + 3NA + 2Kxt+H20 <-> ADP + 3 MAxt + 
ATP + 3 NA + 2 Kxt + H20 <-> AD P + 3 NAxt + 
ATP + 3 NA + 2 Kxt + H20 <-> ADP + 3 NAxt + 
ATP + 3NA + 2Kxt+H20 <-> ADP + 3 NAxt + 
ATP+2CA + H20<-> ADP + PI +2CAxt 



2K + PI 
2K + PI 
2K + PI 
2K + PI 
2K + PI 
2K + PI 
2K + PI 
2K + PI 



1.9.3.1 
1.9.3.1 
19.3.1 

1 ,9,3. 1 

1.9.3.1 
1.9.3.1 
1.9.3.1 

1.9.3.1 

1.9.3.1 

,1.9.3.1 

1.9.3.1 

1.9.3.1 

1.9.3.1 

1.9.3.1 

1.9,3.1 

3.6.1.34 

3,6.1. 34 

3.6.1.34 

3.6.1.34 

3.6.1.34 

3.6.1.34 

3.6.1.34 

3.6.1.34 

3-6.1.34 

3.6.1.34 

3,6,1.34 

3,6.1.34 

3,6.1.34. 
3,6.1 .31 
3,6.1.34 

3.6.1.34 

3.6.1.34 

3.6.1.34 

3.6.1.34 

3,6.1.34. 

3,6.1, 34 

3 ,6 .1 .34 

3.6.1,34- 

3.6.1.34 

3.6.1. 34 

3.6.1.34 

36-1.34 

36.1.34 

3.6.1.34 

3-6.1.34 

3.6.1.34 



3.6.1.34 

3,6.1.34 
3R.1.34 
3.6.1.34 

3.5JLM 
3R1.34 
a.fi.i.34 

3R.1.34 
3 R 1.36 
3 R 1.36 
3.R.1.37 
3R1.37 
3R1.37 
3R.1.37 
3 6.1.37 
3 6.1.37 
3.6.1.37 

3.6.1.37 
3.6.1 .38_ 
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487 ATP2A1 , SERCA1 , ATP2A 

488 ATP2A2, ATP2B, SERCA2, DAR, DD 

489 ATP2A3, SERCA3 
42Q ATP2B1 , PMCA1 
491 ATP2B2, PMCA2 
492ATP2B3, PMCA3 
493 ATP2B4, ATP2B2, PMCA4 
538 ATP7A, MK, MNK, OHS 
540 ATP7B, WND 

5464 PP, SID6-8061 

2.2 Photosynthesis PATH:hsaO0195 

2.3 Carbon fixation PATH:hsa00710 

26Q5GOT1 
2BQ6GOT2 
2SZ5GPT 

2.4 Reductive carboxyiate cycle (C02 fixation) PATH:hsa00720 

2.5 Methane metabolism PATH:hsa00680 

B4Z CAT 
4025 LPO, SPO 
4352 MPO 

£2fia EPX, EPX-PEN, EPO, EPP 
S5Ba KIAA0106, AOP2 
6470 SHMT1. CSHMT 
6472 SHMT2, GLYA, SHMT 
51 004 LOC51004 



ATP + 2 CA + H20 <-> ADP + PI + 2 CAxt 
ATP + 2 CA + H20 <-> ADP + PI + 2 CAxt 
ATP + 2 CA + H20 <-> ADP + PI + 2 CAxt 
ATP + 2 CA + H20 <-> ADP + PI + 2 CAxt 
ATP + 2 CA + H20 <-> ADP + PI + 2 CAxt 
ATP + 2 CA + H20 <-> ADP + PI + 2 CAxt 
ATP + 2 CA + H20 <-> ADP + PI + 2 CAxt 
ATP + H20 + Cu2 -> ADP + PI + Cu2xt 
ATP + H20 + Cu2 -> ADP + Pl + Cu2xt 
PPI->2PI 



OAm + GLUm <-> ASPm + AKGm 
OA + GLU <-> ASP + AKG 
PYR + GLU <-> AKG + ALA 



2H2G2->02 



9420 CYP7B1 

2.6 Nitrogen metabolism PATH:hsa00910 
11238 CA5B 
23532 CA14 
259 CA1 

760 CA2 

761 CA3, CAlll 
2fi2CA4,CA1V 

763 CA5A, CAS, CAV, CAVA 

765 CAS 

ZSSCA7 

ZfiZCA8,CALS,CARP 
768 CA9, MN 
2ZD CA1 1 , CARP2 
771 CA12 
1373 CPS1 

275 AMT 

2Q24HAL,HSTD t HIS 
2Z45 GLUD1, GLUD 

8307 GLUD2 

2Z52 GLUL, GLNS 
22S42 KIAA0838 
2H£5GA 

2Z44GLS 
44Q ASNS 

1491 CTH 

2.7 Sulfur metabolism PATH:hsa00920 
9060 PAPSS2, ATPSK2, SK2 

3061 PAPSS1, ATPSK1, SK1 

1Q2BQ BPNT1 
6799 SULT1A2 
fifijJSULT1A1 t STP1 
6818 SULT1A3, STM 
££22 SULT2A1 , STD 



THF + SER <-> GLY + METTRF 
THFm + SERm <-> GLYm + METTHFm 
20PMPm + 02m -> 20PMBm 
20PMMBm + 02m -> 20MHMBm 
20PMPm + 02m -> 20PMBm 
20PMMBm + 02m -> 20MHMBm 



GLUm + C02m + 2 ATPm -> 2 ADPm + 2 Plm + CAPm 
GLYm + THFm + NADm <-> METTHFm + NADHm + CO 2m + 
NH3m 

HIS->NH3 + URO 

AKGm + NADHm + NH3m <-> NADm + H20m + GLUm 
AKGm + NADPHm + NH3m <-> NADPm + H20m + GLUm 
AKGm + NADHm + NH3m <-> NADm + H20m + GLUm 
AKGm + NADPHm + NH3m <-> NADPm + H20m + GLUm 
GLUm + NH3m + ATPm -> GLNrn + ADPm + Pirn 
GUM -> GLU + NH3 
GLN ■> GLU + NH3 
GLNrn -> GLUm + NH3m 

ASPm + ATPm + GLNrn -> GLUm + ASNm + AMPm + PPIm 
LLCT + H20 -> CYS + HSER 
OBUT + NH3 <-> HSER 

APS + ATP -> ADP + PAPS 
SLF+ATP->PPI + APS 
APS + ATP -> ADP + PAPS 
SLF + ATP -> PPI + APS 
PAP -> AMP + PI 



3.6.1.38 

3.6.1.38 

3.6,1. 33 

3.6.1.38 

3.6.1.38 

3.6.1.38 

3,6,1.38 

3.6.3.4 

3-6.3.4 

3.6.1.1 



2,6,1.1 

2.6.1.1 
2 6.1.2 



1.11.1.6 

1.11.1.7 

1.11.1.7 

1.11.1.7 

111.1.7 

2.1.2.1 

2.1.2.1 

1.14,13,- 

1 ,14 .13.- 



4.2.1.1 
4.2.1.1 
4.2.1.1 
4.2,1,1 
4.2.1.1 
4 .2.1. 1 
4,2,1.1 
4,2.1.1 
4-2.1.1 

4.2.1.1 

42.1.1 
4.2.1.1 

4.2.1.1 

5.3. 4 . 16 

2.1.2.10 

4.3.1,3. 

1.4.1.3 

1.4.1.3. 

6 3.1.2 
&5X2 
3.5X2 
&5X2 
B.3.5.4 
4,4.1.1 
4 4.1.1 

2 7.1.25 

27.7.4 

2.7.1.25 

2.7.7.4 

31-3.7 

2 r B.2.1 

2 , 8.2.1 

2-8.2.1 

Z&2JI 
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6783 STE, EST 
£B21SUOX 
3. Lipid Metabolism 

3.1 Fatty acid biosynthesis (path 1) PATH:hsa00061 

2334 FASN 

3.2 Fatty acid biosynthesis (path 2) PATH:hsa00062 

10449 ACAA2, DSAEC MAACOAm -> ACCOAm + PROPCOAm 

3Q ACAA1 , ACAA MAACOA -> ACCOA + PROPCOA 

2Q32 HADHB MAACOA -> ACCOA + PROPCOA 

3.3 Fatty acid metabolism PATH :hsa00071 

51AC0X1.AC0X 
32ACADU LCAD ' 
2639 GCDH 

2179 FACL1 f LACS ATP + LCCA + COA <-> AMP + PPI + ACOA 

21SQ FACL2, FACL1, LACS2 ATP + LCCA + COA <-> AMP + PPI + ACOA 

21BZ FACL4, ACS4 ATP + LCCA + COA <-> AMP + PPI + ACOA 

ISM CPT1A, CPT1, CPT1-L 
12Z5 CPT1B, CPT1-M 
1376 CPT2, CPT1, CPTASE 
ifi22 DCI 
11283 CYP4F8 
1M3CYP1A1, CYP1 

1544 CYP1A2 

1545 CYP1B1, GLC3A 
1548 CYP2A6, CYP2A3 
1543CYP2A7 
1551CYP3A7 
1553CYP2A13 

1554 CYP2B 

1555 CYP2B6 

1557 CYP2C19, CYP2C, P450UC19 
155BCYP2CB 

1559 CYP2C9, P450UC9, CYP2C10 
1562 CYP2C18, P450IIC17, CYP2C17 
1565 CYP2D6 

1571 CYP2E, CYP2E1, P450C2E 

1572 CYP2F1, CYP2F 
15Z3CYP2J2 

1575 CYP3A3 
15Z6CYP3A4 
1577 CYP3A5, PCN3 
1580 CYP4B1 
15&8CYP19, ARO 
J595CYP51 
134AHHR, AHH 

3.4 Synthesis and degradation of ketone bodies PATH:hsa00072 

3.5 Sterol biosynthesis PATH:hsa00100 



3156 HMGCR 
4598 MVH MVLK 



1Q654 PMVK, PMKASE, PMK, HUMPMKI 
4597 MVD, MPD 
3422 IDI1 
2224 FDPS 

2453 GGPS1, GGPPS 



2222 FDFT1, DGPT 
6Z13SQLE 
4047 LSS, OSC 

iZ2S DIA4, NMOR1, NQ01, NMORl 
4S35 NMOR2, NQ02 

3Z ACADVL, VLCAD, LCACD 
3.6 Bile acid biosynthesis PATH:hsa00120 



MVL + COA + 2 NADP <-> H3MCOA + 2 NADPH 

ATP + MVL-> ADP + PMVL 

CTP + MVL -> CDP + PMVL 

GTP + MVL -> GDP + PMVL 

UTP + MVL -> UDP + PMVL 

ATP + PMVL -> ADP + PPMVL 

ATP + PPMVL -> ADP + PI + IPPP + C02 

IPPP <-> DMPP 

GPP + IPPP -> FPP + PPI 

DMPP + IPPP -> GPP + PPI 

DMPP + IPPP -> GPP + PPI 

GPP + IPPP -> FPP + PPI 

2 FPP + NADPH -> NADP + SQL 
SQL + 02 + NADP -> S23E + NADPH 
S23E -> LNST 



2.8.2.4 
1.8.3.1 



2.3.1.85 

2.3.1.16 
2.3,1.16 
2,3.1.16 



1.3.3.6 

1.3.99.13 

1.3.99.7 

&2X3. 

6.2-1.3 

6.2.1.3 

2.3. 1 . 21 

2.3.1.21 

2.3.1.21 

5.3.3.8 

1.14.14.1 

1.14.14.1 

1.14.14.1 

1.14.14.1 

1.14.14.1 

1.14.14.1 

1.14.14.1 

1.14.14.1 

1.14.14.1 

1.14.14.1 

1,14.14.1 

1.14.14.1 

1.14.14.1 

1.14.14.1 

1.14.14.1 

1.14.14.1 

1.14.14.1 

1.14.14.1 

1.14.14.1 

1/14.14.1 

1.14.14.1 

1.14.14.1. 

1.14.14.1 

1. 14.14 .1 

1.14.14.1 



1.1.1.34 

2.7.1.36 



2.7.4.2 

4.1.1.33 

5.3.3.2- 

2.5.1. 10 

2.5.1.1 

2.5.1.1 

2.5.1.10 

2.5.1.29 

2.5.1.21 

1. 14 .99.7 

5.4.99.7 

1,6,99,2. 

1.6.99,2 

1.3.99.- 
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1055 CEU BSSL, BAL 

32BS LI PA, UL 

S0AT1, ACAT, STAT, SOAT, ACAT1, 

ACACT 
1531CYP7A1.CYP7 
5Z15SRD5A1 
6716 SRD5A2 

6718 AKR1D1, SRD5B1 , 3o5bred 
570 BAAT, BAT 
3.7 C21 -Steroid hormone metabolism PATH:hsa00140 
15B3'CYP11A, P450SCC 

3233 HSD3B1 , HSD3B, HSDB3 IMZYMST -> IIMZYMST + C02 

(MZYMST -> IIZYMST + C02 



3284 HSD3B2 



IMZYMST -> (IMZYMST + C02 
IMZYMST -> IIZYMST + C02 



1589 



CYP21A2, CYP21, P450C21B, 



1 CA21H, CYP21B. P450c21B 
15S6CYP17,P450C17 
1524 CYP11B1, P450C11, CYP11B 
l5fi5CYP11B2,CYP11B 

322Q HSD11B1, HSD11, HSD11L, HSD11B 

3221HSD11B2, HSD11K 
3.8 Androgen and estrogen metabolism PATH:hsa00150 

HSD17B1, EDH17B2, EDHB17, 

HSD17 
3223 HSD17B3, EDH17B3 
3294HSD17B2, EDH17B2 
3295HSD17B4 

oooc HSD17BP1, EDH17B1, EDHB17, 
^HSD^ 
51478 HSD17B7, PRAP 

412 STS, ARSC, ARSC1 , SSDD 
414ARSD 

415 ARSE, CDPX1, CDPXR, CDPX 
11185 INMT 
24140 JM23 

29104 N6AMT1 , PRED28 

2326QFJH1 
32Zfi HRMT1L2, HCP1, PRMT1 

51323 LOC51628 

54743 HASJ4442 

2Z222 HSA9761 
4. Nucleotide Metabolism 
4.1 Purine metabolism PATH:hsa00230 



3.1.1.3 

3,1.1,13 

3,1,1,13 

2.3.1.26 

1.14,13,17 

1.3.99.5 

1.3.99.5 



2^165 

1-14.15.6 
5.3.3.1 

1,1,1.145 

5.3.3.1 

1,1.1.145 

1.14.99.1Q 

1.14.99.9 

1.14.15 .4 

1.14.15.4 

1.1.1.146' 
1.1,1.146 

1.1.1.62 

1.1.1.62 
1.1.1.62 
,1.1.1.62 

1.1.1-62 

.1,1.1.62 
3.1.6.2 
3.1.6.1 
3,1.6.1 
2.1 .1 
2.1. 1 .- 
2.1.1 
2.1.1.- 
2.1.1.- 
2.1. 1.- 
2.1.1- 
2.1. 1.- 



11164 NUDT5, HYSAH1, YSA1H 




3.6.1.13 


54Z1 PPAT, GPAT 


PRPP + GLN -> PPI + GLU + PRAM 


2.4,2.14 


2618 GART, PGFT, PRGS 


PRAM + ATP + GLY <-> ADP + PI + GAR 


6.3.4.13 




FGAM + ATP -> ADP + PI + AIR 


6.3.3.1 




GAR + FTHF -> THF + FGAR 


Z1Z2 


5198 PFAS, FGARAT, KIAA0361, PURL 


FGAR + ATP + GLN -> GLU + ADP + PI + FGAM 




10606 ADE2H1 


CAIR + ATP + ASP <-> ADP + PI + SAICAR 


6.3.2.6 




CAIR <-> AIR + C02 


4.1.1.21 


£Q59 PAICS, AIRC, PAIS 


CAIR + ATP + ASP <-> ADP + PI + SAICAR 


fii3ZG_ 


153 ADSL 


ASUC<->FUM+AMP 


4.3.2.2 


4Z1ATIC, PURH 


AICAR + FTHF <-> THF + PRFICA 


2X2J3 




PRFICA <-> IMP 


3.5.4.10 


3251 HPRT1, HPRT, HGPRT 


HYXAN + PRPP -> PPI + IMP 






GN + PRPP -> PPI + GMP 




3614 IMPDH1 


IMP + NAD -> NADH + XMP 


1.1.1.205 


3615 IMPDH2 


IMP + NAD -> NADH + XMP 


11X205 


3333 GMPS 




6.3.5.2 


14923 






29SZGUK1 


GMP + ATP <-> GDP + ADP 


2.7.4.8 
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2988 GUK2 



10621 RPC39 

1S622RPC32 

JUI623RPC62 

11128 RPC155 

25335 DKFZP58BM01 22 

3QB34ZNRD1 

51Q82LOC51082 

5172B LOC51728 

5430 POLR2A, RPOL2, POLR2, POLRA 

5431 POLR2B, POL2RB 
5432POLR2C 

5433 POLR2D, HSRBP4, HSRPB4 

5434 POLR2E, RPB5, XAP4 

5435 POLR2F, RPB6, HRBP14.4 

5436 POLR2G, RPB7 

5437 POLR2H, RPB8, RPB17 
5435 POLR2I 
5439POLR2J 

5440 POLR2K, RPB7.0 

5441 POLR2L, RPB7.6, RPB10 

5442 POLRMT, APOLMT 
544Z9FU1Q816, Rp 0 1-2 
557 D3 FU1038B 

661 BN51T 
3533 RPA40, RPA39 
10721 POLQ 

H232 POLG2, MTPOLB, HP55, POLB 
2364EPOLA2 

5422 POLA 

5423 POLB 
5424POLD1, POLD 
5425POLD2 
5425 POLE 
542ZPOLE2 
5423 POLG 

59SQ REV3L, POLZ, REV3 
J49SXDH 

9615 GDA, KIAA1258, CYPIN, NEDASIN 
2Z66GMPR 
51292 LOC51292 
Z3ZZUOX 
524QRRM1 



5241 RRM2 



DGMP + ATP <-> DGDP + ADP 
GMP + DATP <-> GDP + DADP 
GMP + ATP<->GDP + ADP 
DGMP + ATP <-> DGDP + ADP 
GMP + DATP <-> GDP + DADP 



4860 NP, PNP 



122Q ECGF1, hPD-ECGF 

353 APRT 
132 ADK 
1633 DCK 



ADP + RTHIO -> DADP + OTHIO 
GDP + RTHIO -> DGDP + OTHIO 
CDP + RTHIO ■> DCDP + OTHIO 
UDP + RTHIO -> DUDP + OTHIO 
ADP + RTHIO -> DADP + OTHIO 
GDP + RTHIO -> DGDP + OTHIO 
CDP + RTHIO -> DCDP + OTHIO 
UDP + RTHIO -> DUDP + OTHIO 
AND + PI <->AD + R1P 
GSN + PI<-> GN + R1P 
DA + PI<->AD + R1P 
DG + PI <-> GN + R1P 
DIN + PK->HYXAN + R1P 
INS + PI<->HYXAN + R1P 
XTSINE + PI <->XAN + R1P 
DU + P1<->URA + DR1P 
DT + p| <-> THY + DR1P 
AD + PRPP->PPI + AMP 
ADN + ATP -> AMP + ADP 



2.7.4.8 



2J.7.6 
2, 7.7.6 
2.7.7.6 
2.7.7.6 
2.7.7.6 

2.7.7.5 

2.7.7.fi 

2JX6 
2,7.7.6 
2.7.7.6 

2.7.7.6 
2.7.7.6 
2.7.7.6 
2.7,7.6 

2.7.7.6 
2.7.7.6 
2.7.7.6 
2.7.7.6 
2.7.7.6 
2.7.7.6 
2.7.76 

2.7.7.6 

2.7.7.6 

2,7.7.6 

2.7.7.6 

2X12 

2.7.7.7 

2.7.7.7 

2.7.7.7 

2.7.7.7 

2.7.7.7 

2.7.7.7 

2.7.7.7 

2.7.7.7 

2.7.7.7 

2-7.7.7 

1.1.3.22 

1.1.1.204 

3.5.4.3 

1.6.6.8 

1.6.6.8 

17.3.3 

.1.17.4.1 



U7.4J 



2-4.2.4 

2.4.2.7 
Z1A2SL 
2.7 A J A 



WO 03/082214 

1716 DGUOK 

203 AK1 

204 AK2 

205 AK3 
26289 AK5 

4830 NME1, NM23, NM23-H1 

4fi31NME2,NM23-H2 
4832 NME3, DR-nm23, DR-NM23 
4S22NME4 

229ZS NT5B, PNT5, NT5B-PENDING 



4877 NT3 



79 



ATP + AMP <-> 2 ADP 

GTP + AMP <-> ADP + GDP 

ITP+AMP<-> ADP + IDP 

ATP + AMP <-> 2 ADP 

GTP + AMP <-> ADP + GDP 

ITP + AMP <-> ADP + IDP 

ATP + AMP <-> 2 ADP 

GTP + AMP <-> ADP + GDP 

ITP + AMP <-> ADP + IDP 

ATP + AMP <->2ADP 

GTP + AMP <-> ADP + GDP 

ITP + AMP <-> ADP + IDP 

UDP + ATP <-> UTP + ADP 

CDP + ATP <-> CTP + ADP 

GDP + ATP <-> GTP + ADP 

IDP + ATP<->ITP + IDP 

DGDP + ATP <-> DGTP + ADP 

DUDP + ATP <-> DUTP + ADP 

DCDP + ATP <-> DCTP + ADP 

DTDP + ATP <-> DTTP + ADP 

DADP + ATP <-> DATP + ADP 

UDP + ATP <-> UTP + ADP 

CDP + ATP <-> CTP + ADP 

GDP + ATP <-> GTP + ADP 

IDP + ATP <-> ITP + IDP 

DGDP + ATP <-> DGTP + ADP 

DUDP + ATP <-> DUTP + ADP 

DCDP + ATP <-> DCTP + ADP 

DTDP + ATP <-> DTTP + ADP 

DADP + ATP <-> DATP + ADP 

UDP + ATP <-> UTP + ADP 

CDP + ATP <-> CTP + ADP 

GDP + ATP <-> GTP + ADP 

IDP + ATP <-> ITP + IDP 

DGDP + ATP <-> DGTP + ADP 

DUDP + ATP <-> DUTP + ADP 

DCDP + ATP <-> DCTP + ADP 

DTDP + ATP <-> DTTP + ADP 

DADP + ATP <-> DATP + ADP 

UDPm + ATPm <-> UTPm + ADPm 

CDPm + ATPm <-> CTPm + ADPm 

GDPm + ATPm <-> GTPm + ADPm 

IDPm + ATPm <-> ITPm + IDPm 

DGDPm + ATPm <-> DGTPm + ADPm 

DUDPm + ATPm <-> DUTPm + ADPm 

DCDPm + ATPm <-> DCTPm + ADPm 

DTDPm + ATPm <-> DTTPm + ADPm 

DADPm + ATPm <-> DATPm + ADPm 

AMP + H20 -> PI + ADN 

GMP->PI + GSN 

CMP->CYTD + PI 

UMP -> PI + URI 

IMP->PI + INS 

DUMP->DU + PI 

DTMP -> DT + PI 

DAMP->DA + PI 

DGMP -> DG + PI 

DCMP->DC + PI 

XMP->PI+XTSINE 

AMP ■> PI + ADN 

GMP->PI + GSN 

CMP->CYTD + PI 

UMP -> PI + URI 

IMP -> PI + INS 

DUMP •> DU + PI 

DTMP -> DT + PI 



PCT/US03/09751 

2.7.1.113 

2.7.4.3 

2.7.4.3 
2.7.4.3 
2.7.4.3 
2.7.4.6. 

27.4.6 

2.7.4.6 

2.7.4.6 

3 1.3.5 



3 1.3.5 
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490? NTS, CD73 



7370 UMPH2 



1Q846 PDE10A 



27115 PDE7B 



513SPDE1A 



513ZPDE1C, HCAM3 



5138 PDE2A 



5139 PDE3A, CGI-PDE 



5140 PDE3B 
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DAMP->DA + PI 
DGMP->DG + PI 
DCMP -> DC + PI 
XMP->PI+XTSINE 

AMP -> PI + ADN 3 ,1 ,3 .5 

GMP->Pl + GSN 

CMP ->CYTD + PI 

UMP->PI + URI 

IMP->Pl + INS 

DUMP->DU + PI 

DTMP->DT + PI 

DAMP->DA + PI 

DGMP->DG + PI 

DCMP -> DC + PI 

XMP->PI + XTSINE 

AMP -> PI + ADN 3 . 1 .3,5 

GMP -> PI + GSN 

CMP -> CYTD + PI 

UMP-> PI + URI 

IMP -> PI + INS 

DUMP->DU + PI 

DTMP -> DT + PI 

DAMP->DA+PI 

DGMP->DG + PI 

DCMP -> DC + PI 

XMP->PI+XTSINE 

cAMP -> AMP 3.1,4 . 17 - 

cAMP -> AMP 

cdAMP -> dAMP 

cIMP -> IMP 

cGMP->GMP 

cCMP -> CMP 

cAMP -> AMP 3 , 1.4.17 

cAMP -> AMP 

cdAMP -> dAMP 

cIMP -> IMP 

cGMP->GMP 

cCMP -> CMP 

cAMP -> AMP 3 . 1.4 . 17 

cAMP -> AMP 

cdAMP -> dAMP 

cIMP -> IMP 

cGMP->GMP 

cCMP -> CMP 

cAMP -> AMP 3.1 . 4 . 17 

cAMP->AMP 

cdAMP -> dAMP 

cIMP-> IMP 

cGMP -> GMP 

cCMP -> CMP 

cAMP -> AMP 3,1,4.17 

cAMP->AMP 

cdAMP -> dAMP 

clMP->IMP 

cGMP -> GMP 

cCMP -> CMP 

cAMP -> AMP 3,1 . 4.17 

cAMP -> AMP 

cdAMP -> dAMP 

dMP-> IMP 

cGMP -> GMP 

cCMP -> CMP 

cAMP->AMP 3 . 1,4,17 . 

cAMP ->AMP 

cdAMP -> dAMP 

cIMP -> IMP 

cGMP->GMP 
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5141 PDE4A, DPDE2 

5142 PDE4B, DPDE4, PDEIVB 

5343 PDE4C, DPDE1 

5144 PDE4D, DPDE3 

5145 PDE6A, PDEA, CGPR^A 

5146 PDE6C, PDEA2 

5147 PDE6D 

5148 PDE6G, PDEG 

5149 PDE6H 
5152 PDE9A 



S152PDES1B 



PDE6B, CSNB3, PDEB 
8654 PDE5A 
ACQ ADA 

2ZQ AMPD1 , MADA 
271 AMPD2 
222AMPD3 
253ENTPD1,CD39 
3704 ITPA 

107 ADCY1 

108 ADCY2, HBAC2 

109 ADCY3, AC3, KIAA0511 
Hi}ADCY4 

111 ADCY5 

112 ADCY6 

113 ADCY7, KIAA0037 

114 ADCY8, ADCY3, HBAC1 
115ADCY9 

2977 GUCY1A2, GUC1A2, GC-SA2 
oqqo GUCY1A3, GUC1A3, GUCSA3, GC- 
2982 SA3 

„„ GUCY1B3, GUC1B3, GUCSB3, GO 
2233 SB3 

22B4 GUCY2C, GUC2C, STAR 

oooe GUCY2F, GUC2F, GC-F, GUC2DL, 

™RETGC-2 



cCMP -> CMP 
cAMP -> AMP 
CAMP -> AMP 
cAMP -> AMP 
cAMP -> AMP 
cGMP->GMP 
cGMP -> GMP 
cGMP -> GMP 
cGMP -> GMP 
cGMP -> GMP 
cAMP -> AMP 
cAMP-> AMP 
cdAMP -> dAMP 
cIMP -> IMP 
cGMP -> GMP 
cCMP -> CMP 
cAMP -> AMP 
cAMP -> AMP 
cdAMP -> dAMP 
cIMP -> IMP 
cGMP -> GMP 
cCMP -> CMP 
cGMP -> GMP 
cGMP -> GMP 
ADN->1NS + NH3 
DA-> D1N + NH3 
AMP-> IMP + NH3 
AMP -> IMP + NH3 
AMP -> IMP + NH3 



3000 



GUCY2D, CORD6, GUC2D, LCA1, 



2 GUC1A4,LCA,retGC 
4881 NPR1 t ANPRA, GUC2A, NPRA 
, NPR2, ANPRB, GUC2B, NPRB, 
' NPRBi 
159ADSS 
318 NUDT2, APAH1 

5167 



ENPP1, M6S1, NPPS, PCA1, PC-1, 
L PDNP1 

53fia ENPP2, ATX, PD-IALPHA, PDNP2 
5169 ENPP3, PD-IBETA, PDNP3 

22Z2FHIT 

4.2 Pyrimidine metabolism PATH:hsa00240 
790 CAD 



1723 DHODH 
7372 UMPS, OPRT 



ATP -> 
ATP-> 
ATP-> 
ATP-> 
ATP -> 
ATP-> 
ATP-> 
ATP-> 
ATP-> 



cAMP + PPI 
cAMP + PPI 
cAMP + PPI 
cAMP + PPI 
cAMP + PPI 
cAMP + PPI 
cAMP + PPI 
cAMP + PPI 
cAMP + PPI 



3.1.4.17 
3-1.4.17 
3,1.4.17 
3.1.4.17 
3.1.4.17 
3.1.4,17 
3.1.4.17 
3.1.4.17 
3.1.4.17 
3.1.4.17 



3.1.4.17 



3.1.4.17 
3.1.4.17 

3.5.4.4 



IMP + GTP + ASP -> GDP + PI + ASUC 



GLN + 2 ATP + C02 -> GLU + CAP + 2 ADP + PI 

CAP + ASP->CAASP + PI 

CAASP <-> DOROA 

DOROA + 02 <-> H202 + OROA 

OMP->C02 + UMP 



3.5.4.6 

3.5.4.6 

3.5.4.6 

3.6.1.5 

3.6.1.19 

4.6.1,1 

4.6.1.1 

4.6.1.1 

4.6.1.1 ' 

4.6.1.1 

4.6.1.1 

4,6,1. 1 

4,6.1.1 

4.6.1.1 

4.6.1.2 

4.6.1.2 

4.6.1.2 
4.6.1.2 
4.6.1.2 

4,6.1.2 
4-61.2 

.4.6.1.2 

6.3.4.4 
3.6.1.17 

3-6.1.9 

3.6.1.9 
3.6.1.9 
• 3.1.4.1 
3 6.1.29 

6.3.5.5 

3*5^3. 
1.3,3,1 
4.1.1.23 
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51Z2ZLOC51727 



50808 AKL3L 
15Q3CTPS 

ZSHUMPK.TSASOS 



7378 UP 

1806 DPYD, DPD 

1B07 DPYS, DHPase, DHPASE, DHP 
51733 LOC51733 
7296 TXNRD1,TXNR 
4654 DUT 

7298 TYMS, TMS, TS 
978 CDA, CDD 

0535 DCTD 

7083 TK1 

7084 TK2 



1841 DTYMK, TYMK, CDC8 
4.3 Nucleotide sugars metabolism PATH:hsa00520 
23463 TDPGD 

1486 CTBS, CTB 
5. Amino Acid Metabolism 
5.1 Glutamate metabolism PATH:hsa00251 

£659 ALDH4, P5CDH 

2QS6 EPRS t OARS, QPRS 



OROA + PRPP <-> PPI + OMP 
ATP + UMP<->ADP + UDP 
CMP + ATP<->ADP + CDP 
DCMP + ATP <-> ADP + DCDP 

UTP + GLN + ATP-> GLU +CTP + ADP + Pl 
ATP + UTP + NH3 -> ADP + PI + OTP 
URl + ATP->ADP + UMP 
URI + GTP -> UMP + GDP 
CYTD + GTP -> GDP + CMP 
URI + Pl<-> URA + R1P 



OTHIO + NADPH -> NADP + RTHIO 
DUTP -> PPI + DUMP 
DUMP + METTHF -> DHF + DTMP 
CYTD->URI + NH3 
DC -> NH3 + DU 
DCMP <-> DUMP + NH3 
DU + ATP -> DUMP + ADP 
DT + ATP -> ADP + DTMP 
DUm + ATPm -> DUMPm + ADPm 
DTm + ATPm -> ADPm + DTMPm 
DTMP + ATP <-> ADP + DTDP 



P5C + NAD + H20 -> NADH + GLU 
GLU + ATP -> GTRNA + AMP + PPI 



F6P + GLN -> GLU + GA6P 
F6P + GLN -> GLU + GA6P 



CYS + GLU + ATP-> GC + PI + ADP 
C YS + G LU + ATP -> GC + PI + ADP 
GLY + GC + ATP -> RGT + PI + ADP 
NADPH + OGT -> NADP + RGT 



2623 GFPT1, GFA, GFAT. GFPT 
9945 GFPT2, GFAT2 
5859 QARS 

2729 GLCLC, GCS, GLCL 
223Q GLCLR 
2932 GSS, GSHS 
2936 GSR 
5188 PET112L, PET112 

5.2 Alanine and aspartate metabolism PATH:hsa00252 

4677 MARS, ASNRS ATP + ASP + TRNA -> AMP + PPI + ASPTRNA 

425 ASL ■ ARGSUCC->FUM+ARG 
189 AGXT, SPAT SERm + PYRm <-> ALAm + 3HPm 

^ ALA + GLX <-> PYR + GLY 

JUaAARS 

1615 DARS 

445 ass, CTLN1 , ASS1 CITR + ASP + ATP <-> AMP + PPI + ARGSUCC 

443 ASPA t ASP,ACY2 
1364 CRAT, CAT1 

ACCOA + CAR -> COA + ACAR 

6526 DDO 

5.3 Glycine, serine and threonine metabolism PATH:hsaQ0260 



9-4.2.10 
2.7.4.14 



2.7.4.10 
6.3.4.2 



5223 PSPH, PSP 
29968 PSA 

26222 PHGDH, SERA, PGDH, PGD, PGAD 
23464 GCAT, KBL 

211 ALAS 1. ALAS 

212 ALAS2, ANH1, ASB 

4128 MAOA 

4129 MAOB 

26 ABP1, AOC1.DAO 
314 AOC2, DA02, RAO 
6639 AOC3, VAP-1 , VAP1 , HPAO 
2231 GLDC 



3PSER + H20->PI + SER 
PHP + GLU <-> AKG + 3PSER 
OHB + GLU <-> PHT + AKG 
3PG + NAD <-> NADH + PHP 

SUCCOA + GLY -> ALAV + COA + C02 
SUCCOA + GLY -> ALAV + COA + C02 
AMA + H20 + FAD -> NH3 + FADH2 + MTHGXL 
AMA + H20 + FAD •> NH3 + FADH2 + MTHGXL 



GLY + UPO <-> SAP + C02 



2.7.1 .48 



2.4.2.3 

1 ,3.1. 2. 

3.5.2.2 

3.5.1.6 

1-6.4.5 

3 6.1.23 

2.1.1.45 

3.5.4.5 



3.5.4.12 
2-7.1.21 

2,7.1.21 

2.7.4.9 

4.2.1.46 

3.2.1. - 



1,5.1.12 
6.1.1.17 
6.1.1.15 

2.6.1.16 

2.6.1.16 

6,1.1.18 

6.3.2.2 

63.2.2 

6^Z3_ 

1.6.4.2 

6.3.5.- 

6X122. 

4.3.2.1 

2.6.1.51 

2.6.1. 44 

6.1.1.7 

6.1.1.12 

fi-3.4.5 

3.51.15 

2.3.1.7 

.1.4.3.1 

3.1.3,3 
Z6JL52. 

1.1.1.95 

2.3.1.29 

2 3.1.37 

23.1.37 

14.3.4 

1A3.4. 

1A3.6 

14,3,6 

1.4.3.6 

14.4.2 
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1510. DAO, DAMOX 
2512 GARS 
2622 GATM 
2523 GAMT 

PISD, PSSC, DKFZP566G2246, 
23761 DJ858B16 

£25 BHMT 
23258 DMGDH 
fiZ5 CBS 

52d1 sars, sers 
10233 sds, sdh 

522ZTARS 
5.4 Methionine metabolism PATH:hsa00271 



PS-> PE + C02 

SER + HCYS -> LLCT + H20 
SER->PYR + NH3 + H20 



4143 MAT1A, MATA1, SAMS1 , MAT, SAMS MET + ATP + H20 -> PPI + PI + SAM 



MET + ATP + H20 ■> PPi + PI + SAM 
SAM + DNA -> SAH + DNA5MC 
SAH + H20 -> HCYS + ADN 
SAH + H2Q-> HCYS + ADN 



4144 MAT2A, MATA2, SAMS2, MATH 
IZfifi DNMT1, MCMT, DNMT 
10768 AHCYL1, XPVKONA 

191 AHCY, SAHH 
4141 MARS, METRS, MTRNS 
4543 MTR 

5.5 Cysteine metabolism PATH:hsaO0272 

833 CARS 
1G36 CDOI 
2509 NDST2, HSST2, NST2 

5.6 Valine, leucine and isoleucine degradation PATH:hsa00280 

526 BCAT1 , BCT1 , ECA39, MECA39 AKG + ILE -> OMVAL + GLU 

AKG + VAL -> OIVAL + GLU 



HCYS + MTHF ->THF + MET 



CYS + 02<->CYSS 



1.4.3.3 

fi.a.1.14 

2.1.4.1 
2X12. 

4X155 

2.1.1 .5 

1^22x2. 

42X22. 

6.1.1.11 

4 .2.1.13 

6.1.1.3 

2.3-1.6 

2.5.1.6 

2.1.1.37 

3.3.1.1 

3.3.1.1 

6.1.1.10 

2.1.1.13 

6.1.1.16 
1.13.11.20 

2.8.2.- 

2.6.1.42 



587 BCAT2, BCT2 

5014 OVD1A 
592 BCKDHA, MSUD1 



594 BCKDHB, E1B 



3712 IVD 
215 AOX1.AO 

4164 MCCC1 

4165 MCCC2 



AKG + LEU -> OICAP + GLU 
OICAPm + GLUm <-> AKGm + LEUm 
OMVALm + GLUm <-> AKGm + JLEm 

OMVALm + COAm + NADm -> MBCOAm + NADHm + C02m 
OIVALm + COAm + NADm -> IBCOAm + NADHm + C02m 
OICAPm + COAm + NADm -> IVCOAm + NADHm + C02m 
OMVALm + COAm + NADm -> MBCOAm + NADHm + C02m 
OIVALm + COAm + NADm -> IBCOAm + NADHm + C02m 
OICAPm + COAm + NADH -> IVCOAm + NADHm + C02m 
IVCOAm + FADm -> MCRCOAm + FADH2m 

MCRCOAm + ATPm + C02m + H20m -> MGCOAm + ADPm + 
Pirn 

MCRCOAm + ATPm + C02m + H20m -> MGCOAm + ADPm + 
Pirn 



5.7 Valine, leucine and isoleucine biosynthesis PATH:hsa00290 
22325 K1AA0028, LARS2 

2225 LARS 
22Z6 IARS, ILRS 
24D6 VARS1, VARS 
Z4QZ VARS2, G7A 

5.8 Lysine biosynthesis PATH:hsaO030Q 

3735 KARS, KIAA0070 ATP + LYS + LTRNA -> AMP + PPI + LLTRNA 

5.9 Lysine degradation PATH:hsa00310 

8424 BBOX, BBH. GAMMA-BBH, G-BBH 

5351 PLOD. LLH 

5352 PLOD2 
2925 PLOD3, LH3 

10157 LKR/SDH, AASS LYS + NADPH + AKG -> NADP + H20 + SAC 

SAC + H20 + NAD -> GLU + NADH + AASA 

5.10 Arginine and proline metabolism PATH:hsa00330 

5QQ9 OTC ORNm + CAPm -> CITRm + Pirn + Hm 

223 ARG1 ARG -> ORN + UREA 

234 ARG2 ARG -> ORN + UREA 

4842 NOS1.NOS 

4242 NOS2A, NOS2 

4346 NOS3, ECNOS 
, 4242 OAT ORN + AKG <-> GLUGSAL + GLU 



2.6.1.42 

1.2.4.4 
1.2.4.4 



1.2.4.4 



1,3.99,10 

1.2.3.1 

6.4.1.4 
6.4.1.4 



6.4.1.4 
6.4.1.4 
6.1.1.5 
6.1.1.9 
6.1.1.9 

6.1.1.6 

1J 4 .1 1 .1 
1.14.11.4 
1-14.11.4 
1.14.11.4 
1.5.1.9 



2X2£ 

3.5.3.1 

3-5.3.1 

1,14.13.39 

1-14.13.39 

1 14.13.39 

2.6.1.13 
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5831 PYCR1, P5C, PYCR 



5023 P4HA1, P4HA 
5912 RARS 
1152 CKB, CKBB 
1156 CKBE 

1158 CKM, CKMM 

1159 CKMT1, CKMT, UMTCK 

1160 CKMT2, SMTCK 
5223 SRM, SPS1,SRML1 

252 AMD1, ADOMETDC 

262 AMDP1, AMD, AMD2 
1725 DHPS 
6611 SMS 
4253 ODC1 
£303 SAT, SSAT 

5.11 Histidine metabolism PATH:hsa00340 
10841 FTCD 

3062 HDC 
1644DDC.AADC ■ 
3126 HNMT 
218 ALDH3 

220 ALDH6 

221 ALDH7, ALDH4 

222 ALDH8 
3035 HARS 

5.12 Tyrosine metabolism PATH:hsa00350 
6898 TAT 

3242 HPD, PPD 
3081 HGD, AKU, HGO 
2954 GSTZ1 , MAAl 

2184 FAH 

2239 TYR, OCAIA 

2G54TH,TYH 

1621 DBH 

5409 PNMT, PENT 

1332 COMT 

2123TP0,TPX 

5.13 Phenylalanine metabolism PATH:hsa00360 

501 ATQ1 

5.14 Tryptophan metabolism PATH:hsa00380 
6999 TD02, TPH2, TRPO, TDO 

8564 KMO 
8942 KYNU 



P5C + NADPH -> PRO + NADP 
P5C + NADH -> PRO + NAD 
PHC + NADPH -> HPRO + NADP 
PHC + NADH -> HPRO + NAD 

ATP + ARG + ATRNA -> AMP + PPI + ALTRNA 
PCRE + ADP -> CRE + ATP 



PTRSC + SAM -> SPRMD + 5MTA 
SAM <-> DSAM + C02 
SAM <-> DSAM + C02 
SPRMD + Qm -> DAPRP + QH2m 
DSAM + SPRMD -> 5MTA + SPRM 
ORN-> PTRSC + C02 



FIGLU + THF -> NFTHF + GLU 



ACAL + NAD •> NADH + AC 
ACAL + NAD -> NAD H + AC 
ACAL + NAD -> NADH + AC 
ACAL + NAD -> NADH + AC 
ATP + HIS + HTRNA -> AMP + PPI + HHTRNA 

AKG + TYR -> HPHPYR + GLU 
HPHPYR + 02 -> HGTS + C02 
HGTS + 02 -> MACA 
MACA -> FACA 

FACA + H20 -> FUM + ACA 



TRP + G2-> FKYN 

KYN + NADPH + 02 -> HKYN + NADP + H20 
KYN -> ALA + AN 
HKYN + H20 -> HAN + ALA 
HAN+02->CMUSA 



ATPm + TRPm + TRNAm -> AMPm + PPlm + TRPTRNAm 
ATP + TRP + TRNA -> AMP + PPI + TRPTRNA 



23498 HAAO, HAO, 3-HAO 
7166 TPH, TPRH 
438 ASMT, HIOMT, ASMTY 
15 AANAT, SNAT 
3620 INDO, IDO 
10352 WARS2 
7453 WARS, IFP53, IF153, GAMMA-2 
4234 NEDD4, KIAA0093 

5.15 Phenylalanine, tyrosine and tryptophan biosynthesis PATH:hsa00400 

5053 PAH, PKU1 PHE + THBP + 02 -> TYR + DHBP + H20 

10662 FARS1 

2133 FARSL, CML33 
10Q56 PheHB 

3565 YARS, TYRRS, YTS, YRS 

5.16 Urea cycle and metabolism of amino groups PATH:hsa00220 
5632 PYCS 

GLUP + NADH -> NAD + PI + GLUGSAL 
GLUP + NADPH -> NADP + PI + GLUGSAL 



1.5.1.2 



1.14.11.2 

6.1.1.19 

2.7.3.2 

2.7.3.2 

2.7.3.2 

2.7.3.2 
2-5.1.16 

4,1.1.50 
4.1.1,50. 

1.5.99.6 
2^L22_ 
4.1.1.17 
2.3.1.57 

4.3.1.4 

4.1.1.22 

4.1.1.28 

2.1.1.8 

1.2.1.5 

1.2.1,5 

1.2.1.5 

.1.2.1.5 

6.1.1.21 



2.6.1.5 

113-11 -27 

1-13.11.5 

5.2.1.2 

2.5.1.18 

3.7.1.2 

1.14.18.1 

1.14.16.2 

1.14.17.1 

2.1.1.28 

1 .11.1.8. 

I.2.1.- 

1.13.11.11 
1,14.13.9 
3.7.1.3 

1-13.11.6 

1.14.16.4 

2.1.1.4 

2.3.1,67 

1.13,11.42 

6X12 

6J.1.2 

£^2w 



1-14.16.1 
6.1.1.20 
6.1.1.20 
61.1.20 
6.1.1.1 



7, 7.2.11 
1-2.1.41 



WO 03/082214 



PCT7US03/09751 



85 



95ACY1 3.5.1.14 
6. Metabolism of Other Amino Adds 

6.1 beta-Alanine metabolism PATH :hsa00410 

6.2 Taurine and hypotaurine metabolism PATH:hsa00430 

25Z8 JJJ 1 ' GTG ' D22S672 « D22S732 « RGT + AU-> CGLY + ALAGLY 2^L 

2679 GGT2, GGT RGT +ALA-> CGLY + ALAGLY 

2680 GGT3 RGT + ALA -> CGLY + ALAGLY 2322. 
268Z GGTLA1 , GGT-REL, DKFZP566O01 1 RGT + ALA -> CGLY + ALAGLY 

6.3 Aminopbosphonate metabolism PATH:hsa00440 

5130 PCYT1A, CTPCT, CT, PCYT1" PCHO + CTP -> CDPCHO + PPI 2,7.7.15 

9791 PTDSS1 , KIAA0024, PSSA CDPDG + SER <-> CMP + PS 2,7.8,- 

6.4 Selenoamino acid metabolism PATH:hsa00450 

22928 SPS2 2,7.9.3. 

22929 SPS, SELD 2.7.9.3. 

6.5 Cyanoamlno acid metabolism PATH:hsa00460 

6.6 D-Glutamine and D-glutamate metabolism PATH:hsa00471 

6.7 D-Arginine and D-omithlne metabolism PATH:hsa00472 
6.9 Glutathione metabolism PATH:hsa00480 

5182 PEPB 3,4.11,4 

2655 GCTG Z2L2A 

28Z5 GPX1.GSHPX1 2 RGT + H202 <-> OGT 1.11.1.9 

2877 GPX2, GSHPX-GI 2 RGT + H202 <-> OGT 1.11.1.9 

2878 GPX3 2 RGT + H202 <-> OGT 1.11.1.9 
^Z9 GPX4 2 RGT + H202 <-> OGT 1.11.1.9 
288Q GPX5 2RGT+H202<->OGT 1.11.1.9 
2881 GPX6 2 RGT + H202 <-> OGT 1.11.1.9 
2330 GSTA1 2,5.1.18 
2333 GSTA2, GST2 2.5.1,18 

2940 GSTA3 2.5.1.18 

2941 GSTA4 2.5.1.18 
2944 GSTM1, GST1, MU 2.5.1.18 

2946 GSTM2, GST4 2.5,1.18 

2947 GSTM3, GST5 2.5,1.18 

2948 GSTM4 2.5.1.18 

2949 GSTM5 2.5.1.18 

2950 GSTP1, FAEES3, DFN7, GST3, PI 2.5,1.18 

2252 GSTT1 2.5.1.18 

2253 GSTT2 2r5.1.18 

4257 MGST1 , GST12, MGST, MGST-1 2.5.1.18 

4258 MGST2, GST2 r MGST-JI 2,5.1.18 
4253 MGST3, GST-Ill 2.5.1.18 

7. Metabolism of Complex Carbohydrates 

7.1 Starch and sucrose metabolism PATH:hsaO05O0 

6476 SI 3.2,1. 1Q 

3.2.1.48 

111B1 TREH, TRE, TREA TRE -> 2 GLC &2X28_ 
2990 GUSB 

2532 GBE1 GLYCOGEN + PI -> G1 P 2.4,118 

5834 PYGB GLYCOGEN + PI -> G1 P 2.4.1.1 

5822 PYGL GLYCOGEN + PI -> G1 P 2.4.1.1 

5S3Z PYGM GLYCOGEN + PI -> G1 P ZAAA 

233Z GYS1, GYS UDPG -> UDP + GLYCOGEN 2.4.1.11 

2398 GYS2 UDPG -> UDP + GLYCOGEN 2.4.1,11 

276 AMY1A, AMY1 3.2.1,1 

277 AMY1B, AMY1 3.2.1,1 
2Z8 AMY1C, AMY1 2u2AA 
2Z9 AMY2A, AMY2 .3,2.1.1 
2BQ AMY2B, AMY2 3.2.1.1 
178 AGL, GDE 2,4.1.25 

3.2.1.33 

10000 AKT3, PKBG, RAC-GAMMA, PRKBG 2.7.1 c 

1017 CDK2 2.7.1.- 

1018 CDK3 2.7.1- 

1019 CDK4, PSK-J3 2J.1- 

1020 CDK5, PSSALRE 2.7. 1.- 
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1021 CDK6, PLSTIRE 

1022 CDK7, CAK1, STK1, CDKN7 
1024 CDK8, K35 

1G25 CDK9, PITALRE, CDC2L4 
1029S PAK4 
10746 MAP3K2, MEKK2 

H11CHEK1.CHK1 
11200 RAD 53, CHK2, CDS1, HUCDS1 

H95CLK1,CLK 

1326 MAP3K8, COT, EST, ESTF, TPL-2 
w> MAPK 14 - CSBP2. CSPB1, PRKM14. 
J£a ^ PRKM15, CSBP1, P38, MXI2 
1452 CSNK1A1 
1452CSNK1D, HCKID 

1454 CSNK1E, HCKJE 

1455 CSNK1G2 

1456 CSNK1G3 
1612 DAPK1, DAPK 

1760 DMPK, DM, DMK, DM1 

1£59 DYRK1A, DYRK1, DYRK, MNB, MMBH 

2QS AKT2, RAC-BETA, PRKBB, P KB BETA 

259 AMHR2, AMHR 
2233Q RPS6KA6, RSK4 

2608 GPRK2L, GPRK4 

28SS GPRK5, GRK5 

2870 GPRK6, GRK6 
29904 HSU93850 
30811 HUNK 

3B11 ILK, P59 

3654 IRAKI, IRAK 

369 ARAF1, PKS2, RAFA1 

370 ARAF2P, PKS1 r ARAF2 
2924 UMK1, LIMK 

3985 LIMK2 

4112 MAK 

4140 MARK3, KP78 

4215 MAP3K3, MAPKKK3, MEKK3 

MAP3K4, MAPKKK4, MTK1, MEKK4, 

KIAA0213 

4217 MAP3K5, ASK1, MAPKKK5, MEKK5 
4292 MAP3K9, PRKE1, MLK1 
4294 MAP3K10, MLK2, MST 
4242 MOS 

4751 NEK2, NLK1 

4752 NEK3 

5053 PAK1, PAKalpfia 
5052 PAK2, PAK65, PAKgamma 
5063 PAK3, MRX30, PAK3beta 
5122 PCTK1, PCTGAIRE 

5128 PCTK2 

5129 PCTK3, PCTA1RE 
5292PIM1,PIM 
5342 PLK, PLK1 
5562 PRKAA1 

5562 PRKAA2, AMPK, PRKAA 

5522 PRKCA, PKCA 

5529 PRKCB1, PKCB, PRKCB, PRKCB2 

55SQ PRKCD 

5581 PRKCE 

5522 PRKCG, PKCC, PKCG 

55S3 PRKCH, PKC-L, PRKCL 

5524 PRKCI, DXS1179E, PKCI 

55S5 PRKCL1, PAK1, PRK1, DBK, PKN 

5526 PRKCL2, PRK2 

55SS PRKCQ 



86 

2.7. 1.- 
2.7.1.n 
2.7.1.- 
2.7.1. - 
2.7. 1.- 
2.7.1,-. 
2.7.1.- 
2.7.1- 
2.7.1.- 
2,7.1,- 

2.7.1. - 

2.7.1 . - 

2,7 .1,- 

2.7.1. - 

2.7.1 .- 

2.7.1- 

2.7.1 

2.7.1 

2.7.1- 

2.7.1 .- 

2.7.1. - 
2.7.1 
2,7.1, - 
2.7.1 
2.7 A.- 
2 .7.1. - 
2.7. t= 
2.7.1. - 
2.7.1. - 
2.7.1.- 
2.7.1.- 
2.7.1- 
2.7.1 
2.7.1. - 
2.7.1. - 
2.7.1. - 

2.7.1.- 

2.7.1- 

2.7.1 

2.7.1 .- 

2.7.1 

2,7.1.- 

2.7.1. - 

2,7.1 r 

2.7.1.- 

2.7.1 .- 

2-7.1.- 

2.7.1 .- 

2.7.1.- 

2-7-1. - 

2.7.1.- 

2-7.1.- 

2.7.1. - 
2.7. 1. - 
2-71- 
2.7.1 . - 
2-7.1.- 
2-7.1.- 
. 2.7.1.- 
2X1= 
2.7.1. - 
2.7.1. - 
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559Q PRKCZ 2.7.1. - 
MAPK1, PRKM1, P41MAPK, 

5524 P42MAPK, ERK2, ERK, MAPK2, 2.7.1. - 
PRKM2 

ccqc MAPK3, ERK1, PRKM3, P44ERK1, . _ , 

5325 P44MAPK 2JJX ^ 

559Z MAPK6, PRKM6, P97MAPK, ERK3 2.7.1. - 

5593 MAPK7, BMK1, ERK5, PRKM7 2.7.1 



MAP KB, JNK, JNK1, SAPK1, PRKM8, 
JNK1A2 

MAPK9, JNK2, PRKM9, P54ASAPK, 



PRKMK6 



2.7.1. - 



5m JUNKINASE 2XLz 
MAPK10, JNK3, PRKM10. P493F12, 

P54BSAPK ILL** 
MAPK13, SAPK4, PRKM13, 

P38DELTA * ,7 ,1 '~ 
MAP2K1 , MAPKK1 , MEK1 , MKK1 , 

PRKMK1 2 , 7 , 1 c 

5505 MAP2K2, MEK2, PRKMK2 2.7.1 

5505 MAP2K3, MEK3, MKK3, PRKMK3 2.7.1 .> 

5607 MAP2K5, MEK5, PRKMK5 2.7.1. - 
. MAP2K6, MEK6, MKK6, SAPKK3, 



2.7.1. - 



ccrsn MAP2K7, MAPKK7, MKK7, PRKMK7, 

5SQ9 JNKK2 2X1=. 

561Q PRKR, EIF2AK1, PKR 2.7.1. - 

5513 PRKX, PKX1 2.7.1 

5m RAF1 2 ,7 ,1 - 

e j 0 BCR, CMU PHL, BCR1 , D22S1 1 , _ _ < 

^ D22S662 . 22A ^' 

M0 - RPS6KA1, HU-1. RSK, RSK1. 

^ MAPKAPK1A 2JJ ^ 
R1QR RPS6KA2, HU-2, MAPKAPK1C, RSK, 

RPS6KA3, RSK2, HU-2, HU-3, RSK, 

MAPKAPK1B, ISPK-1 *'' f1r 

6198 RPS6KB1, STK14A 2.7.1. - 

6199 RPS6KB2, P70-BETA, P70S6KB 2.7.1. > 
ft<anri MAPK1 2, ERK6, PRKM12, SAPK3, _ _ „ 
52252 P38GAMMA, SAPK-3 2lZ ^ 

MAP2K4, JNKK1, MEK4, PRKMK4. 

524112 SERK1 , MKK4 * ,M 

6446 SGK 2.7.1 

555 BMPR1B, ALK-6, ALK6 * . • ' ■ 2.7.1- 

552 BMPR2, BMPR-lf, BMPR3, BRK-3 2.7.1. - 

523 BRAF 2.7.1. - 

6792 STK9 2.7.1.- 

6794 STK1 1 , LKB1 , PJS 2.7.1 

5555 MAP3K7, TAK1 2.7.1 

£92 BUB1 2.7.1. - 

Z01 BUB1B, BUBR1, MAD3L 2.7.1 

7016 TESK1 2.7. 1.- 

J222 TTK, MPS1L1 2.7.1 

286Z MAPKAPK3, 3PK, MAPKAP3 2.7.1 

8408 ULK1 2.7.1. - 

5558 CDK10, PISSLRE 2.7,1. - 

5621 CDC2L5, CDC2L, CHED 2.7.1 

5Z3IRIPK1.RIP 2.7.1. . 

8814 CDKL1, KKIALRE 2.7.1 

8899 PRP4, PR4H 2.7.1. - 

2054 MAP3K6, MAPKKK6 2.7.1. - 

9149 DYRK1B 2.7.1- 

32 ACVR2, ACTRll 2.7.1,, 

2201 DCAMKL1, KIAA0369 2.7.1 

23 ACVR2B 2.7.1 

253 CDC2 2.7. 1.- 

984 CDC2L1 2.7.1. - 
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5205 FIC1, BFUC, PFIC1, PFIC, ATP8B1 



DHPP->DHP + PI 
GTP->GSN + 3PI 
DGTP-> DG + 3PI 



izaa 



4245 



424Z 



7.2 Glycoprotein biosynthesis PATH:hsa00510 
, DPAGT1, DPAGT, UGAT, UAGT, 
! D11S366, DGPT, DPAGT2, GPT 
ALG5 
8813 DPMI 

1650 DDOST, OST, OST48, KIAA0115 

6184 RPN1 

6185 RPN2 
10130 P5 
10954 PDIR 
11008 PDI 

GRP58, ERp57, ERp60, ERp61, 
2223 GRP57, P58, PI- PLC, ERP57, ERP60, 
ERP61 

5Q34 P4HB. PROHB, P04DB, ERBA2L 
7841 GCS1 

4121 MAN1A1, MAN9, HUMM9 
MGAT1, GLYT1, GLCNAC-T7, GNT-l, 

! MGAT 

4122 MAN2A2, MANA2X 

4124 MAN2A1, MANA2 
, MGAT2, CDGS2, GNT-II, GLCNACTU, 
' GKT2 

4246 MGAT3, GNT-W 
6487 SIAT6, ST3GAUI 
6480 S1AT1 

^39 FNTA, FPTA, PGGT1A 
2342 FNTB, FPTB 
£229 PGGT1B, BGGI, GGTI 

5875 RABGGTA 

5876 RABGGTB 
1352 COX10 

7.3 Glycoprotein degradation PATH:hsa0051 1 

4758 NEU1, NEU 

3073 HEXA, TSD 

3074 HEXB 

4123 MAN2C1, MANA, MANA1 , MAM6A8 

4125 MAN2B1, MANB, LAMAN 

4126 MANBA, MANB1 
2517 FUCA1 
2519 FUCA2 • 

125 AGA.AGU 

7.4 Aminosugars metabolism PATH:hsa00530 

6675 UAP1, SPAG2, AGX1 
10020 GNE, GLCNE 
22951 CMAS 

1727 DIA1 

4669 NAGLU, NAG 

7.5 Upopolysaccharide biosynthesis PATH:hsa00540 

6485 SIAT5, SAT3, STZ 
7903 SIAT8D, PST, PST1, ST8SIA-IV 
8128 SIAT8B, STX, ST8SIA-II 
7.7 Glycosaminoglycan degradation PATH:hsaO0531 
3423 IDS, MPS2, SIDS 
3425 IDUA, IDA 
411 ARSB 
2799 GNS, G6S 

2588 GAUMS, MPS4A, GALNAC6S, GAS 
8. Metabolism of Complex Lipids 
8.1 Glycerofipid metabolism PATH:hsa00561 



GDPMAN + DOLP -> GDP + DOLMANP 



UTP + NAGA1P <-> UDPNAG + PP! 



3.6.1. - 



2.7.B.15 

2.4.1.117 

2.4.1. B3 

2.4.1.119 

2,4,1.119 

2.4,1,119 

5.3.4.1 

5-3.4.1 

5.3.4.1 



5.3.4.1 

5.3.4.1 

3.2.1.106 

3.2.1,113 

2.4.1.1Q1 

3.2.1.114 
3.2.1.114 

2.4.1.143 

2.4.1,1 44 

2.4.99.6 

2.4.99.1 

2.5.1 
2.5.1 

2A1= 

2.5.1. - 
2.5.1. - 
?.5.1 - 

3.2.1.18 
3,2.1.52 
3*2JL52_ 
3.2.1.24 
3.2X24 
3.2.1.25 
3.2.1.51 
3.2.1.51 
3.5.1.26 

2.7,7 .2 3 
5.1.3,14 

2.7.7.43 

1.6.2.2 

3.2.1.50 

2.4.99,- 

74.99.- 

3.1,6.13 
3.2.1 ,76 
3.1.6.12 
3.1,6.14 
3.1.6.4 



10554 AGPAT1, LPAAT-ALPHA, G15 



AGUP + 0.017 C100ACP + 0.062 C120ACP + 0.100 C140ACP h 
0.270 C160ACP + 0.169 C161ACP + 0.055 C180ACP + 0235 
C181ACP + 0.093 C182ACP -> PA + ACP 



7-3.1.51 
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on 






AGL3P +0.017 C100ACP + 0.062 C120ACP +0.100 C140ACP + 




10555 AGPAT2, LPAAT-BETA 


0.270 C160ACP + 0.169 C161ACP + 0.055 C180ACP + 0.235 


2.3.1.51 




C181ACP + 0.093 C182ACP -> PA + ACP 




IfiQfi DGKA, DAGK, DAGK1 




2.7,1.107 


16QB DGKG, DAGK3 




2.7.1,107 


1509 DGKQ, DAGK4 




2,7.1.107 


B525 DGKZ, DAGK5, HDGKZETA 




2.7.1.107 


3525 DGKE, DAGK6, DGK 




2.7.1.107 


B52Z DGKD, DGKDELTA. KIAA0145 




2.7.1.107 


112Q CHKL 


ATP+CHO->ADP+PCHO 


2.7.1.32 


EKH 


ATP + ETHM -> ADP + PETHM 


2.7.1.82 


ma chk, cki 


ATP + CHO -> ADP + PCHO 


2.7.1.32 


43 ACHE, YT 




3.1.1.7 


1103 CHAT 




2.3.1.6 


533Z PLD1 




3.1.4.4 


2S2ZS PLA2G2D, SPLA2S 




3.1.1.4 


30814 PLA2G2E 




3.1.1.4 


5319 PLA2G1B, PLA2, PLA2A, PPLA2 




3.1.1.4 


5320 PLA2G2A, M0M1, PIA2B, PLA2L 




3.1.1.4 


5222 PLA2G5 




3.1.1.4 


£398 PLA2G6, 1PLA2 




3.1.1.4 


3329 PLA2G10, SPLA2 




3.1.1.4 


1040 CDS1 


PA + CTP <-> CDPDG + PPl 


2.77.41 


10423 PIS 


CDPDG + MYOI -> CMP + PINS 


2.7.8.1 1 


2710 GK 


GL+ATP->GL3P + ADP 


2.7.1.30 


2320 GPD2 


GL3Pm + FADm -> T3P2m + FAD H 2m 


1.1.99.5 


frQ i j orU 1 


T?P? + NAHM <-> f5L3P + NAD 


1118^ 


OAft At Dl 


AHTH -!> HUP + ^ PI 
Mil \ U Unr T On 




249 ALPL, HOPS, TNSALP 


AHTD -> DHP + 3 PI 


3.1.3.1 


25Q ALPP 


AHTD -> DHP + 3 PI 


3.1.3.1 


251 ALPPL2 


AHTD -> DHP + 3 PI 


3.1.3.1 


439 ASNA1, ARSA-l 




3.6.1.16 








fUtQA nf5AT AR^P1 


+ 0-270 C160ACP + 0 169 C161ACP + 0 055 C180ACP + 0.235 


2.3:1.20 




C181ACP + 0.093 C182ACP -> TAGLY + ACP 




aaaauPB 




3.1.1.3 


3220 LIPC, HL 




3.1.1.3 


5406 PNUP 




3.1.1.3 


54QZ PNUPRP1, PLRP1 




3.1.1.3 


54QS PNLIPRP2, PLRP2 




3.1.1.3 


8513 LIPF. HGL, HLAL 




3.1.1.3 


4023 LPL, LIPD 




3.1.1.34 


8443 GNPAT, DHAPAT, DAP-AT 




2.3.1.42 



8540 



AGPS, ADAP-S, ADAS, ADHAPS, 



ADPS, ALDHPSY 
MDCR, MDS, US1 
5043 PAFAH1B1, US1, MDCR, PAFAH 

5049 PAFAH1B2 

5050 PAFAH1B3 

5051 PAFAH2, HSD-PLA2 
7941 PLA2G7, PAFAH, LDL-PLA2 

8.2 Inositol phosphate metabolism PATH:hsa00562 



2.5.1.26 

3.1.1.47 
3,1.1.47 
3.1.1.47 
3.1.1.47 
3.1.1.47 
3.1,1.47 



5220 PIK3CA 

5221 P1K3CB, PIK3C1 
5233 PIK3CD 

5294 PIK3CG 

5292 PIK4CA, PWK-ALPHA 

5305 PIP5K2A 

533Q PLCB2 

5331 PLCB3 

5333 PLCD1 

5335 PLCG1, PLC1 

5333 PLCG2 

3512 IMPA1, IMPA 

3513 IMPA2 
3623 INPP1 



ATP + PINS -> ADP + PINSP 
ATP + PINS -> ADP + PINSP 
ATP + PINS -> ADP + PINSP 
ATP + PINS •> ADP + PINSP 
ATP + PINS -> ADP + PINS4P 
PINS4P + ATP -> D45PI + ADP 
D45PI-> TPI + DAGLY 
D45PI -> TPI + PAGLY 
D45PI -> TPI + DAGLY 
D45PI -> TPI + DAGLY 
D45PI -> TPI + DAGLY 
MI1P->MYOI + PI 
MI1P->MY0I + PI 



2,7.1,137 

2.7.1.137 

2.7.1.137 

2.7.1.137 

2.7.1.67 

2.7,1.66 

3.1.4.11 

3.1.4.11 

3.1.4.11 

3.1.4.11 

3.1.4.11 

3*1A25_ 

3.1.3.25 

3.1.3.57 
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3632 INPP5A 

3S33 INPP5B 

353S INPPL1.SHIP2 

4252 OCRL, LOCR, OCRL1, INPP5F 

fififiZSYNJl,INPP5G 

3706 ITPKA 

51477 ISYNA1 G6P -> MI1P 

3631 INPP4A, INPP4 
8821 INPP4B 

8.3 SphingophosphoUpid biosynthesis PATH:hsa00570 

fifiQS SMPD1.NPD 

8.4 Phospholipid degradation PATH:hsa00580 

HZfiCLC 

5321 PLA2G4A, CPLA2-ALPHA, PLA2G4 

8.5 Sphingoglycoiipid metabolism PATH:hsa00600 

10558 SPTLC1, LCB1, SPTI PALCOA + SER -> COA + DHSPH + C02 

9517 SPTLC2, K1AA0526, LCB2 PALCOA + SER -> COA + DHSPH + C02 

427 ASAH, AC, PHP32 
7357 UGCG, GCS 
2629 GBA, GLUC 
2583 GALGT, GALNACT 
6489 SIAT8A, SIAT8, ST8SIA-1 
54S1SIAT2 

4668 NAGA, D22S674, GALB 
2514 CST 
410 ARSA, MLD 

8.6 Blood group glycolipid biosynthesis - lact series PATH:hsa00601 

23 ABO 

2525 FUT3, LE 
2S21 PUTS, FUC-TV 
2528 FUT6 

2523 FUT1, H, HH 

2524 FUT2, SE 

8.7 Blood group glycolipid biosynthesis - neolact series PATH:hsa00602 

2651 GCNT2, IGNT, NACGT1 , NAGCT1 

8.8 Prostaglandin and leukotriene metabolism PATH:hsa00590 

239 ALOX12, LOG12 
245 ALOX15 

240 ALOX5 
4056 LTC4S 
4048 LTA4H 

4051 CYP4F3, CYP4F, LTB4H 
8529 CYP4F2 

5742 PTGS1, PGHS-1 • 

5743 PTGS2, COX-2, COX2 
27306 PGDS 

5730 PTGDS 
5740 PTGIS, CYP8, PGIS 
6916 TBXAS1.CYP5 
873 CBR1.CBR 



874 CBR3 
9. Metabolism of Cofactors and Vitamins 

9.2 Riboflavin metabolism PATH:hsa00740 

52 ACP1 

FMN -> RIBOFLAV + PI 

53 ACP2 FMN -> RIBOFLAV + PI 

54 ACP5, TRAP FMN -> RIBOFLAV + PI 

55 ACPP, PAP FMN -> RIBOFLAV + PI 

9.3 Vitamin B6 metabolism PATH:hsa00750 

fi5£fi PDXK, PKH, PNK PYRDX + ATP -> P5P + ADP 

PDLA + ATP -> PDLA5P + ADP 
PL + ATP -> PL5P + ADP 

9.4 Nicotinate and nicotinamide metabolism PATH:hsa00760 

2 3475 QPRT QA + PRPP->NAMN + C02 + PPI 



3.1.3.56 

3.1.3.56 

3.1,3,56. 

3.1.3.56 

2.7.1.127 

5.5.1.4 

3.1.3.66 

3.1,3.66 

3.1,4.12 

3.1.1.5 
3.1 .1.5 

2,3.1.50 
2.3.1.50 

3.5.1.23 
2.4.1.80 
3.2.1.45 
2.4.1.92 
2.4.99.B 
2.4.99.2 
3.2.1.49 
2.8.2.11 
3.1.6.8 

2.4.1.40 
2,4.1.37 
2,4.1.65 
2.4.1.65 
2.4.1.65 
2 .4 .1.69 
2.4.1.69 

2 .4.1.150 

1.13.11.31 

1-13.11.33 

1.13.1 1 .3 4 

2,5.1.37 

3^2J_ 

1.14.13.30 

1. 14.13,30 

1-14.99.1 

1.14.99.1 

5.3.99.2 

5^29-2. 

5.3.99.4 

5&sa5_ 

1.1.1.134 
1,1.1.189 
1.1.1.197 
1.1.1.184 



3.1.3.48 

3.1.3.2 

3.1.3.2 

3.1.3.2 

3.1.3.2 

2.7.1.35 



2,4,2.19 
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4837 NNMT 

£33 BST1, CD157 NAD -> NAM + ADPRIB 

252 CD38 NAD -> NAM + ADPRIB 

2353Q NNT 

9.5 Pantothenate and CoA biosynthesis PATH:hsa00770 

9.6 Biotin metabolism PATH:hsaO078O 

3141 HLCS. HCS 



GTP -> FOR + AHTD 
DHF + NADPH -> NADP + THF 
THF + ATP + GLU <-> ADP + PI + THFG 



NADPH + DHBP -> NADP + THBP 



ATP + FTHF -> ADP + PI + MTHF 



fiSS BTD 

9 J Folate biosynthesis PATH:hsa00790 
2643 GCH1, DYT5, GCH, GTPCH1 
iZlSDHFR 
2256 FPGS 
£336 GGH, GH 
5SQ5 PTS 
fifiSZSPR 

5SBQ QDPR, DHPR, PKU2 
9.8 One carbon pool by folate PATH:hsa00670 
10840 FTHFD 
10588 MTHFS 

9.10 Porphyrin and chlorophyll metabolism PATH:hsa00860 

21fl ALAD 2 ALAV -> PBG 

3145 HMBS, PBGD, UPS 4 PBG -> HMB + 4 NH3 

7390 UROS HMB -> UPRG 

2239 UROD UPRG -> 4 C02 + CPP 

1371 CPO, CPX 02 + CPP -> 2 C02 + PPHG 

5498 PPOX, PPO 02 + PPHGm -> PPIXm 

2235 FECH, FCE PPIXm ± PTHm 

2162 HMOX1.HO-1 

2163 HMOX2, HO-2 

644 BLVRA, BLVR 

645 BLVRB, FLR 



2232 FDXR, ADXR 
3JJ52 HCCS, CCHL 
1356 CP 

9.11 Ubiquinone biosynthesis PATH:hsa00130 

493SOAS1.IFI-4.0IAS 

4939 OAS2, P69 

555Z PRIM1 

£55fi PRIM2A. PR1M2 

5553 PRIM2B. PRIM2 

701 S TERT. EST2, TCS1» TP2, TRT 

5523 OASL, TRIP14 
10. Metabolism of Other Substances 

10.1 Terpenoid biosynthesis PATH:hsa00900 

10.2 Flavonoids, stflbene and lignin biosynthesis PATH:hsa00940 

10.3 Alkaloid biosynthesis I PATH.tisa00950 

10.4 Alkaloid biosynthesis II PATH:hsa00960 

10.6 Streptomycin biosynthesis PATH:hsa00521 

10.7 Erythromycin biosynthesis PATH:hsa00522 

10.8 Tetracycline biosynthesis PATH:hsa00253 

10.14 gamma-Hexachlorocyclohexane degradation PATH:hsa00361 

5444 PON1.ESA, PON 

5445 PON2 



10.18 1,2-Dichloroethane degradation PATH:hsa00631 

10.20 Tetrachloroethene degradation PATH:hsa00625 
2Q52 EPHX1,EPHX, MEH 

2Q53 EPHX2 

10.21 Styrene degradation PATH:hsa00643 
11. Transcription (condensed) 

11.1 RNA polymerase PATH:hsa03020 



2.1.1.1 
3^5_ 
3.2.2.5 
1.6.1.2 



6.3.4.- 

6.3.4.9 

6 .3. 4.10 

.6.3.4.11 

6.3,4.15 

3.5.1.12 

3.5.4.16 

1.5.1.3 

6.3.2.17 

3.4.19,9 

4.6,1.10 

1.1,1.153 

.1.6.99.7 

1.5.1.6 
6.3.3.2 

4.2.1.24 

4.3.1.8 

4.2.1,75 

4.1.1.37 

1.3.3.3 

1.3.3.4 

4.99.1.1 

1-14.99.3 

1-14.99.3 

1-3.1.24 

1.3.1.24 

1.6.99.1 

1.16.1.2 

4.4.1,17 

.1.16.3.1 

2.7.7.- 
2.7.7.- 
2.7.7.- 
2.7.7.- 
2.7,7.- 
2.7.7.- 
2.7,7.- 



3.1.8.1 
3.1.1.2 
3.1.1.2 
3.1.8.1 



3-3.2.3 
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11.2 Transcription factors PATH:hsa03022 

12. Translation (condensed) 
12.1 Ribosome PATH:hsa03010 

12J2 Translation factors PATH:hsa03012 
. EEF1A1, EF1A. ALPHA, EEF-1, 
' EEF1A 
JL21I EEF1A2, EF1A 
1938 EEF2, EF2, EEF-2 

12.3 Aminoacyl-tRNA biosynthesis PATH:hsa00970 

13. Sorting and Degradation (condensed) 
13.1 Protein export PATH:hsa03060 

234ZSSPC18 

13.4 Proteasome PATH:hsa03050 

5687 PSMA6, IOTA PROS27 

5683 PSMA2, HC3. Ml), PMSA2, PSC2 

5685 PSMA4, HC9 

5688 PSMA7, XAPC7 

5686 PSMA5, ZETA, PSC5 
5682 PSMA1, HC2, NU, PROS30 

5684 PSMA3, HC8 

5698 PSMB9, LMP2, RING12 
5695 PSMB7, Z 
5691 PSMB3, HC1CM1 
5690 PSMB2. HC7-I 
5693 PSMB5, LMPX, MB1 
55S9 PSMB1, HC5. PMSB1 
5622 PSMB4, HN3, PROS26 

14. Replication and Repair 

14.1 DNA polymerase PATH:hsa03030 

14.2 Repncation Complex PATH:hsa03032 
23626 SP011 

7153 TOP2A, TOP2 

7155 f OP2B 

7156 TOP3A, TOP3 
8940 TOP3B 

22. Enzyme Complex 

22.1 Electron Transport System, Complex I PATH:hsa03100 

22.2 Electron Transport System, Complex II PATH:hsa03150 

22.3 Electron Transport System, Complex 111 PATH:hsa03140 

22.4 Electron Transport System, Complex IV PATH:hsa03130 

22.5 ATP Synthase PATH:hsa031 10 
22.8 ATPases PATH:hsa03230 

23. Unassigned 

23.1 Enzymes 
5538 PPT1, CLN1, PPT, INCL 

23.2 Non-enzymes 
22934 RPIA, RPI 

5250 SLC25A3, PHC 



C160ACP + H20 -> C160 + ACP 



51166 LOC51 166 
5625 PRODH 
6517SLC2A4, GLUT4 
6513 SLC2A1, GLUT1, GLUT 

26275 HIBCH, HIBYL-COA-H 

23305 KIAA0837, ACS2, LACS5, LACS2 

8611 PPAP2A PAP-2A 

8612 PPAP2C, PAP-2C 

8613 PPAP2B, PAP-2B 
56994 LOC56994 
10400 PEMT, PEMT2 

5833 PCYT2, ET 
10390 CEPT1 

8394 PIP5K1A 

8395 P1P5K1B, STM7. MSS4 



RL5P <-> R5P 
PI + H <-> Hm + Pirn 
CIT + MALm <-> CITm + MAL 
AADP + AKG -> GLU + KADP 
PRO + FAD -> P5C + FADH2 
GLCxt->GLC 
GLCxt-> GLC 

HIBCOAm + H20m -> HIBm + COAm 
C160 + COA + ATP -> AMP + PPI + C160COA 
PA+H20->DAGLY + PI 
PA+H20->DAGLY + Pl 
PA + H20 -> DAGLY + PI 
CDPCHO + DAGLY -> PC + CMP 
SAM + PE -> SAH + PMME 
PETHM + CTP-> COPETN + PPI 
CDPETN + DAGLY <-> CMP + PE 
PINS4P + ATP -> D45P1 + ADP 
PINS4P + ATP -> D45PI + ADP 
PINS4P + ATP -> D45P1 + ADP 



3 ,5 ,1,48 , 

3.6.1.46 

3.6.1.48 



3 4.21.89 

3-4.99.46 
3.4.99.46 
3.4.99.46 
,3.4,99.46 
3.4.99.46 
3.4. 99 .46 
3.4.99.46 
3.4.99.46 



3.4.99.46 
3.4.99.46 
3.4.99.46 
3,4,99,46 

3.4.99.46 



5,99.1.3 
5,99.1,3 
5.99,1,3 

5.99.1.2 
5.99.1.2 



3.1.2.22 

5.3.1.6 



2.6.1.39 
1.5.3." 



3.1.2.4 



8396 PIP5K2B 

23396 PIP5K1C, KIAA0589, PIP5K-GAMMA PINS4P + ATP •> D45PI + ADP 
24. Our own reactions which need to be found In KEGG 
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amino acid metabolism 



GL3P <-> GL3Pm 
T3P2 <->T3P2m 
PYR <-> PYRm + Hm 

ADP + ATPm + PI + H -> Hm + ADPm + ATP + Ptm 

AKG + MALm <-> AKGm + MAL 

ASPm + GLU + H -> Hm + GLUm + ASP 

GDP + GTPm + PI + H -> Hm + GDPm + GTP + Pirn 

C160Axt + FABP -> C160FP + ALBxt 

C160FP -> C160 + FABP 

C180Axt + FABP -> C180FP + ALBxt 

C180FP -> C180 + FABP 

C161Axt + FABP -> C161FP + ALBxt 

C161FP->C161 + FABP 

C181Axt + FABP -> C181FP + ALBxt 

C181FP-> C181 + FABP 

C182Axt + FABP -> C182FP + ALBxt 

C182FP->C182 + FABP 

C204Axt + FABP -> C204FP + ALBxt 

C204FP -> C204 + FABP 

02xt -> 02 

02<-> 02m 

ACTACm + SUCCOAm -> SUCCm + AACCOAm 
3HB -> 3HBm 

MGCOAm + H20m -> H3MCOAm 
OMVAL -> OMVALm 
OIVAL -> OIVALm 
OICAP -> OICAPm 
C160CAR <-> C160CARm 
CAR <-> CARm 
DMMCOAm -> LMMCOAm 

THR -> NH3 + H20 + OBUT 
THR + NAD -> C02 + NADH + AMA 
THR + NAD + COA -> NADH + ACCOA + GLY 
AASA + NAD -> NADH + AADP 
FKYN + H20 -> FOR + KYN 
CMUSA -> C02 + AM6SA 
AM6SA + NAD -> AMUCO + NADH 
AMUCO + NADPH -> KADP + NADP + NH4 
CYSS + AKG <-> GLU + SPYR 
URO + H20->4l5P 
4I5P + H20 -> FIGLU 
GLU <-> GLUm + Hm 
. ORN + Hm -> ORNm 
ORN + Hm + CITRm <-> CITR + ORNm 
GLU + ATP + NADPH -> NADP + ADP + PI + GLUGSAL 
GLYAm + ATPm -> ADPm + 2PGm 



4.2.1.18 



5.1.99.1 

4.2.1.16 
1.T.1.1Q3 

1.2.1.31 
3.5.1.9 

4.1.1.45 

1-2.1.32 

1.5.1.- 

4.2.1.49 
&5J2J 



fatty acid synthesis 



AM6SA->PlC 

SPYR + H20 -> H2S03 + PYR 
P5C <-> GLUGSAL 

MALCOA + ACP <-> MALACP + COA 
ACCOA + ACP <-> ACACP + COA 

ACACP + 4 MALACP + 8 NADPH -> B NADP + C100ACP + 4 
C02 + 4ACP 

ACACP + 5 MALACP + 10 NADPH -> 10 NADP + C 120 ACP + 5 
C02 + 5ACP 

ACACP + 6 MALACP + 12 NADPH -> 12 NADP + C140ACP + 6 
C02 + 6ACP 

ACACP + 6 MALACP + 11 NADPH -> 11 NADP + C141ACP + 6 
C02 + 6 ACP 

ACACP + 7 MALACP + 14 NADPH -> 14 NADP + C160ACP + 7 
C02 + 7ACP 

ACACP + 7 MALACP + 13 NADPH -> 13 NADP + C161ACP + 7 
C02 + 7 ACP 



2.3.1.39 
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fatty acid degradation 



Phospholipid metabolism 



Muscle contraction 



ACACP + 8 MALACP + 16 NADPH -> 16 NADP + C18QACP + 8 
C02 + 8ACP 

ACACP + 8 MALACP + 15 NADPH -> 15 NADP + C181ACP + 8 
C02 + 8ACP 

ACACP + 8 MALACP + 14 NADPH -> 14 NADP + C182ACP + 8 
C02 + 8ACP 

C160COA + CAR -> C160CAR + COA 
C160CARm + COAm -> C160COAm + CARm 

GL3P + 0.017 C100ACP + 0.062 C120ACP + 0.1 C140ACP + 
0.27 C160ACP + 0.169 C161ACP + 0.055 C180ACP + 0.235 
C181ACP + 0.093 C182ACP -> AGL3P + ACP 
TAGLYm + 3 H20m -> GLm + 3 C160m 

SAM + PMME -> SAH + PDME 
PDME + SAM -> PC + SAH 
PE + SER<->PS + ETHM 

MYOACT + ATP -> MYOATP + ACTIN 

MYOATP + ACTIN -> MYOADPAC 

MYOADPAC -> ADP + PI + MYOACT + CONTRACT 
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Table 2 

// Homo Sapiens Core Metabolic Network // 

i 

// Glycolysis // 

-1 GLC -1 ATP +1 G6P +1 ADP 0 HK1 
-1 G6P -1 H20 +1 GLC +1 PI 0 G6PC 
-1 G6P +1 F6P 0 GPIR 
-1 F6P -1 ATP +1 FDP +1 ADP 0 PFKL 
-1 FDP -1 H20 +1 F6P +1 PI 0 FBP1 
-1 FDP +1 T3P2 +1 T3P1 0 ALDOAR 
-1 T3P2 +1 T3P1 0 TPI1R 

-1 T3P1 -1 PI -1 NAD +1 NADH +1 13PDG 0 GAPDR 

-1 13PDG -1 ADP +1 3PG +1 ATP 0 PGK1R 

-1 13PDG +1 23PDG 0 PGAM1 

-1 23PDG -1 H20 +1 3PG +1 PI 0 PGAM2 

-1 3PG +1 2PG 0 PGAM3R 

-1 2PG +1 PEP +1 H20 0 ENOIR 

-1 PEP -1 ADP +1 PYR +1 ATP 0 PKLR 

-1 PYRm -1 COAm -1 NADm +1 NADHm +1 C02m +1 ACCOAm 0 PDHA1 
-1 NAD -1 LAC +1 PYR +1 NADH 0 LDHAR 
-1 G1P +1 G6P 0 PGM1R 

// TCA // 

-1 ACCOAm -1 OAm -1 H20m +1 COAm +1 CITm 0 CS 
-1 CIT +1 ICIT 0 ACOIR 
-1 CITm +1 ICITm 0 AC02R 

-1 ICIT -1 NADP +1 NADPH +1 C02 +1 AKG 0 IDH1 

-1 ICITm -1 NADPm +1 NADPHm +1 C02m +1 AKGm 0 IDH2 

-1 ICITm -1 NADm +1 C02m +1 NADHm +1 AKGm 0 IDH3A 

-1 AKGm -1 NADm -1 COAm +1 C02m +1 NADHm +1 SUCCOAm 0 OGDH 

-1 GTPm -1 SUCCm -1 COAm +1 GDPm +1 Pirn +1 SUCCOAm 0 SUCLG1R 

-1 ATPm -1 SUCCm -1 COAm +1 ADPm +1 Pirn +1 SUCCOAm 0 SUCLA2R 

-1 FUMm -1 H20m +1 MALm 0 FHR 

-1 MAL -1 NAD +1 NADH +1 OA 0 MDH1R 

-1 MALm -1 NADm +1 NADHm +1 OAm 0 MDH2R 

-1 PYRm -1 ATPm -1 C02m +1 ADPm +1 OAm +1 PIm 0 PC 

-1 OA -1 GTP +1 PEP +1 GDP +1 C02 0 PCK1 

-1 OAm -1 GTPm +1 PEPm +1 GDPm +1 C02m 0 PCK2 

-1 ATP -1 CIT -1 COA -1 H20 +1 ADP +1 PI +1 ACCOA +1 OA 0 

ACLY 
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// PPP // 

-1 G6P -1 NADP +1 D6PGL +1 NADPH 0 G6PDR 

-1 D6PGL -1 H20 +1 D6PGC 0 PGLS _ • 

-1 D6PGC -1 NADP +1 NADPH +1 C02 +1 RL5P 0 PGD 

-1 RL5P +1 X5P 0 RPER 

-1 R5P -1 X5P +1 T3P1 +1 S7P 0 TKT1R 

-1 X5P -1 E4P +1 F6P +1 T3P1 0 TKT2R 

-1 T3P1 -1 S7P +1 E4P +1 F6P 0 TALD01R 

-1 RL5P +1 R5P 0 RPIAR 

// Glycogen // 

-1 G1P -1 UTP +1 UDPG +1 PPI 0 UGP1 
-1 UDPG +1 UDP +1 GLYCOGEN 0 GYS1 
-1 GLYCOGEN -1 PI +1 G1P 0 GBE1 

I i // 

-1 MALra -1 NADPm +1 C02m +1 NADPHm +1 PYRm 0 ME 3 

-1 MALm -1 NADm +1 C02m +1 NADHm +1 PYRm 0 ME2 

-1 MAL -1 NADP +1 C02 +1 NADPH +1 PYR 0 ME1 

-1 NADHm -1 Qm -4 Hm +1 QH2m +1 NADm +4 H 0 MTND1 

-1 SUCCm -1 FADm +1 FUMm +1 FADH2m 0 SDHC1R 

-1 FADH2m -1 Qm +1 FADm +1 QH2m 0 SDHC2R 

-1 02m -4 FEROm -4 Hm +4 FERIm +2 H20m +4 H 0 UQCRFS1 

-1 QH2m -2 FERIm -4 Hm +1 Qm +2 FEROm +4 H 0 COX5BL4 

-1 ADPm -1 PIm -3 H +1 ATPm +3 Hm +1 H20m 0 MTAT 

-1 ADP -1 ATPm -1 PI -1 H +1 Hm +1 ADPm +1 ATP +1 PIm 0 AT PMC 

-1 GDP -1 GTPm -1 PI -1 H +1 Hm +1 GDPm +1 GTP +1 PIm 0 GTPMC 

-1 PPI +2 PI 0 PP 

-1 ACCOA -1 ATP -1 C02 +1 MALCOA +1 ADP +1 PI 0 ACACAR 
-1 GDP -1 ATP +1 GTP +1 ADP 0 GOT3R 



// Transporters // 

-1 CIT -1 MALm +1 CITm +1 MAL 0 CITMCR 
-1 PYR -1 H +1 PYRm +1 Hm 0 PYRMCR 

// Glycerol Phosphate Shuttle // 

-l'GL3Pm -1 FADm +1 T3P2m +1 FADH2m 0 GPD2 

-1 T3P2 -1 NADH +1 GL3P +1 .NAD 0 GPD1 

-1 GL3P +1 GL3Pm 0 GL3PMCR 

-1 T3P2 +1 T3P2m 0 T3P2MCR 

// Malate/Aspartate Shuttle // 

-1 OAm -1 GLUm +1 ASPm +1 AKGm 0 GOT1R 

-1 ASP -1 AKG +1 OA +1 GLU 0 GOT2R 

-1 AKG -1 MALm +1 AKGm +1 MAL 0 MALMCR 

-1 ASPm -1 GLU -1 H +1 Hm +1 GLUm +1 ASP 0 ASPMC 
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// Exchange Fluxes // 

+1 GLC 0 GLCexR 

+1 PYR 0 PYRexR 

+1 C02 0 C02exR 

+1 02 0 02exR 

+1 PI 0 PlexR 

+1 H20 0 H20exR 

+1 LAC 0 LACexR 

+1. C02m 0 C02min 
-1 C02m 0 C02mout 
+1 02m 0 02min 
-1 02m 0 02mout 
+1 H20m 0 H20min 
-1 H20m 0 H20mout 
+1 PIm 0 PImin 
-1 PIm 0 PImout 



// Output // 

-1 ATP +1 ADP +1 PI 0 Output 
0.0 end 
end E 0 
max 

1 Output 
0 end 

0 GLCexR 1 
-1000- PYRexR 0 
-1000 LACexR 0 

0 end 0 
rev. rxn 33 
nonrev. rxn 31 
total rxn 64 
matrix columns 97 
unique enzymes 52 
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Table 3 

Abbrev. Reaction 
Glycolysis 
HK1 

G6PC, G6PT 
GPI 
PFKL 
FBP1.FBP 
ALDOA 
TPI1 

GAPD, GAPDH 
PGK1, PGKA 
PGAM1, PGAMA 



EN01.PPH. EN01L1 
PKLR, PK1 

PDHA1, PHE1A, PDHA 
LDHA, LDH1 
PGM1 
TCA 
CS 

AC01, IREB1, IRP1 
AC02 
IDH1 
IDH2 
IDH3A 
OGDH 

SUCLG1 , SUCLA1 
SUCLA2 
FH 

MDH1 
MDH2 
PC, PCB 

ACLY, ATPCL, CLATP 
PCK1 
PPP 

G6PD. G6PD1 
PGLS, 6PGL 
PGD 
RPE 
TKT 

TALD01 
UGP1 

ACACA, ACAC, ACC 
ETS 
ME3 
MTND1 
SDHC 

UQCRFS1, R1S1 
COX5BL4 
MTATP6 
PP, SID6-8061 
Matate Aspartate shunttie 
GOT1 
GOT2 



GLC + ATP->G6P+ADP 
G6P + H20->GLC + PI 
G6P <-> F6P 
F6P+ATP->FDP + ADP 
FDP + H20->F6P + PI 
FDP<->T3P2 + T3P1 
T3P2^>T3P1 

T3P1 + PI + NAD <-> NADH + 13PDG 
13PDG + ADP <-> 3PG + ATP 
13PDG<~>23PDG 
23PDG + H20 -> 3PG + PI 
3PG <-> 2PG 
2PG <-> PEP + H20 
PEP + ADP -> PYR + ATP . 

PYRm + CO Am + NADm -> + NADHm + C02m + ACCOAm 
NAD + LAC <-> PYR + NADH 
G1P<->G6P 



ACCOAm + OAm + H20m -> COAm + CITm 

CIT<->ICIT 

CITm <-> lCITm 

ICIT + NADP -> NADPH + C02 + AKG 

ICITm + NADPm -> NADPHm + C02m + AKGm 

ICITm + NADm -> C02m + NADHm + AKGm 

AKGm + NADm + COAm -> C02m + NADHm + SUCCOAm 

GTPm + SUCCm + COAm <-> GDPm + Pirn + SUCCOAm 

ATPm + SUCCm + COAm <-> ADPm + Pirn + SUCCOAm 

FUMm + H20m <-> MALm 

MAL + NAD <-> NADH + OA 

MALm + NADm <■> NADHm + OAm 

PYRm + ATPm + C02m -> ADPm + OAm + Pirn 

ATP + C!T + COA + H20 ^> ADP + PI + ACCOA + OA 

OA + GTP -> PEP + GDP + C02 

G6P + NADP <-> D6PGL + NADPH 

D6PGL + H20 -> D6PGC 

D6PGC + NADP r> NADPH + C02 + RL5P 

RL5P<->X5P 

R5P + X5P<->T3P1 +S7P 

X5P + E4P <-> F6P + T3P1 

T3P1 + S7P <-> E4P + F6P 

G1P + IH*P->UDPG + PPI 

ACCOA + ATP + C02 <-> MALCOA + ADP + PI + H 

MALm + NADPm -> C02m + NADPHm + PYRm 
NADHm + Qm + 4 Hm -> QH2m + NADm + 4 H 
SUCCm + FADm <■> FUMm + FADH2m 
FADH2m + Qm <-> FADm + QH2m 
02m + 4 FEROm + 4 Hm -> 4 FERlm + 2 H20m + 4 H 
QH2m + 2 FERlm + 4 Hm ■> Qm + 2 FEROm + 4 H 
ADPm + Pirn + 3 H -> ATPm + 3 Hm + H20m 
PPl -> 2 PI 

OAm + GLUm <-> ASPm + AKGm 
OA + GLU <-> ASP + AKG 
GDP + ATP <-> GTP + ADP 



RxnName 

HK1 

G6PC 

GPI 

PFKL 

FBP1 

ALDOA 

TPI1 

GAPD 

PGK1 

PGAM1 

PGAM2 

PGAM3 

EN01 

PKLR 

PDHA1 

LDHA 

PGM1 

CS 

AC01 

AC02 

IDH1 

IDH2 

IDH3A 

OGDH 

SUCLG1 

SUCLA2 

FH 

MDH1 

MDH2 

PC 

ACLY 

PCK1 

G5PD 

PGLS 

PGD 

RPE 

TKT1 

TKT2 

TALD01 

UGP1 

ACACA 

ME3 

MTND1 

SDHC1 

SDHC2 

UQCRFS1 

C0X5BL4 

MTAT 

PP 

G0T1 
G0T2 
G0T3 
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Glycogen 

GBE1 GLYCOGEN + PI -> G1P GBE1 

GYS1 , GYS UDPG -> UDP + GLYCOGEN GYS1 

Glycerol Phosphate Shunttie 

GPD2 GL3Pm + FADm -> T3P2m + FADH2m GPD2 

GPD1 T3P2 + NADH->GL3P + NAD GPD1 

RPIA, RPI RL5P<->R5P RPIA 

Mitochondria Transport 

CIT + MALm <-> Cljm + MAL CITMC 

GL3P <-> GL3Pm GL3PMC 

T3P2 <-> T3P2m T3P2MC 

PYR <-> PYRm + Hm PYRMC 

ADP + ATPm + PI + H -> Hm + ADPm + ATP + Plm ATPMC 

AKG + MALm <-> AKGm + MAL MALMC 

ASPm + GLU + H -> Hm + GLUm + ASP ASPMC 

GDP + GTPm + PI + H -> Hm + GDPm + GTP + Plm GTPMC 
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TABLE 4 

Metabolic Reaction for Muscle Cells 



GLC + ATP-> G6P + ADP 
G6P<-»F6P 
F6P + ATP -> FDP + ADP 
FDP + H20->F6P + P1 
FDP <-> T3P2 + T3P1 
T3P2<->T3P1 

T3P1 + PI + NAD <-> NADH + 13PDG 
13PDG + ADP<-> 3PG + ATP 
3PG <-> 2PG 
2PG <-> PEP + H20 
PEP + ADP -> PYR + ATP 

PYRm + COAm + NADm -> + NADHm + C02m + ACCGAm 
NAD + LAC <-> PYR + NADH 
G1P<->G6P 

ACCOAm + OAm + H20m -> COAm + CITm 

crr<-> icrr 
CITm «■> icrrm 

1CTT + NADP -> NADPH + C02 + AKG 
iOTm + NADPm -> NADPHm + C02m + AKGm 
ICTTm + NADm ■> C02m + NADHm + AKGm 
AKGm + NADm + COAm •> C02m *■ NADHm + SUCCOAm 
GTPm + SUCCm + COAm <•> GDPm + Pirn + SUCCOAm 
ATPm + SUCCm + COAm o ADPm + Plm + SUCCOAm 
RJMm + H20m <-> MALm 
MAL + NAD <■> NADH + O A 
MALm + NADm <-> NADHm + OAm 
PYRm + ATPm + C02m -> ADPm + OAm + Plm 
ATP + CIT + CO A + H20 •> ADP + PI + ACCOA + OA 

OA + GTP -> PEP + GDP + C02 
OAm + GTPm -> PEPm + GDPm + C02m 
G6P + NADP <■> D6PGL + NADPH 
D6PGL+H20->D6PGC 
D6PGC + NADP -> NADPH + C02 + RL5P 
RL5P<~>X5P 
R5P + X5P <-> T3P1 + S7P 
X5P + E4P <->F6P + T3P1 
T3P1 + S7P <-> E4P + F6P 
RL5P<->R5P 
G1P + UTP •> UDPG + PPI 
GLYCOGEN + PI -> G1P 

UDPG -> UDP + GLYCOGEN 

MALm + NADm ■> C02m + NADHm + PYRm 

MALm + NADPm -> C02m + NADPHm + PYRm 

MAL + NADP -> C02 + NADPH + PYR 

NADHm + Can + 4 Hm -> QH2m + NADm + 4 H 

SUCCm + FADm <-> FUMm + FADH2m 

FADH2m + Qm <-> FADm + QH2m 

02m + 4 FEROm + 4 Hm -> 4 FERIm + 2 H20m + 4 H 

QH2m + 2 FERIm + 4 Hm -> Qm + 2 FEROm + 4 H 

ADPm + Plm + 3 H -> ATPm + 3 Hm + H20m 

ADP + ATPm + PI + H -> Hm + ADPm + ATP + Plm 

GDP + GTPm + PI + H •> Hm + GDPm + GTP + Plm 

PPI -> 2 PI 

GDP + ATP <-> GTP + ADP 
ACCOA + ATP + C02 <-> MALCOA + ADP + PI + H 
MALCOA + ACP <-> MALACP + COA 
ACCOA + ACP <^> ACACP + COA 

ACACP + 4 MALACP + 8 NADPH -> 8 NADP + C100ACP + 4 C02 + 4 ACP 
ACACP + 5 MALACP + 10 NADPH -> 10 NADP + C12QACP + 5 C02 + 5 
ACP 

ACACP ♦ 6 MALACP + 12 NADPH -> 12 NADP + C140ACP + 6 C02 + 6 
ACP 

ACACP + 6 MALACP + 1 1 NADPH -> 1 1 NADP + C141 ACP + 6 C02 + 6 



ACP 



ACACP + 7 MALACP + 14 NADPH -> 14 NADP + C160ACP + 7 C02 + 7 



ACP 



ACACP + 7 MALACP + 13 NADPH -> 13 NADP + C161 ACP + 7 C02 + 7 



ACP 



ACACP + 8 MALACP + 16 NADPH -> 16 NADP + C18QACP + 8 C02 + 8 



ACP 



ACACP + 8 MALACP + 15 NADPH -> 15 NADP + C181ACP + 8 C02 + 8 
ACP 

ACACP + 8 MALACP + 14 NADPH -> 14 NADP + C182ACP + 8 C02 + 8 
ACP 

C160ACP + H20 -> C160 + ACP 

C160 + COA + ATP*> AMP + PPI + C160COA 



RxtName 
0HK1 
0GPI 
0PFKL1 
0FBP1 
0ALDOA 
0TPI1 
0GAPD 
0PGK1 
0PGAM3 
0BIO1 
0PK1 
0PDHA1 
0LDHA 
0PGM1 
OCS 
0ACO1 
0ACO2 
0IDH1 
01DH2 
0IDH3A 
0OGDH 
0SUCLG1 
0 SUCLA2 
0FH 
0MDH1 
0MDH2 
0PC 
0ACLY 
0PCK1 
0PCK2 
0G6PD 
0H6PD 
0PGD 
ORPE 
0TKT1 
0TKT2 
0TALDO1 
0RPIA 
0UGP1 
0GBE1 
0GYS1 
0ME2 
0ME3 

0 HUMNDME 
0MTND1 
0 SDHC1 
0SDHC2 
0 UQCRFS1 
0 COX5BL4 
0MTAT1 
0 ATPMC 
0GTPMC 
0PP 
0NME1 
0ACACA 
0 FAS1JI 
0 FAS1_2 

0C100SY 

0 C120SY 

0 C140SY 

0 C141SY 

0C160SY 

0 C161SY 

0 C180SY 

0 C181SY 



0 C182SY 
0 PPT1 
0KIAA 
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C160COA + CAR-c C160CAR + COA 0 C160CA 

C160CARm + COAm •> C160COAm + CARm 0 C160CB 
C160CARm + COAm + FADm + NADm -> FADH2m + NADHm + 

C140COAm + ACC0Am 0 HADHA 
C140COAm + 7 COAm + 7 FADm + 7 NADm -> 7 FADH2m + 7 NADHm + 7 

ACCOAm 0HADH2 

TAGLYm + 3 H20m -> GLm + 3 C160m 0 TAGRXN 
GL3P + 0.017 C10DACP + 0.062 C120ACP + 0.1 C140ACP + 0.27 
C160ACP + 0.169 C161ACP + 0.055 C180ACP + 0.235 C181ACP + 0.093 

C182ACP->AGL3P + ACP 0GAT1 
AGL3P + 0.017 C100ACP + 0.062 C120ACP + 0.100 C140ACP + 0.270 
C160ACP + 0.169 C161ACP + 0.055 C1B0ACP + 0.235 C181ACP + 0.093 

C182ACP -> PA + ACP 0 AGPAT1 

ATP + CHO -> ADP + PCHO 0 CHKL1 

PCHO + CTP -> CDPCHO + PPI 0PCYT1A 

CDPCHO + DAGLY -> PC + CMP 0 LOC 

SAM+PE->SAH + PMME 0PEMT 

SAM + PMME -> SAH + PDME 0MFPS 

PDME + SAM -> PC + SAH 0PNMNM 

GSPo-MHP 0ISYNA1 

M!1P->MYOI + PI 0IMPA1 

PA + CTP«w CDPDG + PPI 0CDS1 

CDPDG + MYOI •> CMP + PINS 0PIS 

ATP + PINS -> ADP + PINSP 0 PIK3CA 

ATP + PINS -> ADP + PINS4P 0PIK4CA 

PINS4P + ATP -> D45PI + ADP 0 PIP5K1 

D45PI -> TPi + DAGLY 0 PLCB2 

PA + H20->DAGLY+PI 0PPAP2A 
DAGLY + 0.017 C100ACP + 0.062 C120ACP + 0.100 C140ACP + 0.270 
C160ACP + 0.169 C161ACP + 0.055 C180ACP + 0.235 C181ACP + 0.093 

C182ACP->TAGLY + ACP - 0 DGAT 

CDPDG + SER <-> CMP + PS 0 PTDS 

CDPETN + DAGLY o CMP + PE 0CEPT1 

PE + SER <-> PS + ETHM 0PESER 

ATP + ETHM -> ADP + PETHM 0 EK11 

PETHM + CTP-> CDPETN + PPI 0 PCYT2 

PS->PE + C02 0PISD 

3HBm + NADm -> NADHm + Hm + ACTACm 0BDH 

ACTACm + SUCCOAm -> SUCCm + AACOAm 0 30CT 

THF + SER <-> GLY + METTHF 0SHMT1 

THFm + SERm <-> GLYm + METTHFm 0 SHMT2 

SERm + PYRm <-> ALAm + 3HPm 0 AGXT 

3PG + NAD <-> NADH + PHP 0 PHGDH 

PHP + GLU <->AKG + 3PSER 0 PSA 

3PSER+H20->PI + SER 0 PSPH 

3HPm + NADHm -> NADm + GliYAm 0 GLYD 

SER-> PYR + NH3 + H20 0SDS 

GLYAm + ATPm ■> ADPm + 2PGm 0GLTK 

PYR + GLU <-> AKG + ALA 0GPT 

GLUm + C02m + 2 ATPm -> 2 ADPm + 2 Plm + CAPm 0CPS1 

AKGm + NADHm + NH3m <-> NADm + H20m + GLUm 0GLUD1 

AKGm + NADPHm + NH3m <-> NADPm + H20m + GLUm 0 GLUD2 

GLUm + NH3m + ATPm •> GLNm + ADPm + Plm 0GLUL 

ASPm + ATPm + GLNm -> GLUm + ASNm + AMPm + PPIm 0 ASNS 

ORN + AKG <-> GLUGSAL + GLU 0OAT 

GLU <-> GLUm + Hm 0 GLUMT 

GLU + ATP + NADPH -> NADP + ADP + PI + GLUGSAL 0 P5CS 
GLUP + NADH -> NAD + PI + GLUGSAL i = = 0 RYCS 

P5C <-> GLUGSAL 0 SPTC 

H1SO.NH3 + URO 0HAL 

UR0 + H20->4I5P 0 UROH 

4I5P + H20->FIGLU 0 IMPR 

F1GLU + THF -> NFTHF * GLU 0 FTCD 

MET + ATP + H20 -> PPI + PI + SAM 0 MAT1A 

SAM + DNA -> SAH + DNA5MC 0 DNMT1 

SAH + H20 -> HCYS + ADN 0AHCYL1 

HCYS + MTHF-> THF + MET 0MTR 

SER + HCYS -> LLCT + H20 0CBS 

LLCT + K20-> CYS + HSER 0 CTH1 

OBUT + NH3<->HSER 0 CTH2 

CYS + 02 <-> CYSS 0CDO1 

CYSS + AKG <-> GLU + SPYR 0 CYS AT 

SPYR + H20 -> H2S03 + PYR 0 SPTB 

LYS + NADPH + AKG -> NADP + H20 + SAC 0 LKR1 

SAC + H20 + NAD •> GLU + NADH ♦ AASA 0 LKR2 

AASA+ NAD -> NADH + AADP 0 2ASD 

AADP + AKG -> GLU + KADP 0LOC5 

TRP + 02 -> FKYN 0TDO2 

FKYN ♦ H20 -> FOR + KYN 0 KYNF 

KYN + NADPH + 02 ->HKYN + NADP + H20 0KMO 

HKYN + H20 -> HAN + ALA 0 KYNU2 
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HAN + 02 -> CMUSA 
CMUSA->C02 + AM6SA 
AM6SA->PIC 

AM6SA + NAD -> AMUCO + NADH 

AMUCO + NADPH -> KADP + NADP + NH4 

ARG -> ORN + UREA 

ORN + Hm -> ORNm 

ORN + Hm + CTTRm <-> CfTR + ORNm 

ORNm + CAPm-> OTRm + Pirn + Hm 

CUR + ASP + ATP <-> AMP + PPi + ARGSUCC 

ARGSUCC->FUM + ARG 

PRO + FAD->P5C + FADH2 

P5C + NADPH -> PRO + NADP 

THR -> NH3 + H20 + OBUT 

THR + NAD -> C02 + NADH + AMA 

AMA + H20 + FAD •> NH3 + FADH2 + MTHGXL 

GLYm + THFm+ NADm <-> METTHFm + NADHm + C02m + NH3m 

PHE + THBP + 02 -> TYR + DHBP ♦ H20 

NADPH + DHBP -> NADP + THBP 

AKG + TYR -> HPHPYR + GLU 

HPHPYR + 02 ■> HGTS + C02 

HGTS + 02 -> MACA 

MACA •> FACA 

FACA + H20 -> FUM + ACA 

AKG + flJE-> OMVAL + GLU 

OMVALm + COAm + NADm -> MBCOAm + NADHm + C02m 

MBCOAm + FADm -> MCCOAm + FADH2m 

MCCOAm + H20m -> MHVCOAm 

MHVCOAm + NADm -> MAACOAm + NADHm 

MAACOAm -> ACCOAm + PROPCOAm 

2 ACCOAm <-> COAm + AACCOAm 

AKG + VAL->01VAL + GLU 

OrVALm + COAm + NADm -> IBCOAm + NADHm + C02m 
IBCOAm + FADm -> MACOAm + FADH2m 
MACOAm + H20m -> HIBCOAm 
H IBCOAm + H20m -> HIBm + COAm 
HIBm + NADm -> MMAm + NADHm 

MMAm + COAm + NADm •> NADHm + C02m + PROPCOAm 
PROPCOAm + C02m + ATPm -> ADPm + Plm + DMMCOAm 
DMMCOAm •> LMMCOAm 
LMMCOAm -> SUCCOAm 
AKG + LEU -> OICAP + GLU 

OICAPm + COAm + NADm •> IVCOAm + NADHm + C02m 

OICAPm + COAm + NADH *> IVCOAm + NADHm + C02m 

OICAPm + COAm + NADHm -> IVCOAm + NADHm + C02m 

IVCOAm + FADm -» MCRCOAm + FADH2m 

MCRCOAm + ATPm + C02m + H20m -> MGCOAm + ADPm + Pirn 

MGCOAm + H20m -> H3MCOAm 

H3MCOAm -> ACCOAm + A,CTACm 

MYO ACT + ATP •> MYOATP + ACTIN 

MYOATP + ACTIN -> MYOADPAC 

MYOADPAC -> ADP + PI + MYOACT + CONTRACT 

PCRE + ADP -> CRE + ATP 

AMP + H20->PI + ADN 

ATP + AMP <-> 2 ADP 

02<->02m 

3HB -> 3HBm 

CfT + MALm <-> CTTm + MAL 

PYR <^> PYRm + Hm 

C160CAR + COAm -> C160COAm + CAR 

OMVAL •> OMVALm 

OIVAL-> OIVALm 

OICAP -> OICAPm 

GL<->GLm 

GL3Pm + FADm -> T3P2m + FADH2m 

T3P2 + NADH <-> GL3P «• NAD 

GL3P <-> GL3Pm 

T3P2 <-> T3P2m 

OAm + GUJm <-> ASPm + AKGm 

OA + GLU <-> ASP + AKG 

AKG + MALm <-> AKGm + MAL 

ASPm + GLU + H -> Hm + GLUm + ASP 

GLCxt->GLC 

02x1 -> 02 

C160Axt + FABP -> C160FP + ALBxt 

C150FP -> C160 + FABP 

C180Axt + FABP •> C180FP + ALBxt 

C180FP->C180 + FABP 

C16lAxt + FABP ■> C161FP + ALBxt 

C161FP -> C161 + FABP 

C161Axt + FABP -> C181FP + ALBxt 



DHAAO 
OACSD 
OSPTA 
OAMSD 
0 2AMR 
0ARG2 
OORNMT 

ooRNcrrr 

OOTC 

OASS 

0 ASL 

OPRODH 

0PYCR1 

OWTDH 

OTDH 

OMAQA 

DAMT 

OPAH 

OQDPR 

OTAT 

OHPD 

OHGD 

0GST21 

OFAH 

0BCAT1A 

0 BCKDHAA 

OACADMA 

0ECHS1B 

0 EHHADHA 

0ACAA2 

OACATml 

0BCAT1B 

0 BCKDHAB 

OACADSB 

0 EHHADHC 

0 HIBCHA 

OEHHADHB 

0 MMSDH 

OPCCA 

0 HIBCHF 

OMUT 

0BCAT1C 

OBCKDHAC 

0 BCKDHBC 

ODBTC 

OIVD 

0MCCC1 

OHIBCHB 

OHMGCL 

OMYOSA 

OMYOSB 

OMYOSC 

OCREATA 

OCREATB 

OCREATC 

0O2MT 

OHBMT 

OCITMC 

OPYRMC 

0 C160CM 

0 HIBCHC 

0 HBCHD 

0 HIBCHE 

OGLMT 

0GPD2 

0GPD1 

0GL3PMC 

0T3P2MC 

0GOT1 

0GOT2 

OMALMC 

OASPMC 

0GLUT4 

0O2UP 

0FAT1 

0FAT2 

0FAT3 

0FAT4 

OFATC 

OFATB 

0FAT7 
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C181FP->C181 + FABP 

C182Axt + FABP -> C182FP + ALBxt 

C182FP->C182 + FABP 

C204Axt + FABP -> C204FP + ALBxt 

C204FP->C204 + FABP' 

PYRxt + HEXT<->PYR+H 

LACxt + HEXT <-> LAC + HEXT 

H<-> HEXT 

C02 <-> C02m 

H20<->H20m 

ATP + AC + COA •> AMP + PPI + ACCOA 

C160CAR <-> C160CARm 

CARm<->CAR 

C02xt<->C02 

H20xt<->H20 

Pbct + HEXT o HEXT + PI 

<-> GLCxt 

<-> PYRxt 

<-> C02xt 

<-> 02xt 

<-> PIxt 

<-> H20xt 

<-> LACxt 

<-> C16GAxt 

<-> C161Axt 

<-> ClBQAxt 

<-> C181Axt 

<-> C182Axt 

<-> C204Axt 

<-> ALBxt 

<-> 3HB 

<-> GLYCOGEN 

<->PCRE 

<-> TAGLYm 

<->ILE 

<->VAL 

<^CRE 

<->ADN 

<->PI 
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OFATB 

0FAT9 

0 FAT10 

0 FAT11 

0 FAT12 

OPYRUP 

OLACUP 

OHextUP 

0CO2MT 

0 H20MT 

0FU2 

0C160MT 

OCARMT 

0CO2UP 

0H2OUP 

OPIUP 

OGLCexR 

0 PYRexR 

0CO2exR 

0O2exR 

OPIexR 

0H2OexR 

OLACexR 

0C160AexR 

0 C161AexR 

OC160AexR 

0 C181AexR 

0 C182AexR 

0 C204AexR 

OALBexR 

OHBexR 

OGLYex 

OPCREex 

OTAGmex 

OtLEex 

OVALex 

OCREex 

DAONex 

OPlex 
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What is claimed is: 

1. A computer readable medium or media, 
comprising: 

(a) a data structure relating a plurality of 
5 Homo sapiens reactants to a plurality of Homo sapiens 

reactions, 

wherein each of said Homo sapiens reactions 
comprises a reactant identified as a substrate of the 
reaction, a reactant identified as a product of the 
10 reaction and a stoichiometric coefficient relating said 
substrate and said product, 

wherein at least one of said Homo sapiens 
reactions is annotated to indicate an associated gene; 

(b) a gene database comprising information 
15 characterizing said associated gene; 

(c) a constraint set for said plurality of 
Homo sapiens reactions, and 

(d) commands for determining at least one 
flux distribution that minimizes or maximizes an 

.20 objective function when said constraint set is applied 
to said data representation, wherein said at least one 
flux distribution is predictive of a Homo sapiens 
physiological function. 

2. The computer readable medium or media of 
25 claim 1, wherein said plurality of Homo sapiens 
reactions comprises at least one reaction from a 
peripheral metabolic pathway. 



3. The computer readable medium or media of 
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claim 2, wherein said peripheral metabolic pathway is 
selected from the group consisting of amino acid 
biosynthesis, amino acid degradation, purine 
biosynthesis, pyrimidine biosynthesis, lipid 
5 biosynthesis, fatty acid metabolism, cof actor 
biosynthesis and transport processes. 

4 . The computer readable medium or media of 
claim 1, wherein said Homo sapiens physiological 
function is selected from the group consisting of 

10 growth, energy production, redox equivalent production, 
biomass production, production of biomass precursors, 
production of a protein, production of an amino acid, 
production of a purine, production of a pyrimidine, 
production of a lipid, production, of a fatty acid, 

15 production of a cofactor, transport of a metabolite, 

and consumption of carbon, nitrogen, sulfur, phosphate, 
hydrogen or oxygen. 

5. The computer readable medium or media of 
claim 1, wherein said Homo sapiens physiological 

20 function is selected from the group consisting of 

degradation of a protein, degradation of an amino acid, 
degradation of a purine, degradation of a pyrimidine, 
degradation of a lipid, degradation of a fatty acid and 
degradation of a cof actor. 

25 6. The computer readable medium or media of 

claim 1, wherein said data structure comprises a set of 
linear algebraic equations. 

7. The computer readable medium or media of 
claim 1, wherein said data structure comprises a 
30 matrix. 
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The computer readable medium or media of 



claim 1, wherein said commands comprise an optimization 
problem. 

9. The computer readable medium or media of 
5 claim 1, wherein said commands comprise a linear 

program. 

10. The computer readable medium or media of 
claim 1, wherein at least one reactant in said 
plurality of Homo sapiens reactants or at least one 

10 reaction in said plurality of Homo sapiens reactions is 
annotated with an assignment to { a subsystem or 
compartment . 

11. The computer readable medium or media of 
claim 10, wherein a first substrate or product in said 

15 plurality of Homo sapiens reactions is assigned to a 
first compartment and a second substrate or product in 
said plurality of Homo sapiens reactions is assigned to 
a second compartment. 

20 12. The computer readable medium or media of 

claim 1, wherein a plurality of said Homo sapiens 
reactions is annotated. to indicate a plurality of 
associated genes and wherein said gene database 
comprises information characterizing said plurality of 

25 associated genes. 

13. A computer readable medium or media, 
comprising: 

(a) a data structure relating a plurality of 
Homo sapiens reactants to a plurality of Homo sapiens 
30 reactions, 
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wherein each of said Homo sapiens reactions 
comprises a reactant identified as a substrate of the 
reaction, a reactant identified as a product of the 
reaction and a stoichiometric coefficient relating said 
5 substrate and said product, 

wherein at least one of said Homo sapiens 
reactions is a regulated reaction; 

(b) a constraint set for said plurality of 
Homo sapiens reactions, wherein said constraint set 

10 includes a variable constraint for said regulated 
reaction, and 

(c) commands for determining at least one 
flux distribution that minimizes or maximizes an 
objective function when said constraint set is applied 

15 to said data representation, wherein said at least one 
flux distribution is predictive of a Homo sapiens 
physiological function. 

14. The computer readable medium or media of 
claim 13, wherein said variable constraint is dependent 
20 upon the outcome of at least one reaction in said data 
structure. 

•15. The computer readable medium or media of 
claim 13, wherein said variable constraint is dependent 
upon the outcome of a regulatory event. 

25 16. The computer readable medium or media of 

claim 13, wherein said variable constraint is dependent 
upon time. 



17. The computer readable medium or media of 
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claim 13 r wherein said variable constraint is dependent 
upon the presence of a biochemical reaction network 
participant. 

5 18. The computer readable medium or media of 

claim 17, wherein said participant is selected from the 
group consisting- of a substrate, product, reaction, 
protein, macromolecule, enzyme and gene. 

IQ 19. The computer readable medium or media of 

claim 13, wherein a plurality of said reactions are 
regulated reactions and said constraints for said 
regulated reactions comprise variable constraints. 

20. A computer readable medium or media, 

15 comprising: 

(a) a data structure relating a plurality of 
Homo sapiens skeletal muscle cell reactants to a 
plurality of Homo sapiens skeletal muscle cell 
reactions, wherein each of said Homo sapiens 

20 reactions comprises a reactant identified as a 

substrate of the reaction, a reactant identified as a 
product of the reaction and a stoichiometric 
coefficient relating said substrate and said product; 

(b) a constraint set for said plurality of 
25 Homo sapiens reactions, and 

(c) commands for determining at least one 
flux distribution that minimizes or maximizes an 
objective function when said constraint set is applied 
to said data representation, wherein said at least one 

30 flux distribution is predictive of Homo sapiens 
skeletal muscle cell energy production. 



21. A method for predicting a Homo sapiens 
physiological' function, comprising: 
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(a) providing a data structure relating a 
plurality of Homo sapiens reactants to a plurality of 
Homo sapiens reactions , 

wherein each of said Homo sapiens reactions 
5 comprises a reactant identified as a substrate of the 
reaction, a reactant identified as a product of the 
reaction and a stoichiometric coefficient relating said 
substrate and said product , 

wherein at least one of said Homo sapiens 
10 reactions is annotated- to indicate an associated gene; 

(b) providing a constraint set for said 
plurality of Homo sapiens reactions; 

(c) providing an objective function, and 

(d) determining at least one flux 

15 distribution that minimizes or maximizes said objective 
function when said constraint set is applied to said 
data structure, thereby predicting a Homo sapiens . 
physiological function related to said gene. 

22. The method of claim 21, wherein said 

20 plurality of Homo sapiens reactions comprises at least 
one reaction from a peripheral metabolic pathway. 

23. The method of claim 22, wherein said 
peripheral metabolic pathway is selected from the group 
consisting of amino acid biosynthesis, amino acid 

25 degradation, purine biosynthesis, pyrimidine 
biosynthesis, lipid biosynthesis, fatty acid 
metabolism, cofactor biosynthesis and transport 
processes . 



24. The method of claim 21, wherein said 
30 Homo sapiens physiological function is selected from 
the group consisting of growth, energy production, 
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redox equivalent production, biomass production, 
production of biomass precursors, production of a 
protein, production of an amino acid, production of a 
purine, production of a pyrimidine, production of a 
5 lipid, production of a fatty acid, production of a 

cofactor, transport of a metabolite, and consumption of 
carbon, nitrogen, sulfur, phosphate, hydrogen or 
oxygen. 

25. The method of claim 21, wherein said 
Homo sapiens physiological function is selected from 
the group consisting of glycolysis, the TCA cycle, 
pentose phosphate pathway, respiration, biosynthesis of 
an amino acid, degradation of an amino acid, 
biosynthesis of a purine, biosynthesis of a pyrimidine, 
biosynthesis of a lipid, metabolism of a fatty acid, 
biosynthesis of a cofactor, transport of a metabolite 
and metabolism of a carbon source, nitrogen source, 
oxygen source, phosphate source, hydrogen source or 
sulfur source. 

20 26. The method of claim 21, wherein said 

data 

structure comprises a set of linear algebraic 
equations . 

27. The method of claim 21, wherein said 
25 data 

structure comprises a matrix. 

28. The method of claim 21, wherein said 
flux 

distribution is determined by linear programming. 

30 



10 



15 
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29. The method of claim 21, further 
comprising: 

(e) providing a modified data structure, 
wherein said modified data structure comprises at least 

5 one added reaction, compared to the data structure of 
part (a) , and 

(f) determining at least one flux 
distribution 

that minimizes or maximizes said objective function 
10 when said constraint set is applied to said modified 
data structure, thereby predicting a Homo sapiens 
phy s iological function . 

30. The method of claim 29, further 
comprising 

15 identifying at least one participant in said at least 
one added reaction. 

31. The method of claim 30, wherein said 
identifying at least one participant comprises 
associating a Homo sapiens protein with said at least 

20 one reaction. 

32. The method of claim 31, further 
comprising 

identifying at least ^one gene "that encodes said 
protein. 

25 33. The method of claim 30, further 

comprising 

identifying at least one compound that alters the 
activity or amount of said at least one participant, 
thereby identifying a candidate drug or agent that 
30 alters a Homo sapiens physiological function. 
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34. The method of claim 21, further 

comprising: 

•(e) providing a modified data structure, 
wherein said modified data structure lacks at least one 
5 reaction compared to the data structure of part (a), 
and 

(f) determining at least one flux 
distribution 

that minimizes or maximizes said objective function 
10 when said constraint set is applied to said modified 
data structure, thereby predicting a Homo sapiens 
physiological function . 

35. The method of claim 34, further 
comprising 

15 identifying at least one participant in said at least 
one reaction. 

36. The method of claim 35, wherein said 
identifying at least one participant comprises 
associating a Homo sapiens protein with said at least 

20 one reaction. 

37. The method of claim 36, further 
comprising 

identifying at least one gene that encodes said protein 
that performs said at least one reaction. 

25 38. The method of claim 35, further 

comprising 

identifying at least one compound that alters the 
activity or amount of said at least one participant, 
thereby identifying a candidate drug or agent that 
30 alters a Homo sapiens physiological function. 
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39. The method of claim 21, further 

comprising: 

(e) providing a modified constraint set, 
wherein said modified constraint set comprises a 

5 changed constraint for at least one reaction compared 
to the constraint. for said at least one reaction in the 
data structure of part (a) , and 

(f) determining at least one flux 
distribution 

10 that minimizes or maximizes said objective function 
when said modified constraint set is applied to said 
data structure, thereby predicting a Homo sapiens 
physiological function. 



identifying at least one participant in said at least 
one reaction . 

41. The method of claim 40 , wherein said 
identifying at least one participant comprises 

20 associating a Homo sapiens protein with said at least 
one reaction. 

42. The method of claim 41, further 
comprising 

identifying at least one gene that encodes said 
25 protein. 



identifying at least one compound that alters the 
activity or amount of said at least one participant, 
30 thereby identifying a candidate drug or agent that 
alters a Homo sapiens physiological function. 



15 



40. 



The method of claim 39, further 
comprising 



43. The method of claim 40, further 
comprising 
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44. The method of claim 21 , further 
comprising 

providing a gene database relating one or more 
reactions in said data structure with one or more genes 
5 or proteins in Homo sapiens. 

45. A method for predicting a Homo sapiens 
physiological function, comprising: 

(a) providing a data structure relating a 
plurality of Homo sapiens reactants to a plurality of 
10 Homo sapiens reactions, 

wherein each of said Homo sapiens reactions 
comprises a reactant identified as a substrate of the 
reaction, a reactant identified as a product of the 
reaction and a stoichiometric coefficient relating said 
15 substrate and said product, 

wherein at least one of said Homo sapiens 
reactions is a regulated reaction; 

(b) providing a constraint set for said 
plurality of Homo sapiens reactions, wherein said 

20 constraint set includes a variable constraint for said 
regulated reaction; 

(c) providing a condition-dependent value to 

said variable constraint; 

(d) providing an objective function, and 
25 (e) determining .at least one flux 

distribution that minimizes or maximizes said objective 
function when said constraint set is applied to said 
data structure, thereby predicting a Homo sapiens 
physiological function . 



30 



46. The method of claim 45, wherein said 
value 
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provided to said variable constraint changes in 
response to the outcome of at least one reaction in 
said data structure. 

47. The method of claim 45, wherein said 
5 value 

provided to said variable constraint changes in 
response to the outcome of a regulatory event. 

48. The method of claim 45, wherein said 
value 

10 provided to said variable constraint changes in 
response to time. 

49. The method of claim 45, wherein said 
value 

provided to said variable constraint changes in 
15 response to the presence of a biochemical reaction 
network participant. 

50. The method of claim 4 9, wherein said 
participant is selected from the group consisting of a 
substrate, product, reaction, enzyme, protein, 

20 macromolecule and gene. 

51. The method of claim 45, wherein a 
plurality of said reactions are regulated reactions and 
said constraints for said regulated reactions comprise 
variable constraints. 

25 52. A method for predicting Homo sapiens 

growth, comprising: 

(a) providing a data structure relating a 
plurality of Homo sapiens skeletal muscle cell 
reactants to a plurality of Homo sapiens skeletal 

30 muscle cell reactions, 
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wherein each of said Homo sapiens reactions 
comprises a reactant identified as a substrate of the 
reaction, a reactant identified as a product of the 
reaction and a stoichiometric coefficient relating said 
*5 substrate and said product; 

(b) providing a constraint set for said 
plurality of Homo sapiens reactions; 

(c) providing an objective function, and 

(d) determining at least one flux 

10 distribution that minimizes or maximizes said objective 
function when said constraint set is applied to said 
data structure, thereby predicting Homo sapiens 
skeletal muscle cell energy production. 

53. A method for making a data structure 
15 relating a plurality of Homo sapiens reactants to a 
plurality of Homo sapiens reactions in a computer 
readable medium or media, comprising: 

(a) identifying a plurality of Homo sapiens 
reactions and a plurality of Homo sapiens reactants 

20 that are substrates and products of said Homo sapiens 
reactions; 

(b) relating- said plurality of Homo sapiens 
reactants to said plurality of Homo sapiens reactions 
in a data structure, 

25 wherein each of said Homo sapiens reactions 

comprises a reactant identified as a substrate of the 
reaction, a reactant identified as a product of the 
reaction and a stoichiometric coefficient relating said 
substrate and said product; 

30 (c) determining a constraint set for said 

plurality of Homo sapiens reactions; 

(d) providing an objective function; 

(e) determining at least one flux 
distribution that minimizes or maximizes said objective 
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function when said constraint set is applied to said 
data structure, and 

(f) if said at least one flux distribution 
is not predictive of a Homo sapiens physiological 
5 function, then adding a reaction to or deleting a 
reaction from said data structure and repeating step 
(e), 

if said at least one flux distribution is 
predictive of a Homo sapiens physiological function/ 
10 then storing said data structure in a computer readable 
medium or media. 



54. The method of claim 53, wherein a 
reaction 

in said data structure is identified from an annotated 
15 genome. 

55. The method of claim 54, further 
comprising 

storing said reaction that is identified from an 
20 annotated genome in a gene database. 

i 

56. The method of claim 53, further 
comprising 

annotating a reaction in said data structure. 

57. The method of claim 56, wherein said 
annotation is selected from the group consisting of 
assignment of a gene, assignment of a protein, 
assignment of a subsystem, assignment of a confidence 
rating, reference to genome annotation information and 
reference to a publication. 

58. The method of claim 53, wherein step (b) 



25 



30 
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further comprises identifying an unbalanced reaction in 
said data structure and adding a reaction to said data 
structure, thereby changing said unbalanced reaction to 
a balanced reaction. 

59. The method of claim 53, wherein said 
adding a reaction comprises adding a reaction selected 
from the group consisting of an intra-system reaction, 
an exchange reaction, a reaction from a peripheral 
metabolic pathway, reaction from a central metabolic 
pathway, a gene associated reaction and a non-gene 
associated reaction. . 



60. The method of claim 59, wherein said 
peripheral metabolic pathway is selected from the group 
consisting of amino acid biosynthesis, amino acid 
15 degradation'; purine biosynthesis, pyrimidine 
biosynthesis, lipid biosynthesis, fatty acid 
metabolism, cofactor biosynthesis and transport 
processes. 



20 



61. The method of claim 53, wherein said 
Homo sapiens physiological function is selected from 
the group consisting of growth, energy production, 
redox equivalent production, biomass production, 
production of biomass precursors, production of a 
protein, production of an amino acid, production of a 
25 purine, production of a pyrimidine, production of a 
lipid, production of a fatty acid, production of a 
cofactor, transport of a metabolite, development, 
intercellular signaling, and consumption of carbon 
nitrogen, sulfur, phosphate, hydrogen or oxygen. 



30 



62. The method of claim 53, wherein said 
Homo sapiens physiological function is selected from 
the group consisting of degradation of a protein, 
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degradation of an amino acid, degradation of a purine, 
degradation of a pyrimidine, degradation of a lipid, 
degradation of a fatty acid and degradation of a 
cof actor. 

5 63. The method of claim 53, wherein said 

data 

structure comprises a set of linear algebraic 
equations. 

64. The method of claim 53, wherein said 
10 data 

structure comprises a matrix. 



65. The method of claim 53, wherein said 
flux 

distribution is determined by linear programming. 



15 
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30 



66. A data structure relating a plurality of 
Homo sapiens reactants to a plurality of Homo sapiens 
reactions, wherein said data structure is produced by a 
process comprising: 

(a) identifying a plurality of Homo sapiens 
reactions and a plurality of Homo sapiens reactants 
that are substrates and products of said Homo sapiens 
reactions; 

(b) relating said plurality of Homo sapiens 
reactants to said plurality of Homo sapiens reactions 
in a data structure, 

wherein each of said Homo sapiens reactions 
comprises a reactant identified as a substrate of the 
reaction, a reactant identified as a product of the 
reaction and a stoichiometric coefficient relating said 
substrate and said product; 
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(c) determining a constraint set for said 
plurality of Homo sapiens reactions; 

(d) providing an objective function; 

(e) determining at least one flux 

5 distribution that minimizes or maximizes said objective 
function when said constraint set is applied to said 
. data structure, and 

(f) if said at least one flux distribution 
is not predictive of Homo sapiens physiology, then 

10 adding a reaction to or deleting a reaction from said 
data structure and repeating step (e) , 

if said at least one flux distribution is 
predictive of Homo sapiens physiology, then storing 
said data structure in a computer readable medium or 

15 media. 
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